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Abstract The EFqFi-ATP synthase mutants bQ64C and 
yri06C were labelled selectively with the fluorophores tetrame- 
thylrhodamine (TMR) at the b-subunit and with a cyanine (Cy5) 
at the T^subimit. After reconstitution into Uposomes, these dou- 
ble-labelled enzymes catalyzed ATP synthesis at a rate of 33 
s~*. Fluorescence of TMR and Cy5 was measured with a con- 
focal set-up for single-molecule detection. Photon bursts were 
detected, when Uposomes containing one enzyme traversed the 
confocal volume. Three states with different fluorescence reso- 
nance energy transfer (FRET) efficiencies were observed. In the 
presence of ATP, repeating sequences of those three FRET- 
states were identified, indicating stepwise rotation of the T^sub- 
unit of EFqFi. © 2002 Federation of European Biochemical 
Societies. Published by Elsevier Science B.V. All rights re- 
served. 

Key words: H+-ATP synthase; EFqFi; 
Inter-subumt rotation; Single-molecule fluorescence 
resonance energy transfer 



1. Introduction 

Membrane-bound H+-ATP synthases catalyze the forma- 
tion of ATP from ADP and inorganic phosphate in mitochon- 
dria, chloroplasts and bacteria. Endergonic ATP synthesis is 
coupled to proton translocation across the membrane. The 
enzyme consists of two major parts. The hydrophobic mem- 
brane-integrated Fo part is involved in proton transport 
across the membrane, whereas the hydrophilic Fi part con- 
tains nucleotide and phosphate binding sites. In Escherichia 
coli, the Fi part consists of five different subunits with the 
stoichiometry aaPsTSe and the Fo part is built of three differ- 
ent subunits with a likely stoichiometry ab2Cio-i4 [!]• The 
three catalytic nucleotide binding sites on the ^-subunits 
undergo conformational changes [2,3], adopting in sequential 
order the 'open', 'tight' smd 'loose' conformation during ATP- 
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synthesis [4]. These changes are induced by 'docking-undock- 
ing* steps of the y-subunit to the three different ap-pairs, i.e. 
by a rotation of the y-subunit (for reviews, see [4—6]). 

ATP-induced rotation of the y-subimit within the asps bar- 
rel has been deduced from ensemble measurements [7,8] and 
was convincingly visualized by a single-enzyme approach [9]. 
In these experiments, the as Pay-fragment of the H+-ATPase 
from Bacillus PS3 was attached with three His-tags onto the 
glass surface of a microscopic coverslip, and, as an indicator, 
a fluorescent pm-long actin filament was connected to the 
T^subunit. The orientation of the actin filament with respect 
to the asPs hexamer was monitored by videomicroscopy. At 
very low ATP concentrations, rotation occurred in discrete 
120*^ steps consistent with sequential ATP hydrolysis at the 
three P-subunits [10]. Recently, substeps during ATP hydro- 
lysis have been identified [1 1]. Rotation of the y-subunit was 
shown by the fluorescent actin filament method with Fi parts 
from chloroplasts and E. coli [6] and was used to demonstrate 
rotation of the e-subunit and the c-ring in FqFi [12,13]. Ac- 
cording to these results, subunits TECio-u ('rotor') rotate 
against subunits a3P35ab2 ('stator') counter-clockwise during 
ATP hydrolysis, when viewed from the membrane to Fi. 
However, it was not possible to demonstrate with this method 
that proton transport is coupled with the movement of the 
rotor subunits. 

To ensure coupling during ATP hydrolysis and to maintain 
full enzymatic activity for ATP synthesis, we incorporated the 
holoenzyme FoFi into a Kpid membrane and labelled specif- 
ically the y-subunit and b-subunit with two different fluoro- 
phores chosen for efficient fluorescence resonance energy 
transfer (FRET). These indicator molecules were small 
enough to allow for an unconstrained subimit movement. In- 
tramolecular subunit movement resulting in a change in the 
distances between the fluorophores were detected by changes 
in FRET efficiencies. After reconstitution of the labelled FqFi 
into liposomes, ATP synthesis was observed by generating a 
pH difference and an additional electric potential difference 
across the membrane. Using confocal single-molecule fluores- 
cence detection, we measured subunit movements during ATP 
hydrolysis by intramolecular FRET in fully functional single 
EFoFi reconstituted in liposomes. 

2. Materials and methods 

The plasmid pRA114 carrying the yri06C mutation was con- 
structed by Aggeler and Capaldi [14] and was expressed in strain 
RAl [15]. EFi was isolated according to [16]. The y-subunit of EFi 
was selectively labeUed with carbocyanine dye Cy5®-maleimide in 50 
mM MOPS/HCI (pH 7.0) as described previously [17,18]. Cy5-malei- 
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mide was kindly provided by E. Schweinberger and C.A.M. Seidel 
(MPI fur Biophysikalische Chemie, Grottingen, Germany). EFi con- 
centrations were determined by UV absorption using an extinction 
coeSiaent e= 191 000 cm~^ at 280 nm. For Cy5, an extinction 
coeffident of 250 000 an~^ at 649 nm was used. Unbound dye 
was removed after 4 min by passing twice through Sephadex GSO 
centrifuge columns. Spedfid^ of y-subunit labelling was diecked by 
fluorescence detection after SDS-PAGE. Approximately 65% of EFi 
was labelled with Cy5. EFi-yri06C-Cy5 was stored at -80**. 

In order to introduce a cysteine in the b-subunit, the new plasmid 
pRR76A)C21S/bQ64C was constructed from pRAlOO [19] with 
pNOC/bC21S obtained from R. FiUingame [20]. The mutation 
bQ64C was introduced into pNOC by a two -stage PCR mutagenesis 
procedure [2i] to yield pRR36. The fragment with the two mutations 
was ligated to pRAlOO. The resulting plasmid pRR76 was expressed 
in strain RAl [15]. EFqFi was isolated according to [22]. The b-sub- 
units in this EFqFi were labelled with tetramethylrhodamine (TMR)- 
maleimide (Molecular Probes) as described above. The degree of la- 
belling was adjusted to approximately 20% to avoid labelling of both 
b-subunits at the same enzyme. EFqFi -bQ64C-TMR was reconsti- 
tuted in liposomes [22] with a stoichiometry of one EFqFi per lipo- 
some. The Fi parts were removed [23] and the Eposomes containing 
EFo-bQ64C-TMR were reassembled with EFi-yri06C-Cy5 to yield 
Kposomes with the double-labeUed enzyme EFo-bQ64C-TMR-Fi- 
yri06C-Cy5 according to [24,25]. The efficiencies of removal of Fi 
and reassembling with labelled Fi were controlled by measuring ATP 
synthesis activities. The initial rate of ATP synthesis catalyzed by 
EFqFi was measured after an acidr-base transition at 23°C as de- 
scribed previously [22]. 

Single-molecule fluorescence measurements of EFqFi in liposomes 
were performed with a confocal set-up of local design [25]. A fre- 
quency doubled Nd:YAG laser (532 nm, 50 mW, Coherent, Ger- 
many) was used for excitation. The laser beam was attenuated to 
120 |iW and focussed by a water immersion objective (UAPO 40 X, 
N.A. 1.15, Olympus) into the buffer solution placed on a microscope 
slide. Fluorescence was measured in two spectral regions. Single pho- 
tons were detected with avalanche photodiodes (SPCM-AQR-15, 
EG&G, Canada) after passing an interference filter HQ 575nm/ 
65nm for TMR and a HQ 665nm LP for Cy5 (AHF, Germany). 
Photons in two channels were counted in parallel by a multidiannel 
scaler (PM S 300, Becker & Hickl, Berlin, Germany). The actual con- 
focal detection volume V=1J fl was calculated from fluorescence 
correlation spectroscopy using rhodamine 6G according to [18]. AH 
fluorescence measurements were performed in buflJer A (50 mM 
HEPES/NaOH, 2.5 mM MgCb, 400 jiM sodium ascorbate, pH 
8.0). Fluorescent impurities of the buffer were removed by activated 
charcoal granula (Merck, Germany) with subsequent sterile filtration. 
ATP and adenosine- 5 '-(p,Y-imido)triphosphate (AMPPNP) were ob- 
tained from Bohringer (Mannheim, Germany). The concentration of 
proteolipo somes was adjusted such that only one liposome was de- 
tected in the confocal volume at any time. 



3. Results 

EFqFi was labelled with two fluorophores. In order to ob- 
serve inter-subunit movement in functionally intact EFqFi by 
intramolecular FRET, it is necessary to label spedficaUy one 
subunit of the *stator' part and one of the 'rotor' part. For 
labelling of the b-subunit with the FRET donor TMR, a 
cysteine in the b-subunit at position 64 was introduced. For 
labelling of the y-subunit with the FRET acceptor Cy5, we 
used the mutant EFi-'yri06C [141. Selective labelling was 
achieved by the following procedure. 

Step 1 : Isolated EFqFi -bQ64C was labelled with TMR at a 
substoichiometric concentration. A ratio of 0.2 TMR per 
EFqFi was chosen to avoid labelling of both b-subunits. 
EFqFi -bQ64C-TMR was reconstituted into liposomes with 
approximately one EFqFi per liposome. The rate of ATP syn- 
thesis was measured after energization by ApH/A^ with a rate 
of 61 s~^. This rate is comparable to that of wild-type EFqFi 
[22]. Therefore, labelling of the b-subunit with TMR did not 
change ATP synthesis activity (see Table 1). Step 2: Fi parts 
of this enzymes were removed and liposomes with incorpo- 
rated EFQ-bQ64C-TMR were retained. EFq -Kposomes showed 
no ATP synthesis activity, neither with TMR-labelled EFq nor 
with unlabeUed EFq, indicating an efficient removal of the Fi 
parts. Step 3: EFi-yri06C was isolated and labelled with Cy5. 
Step 4: Finally, we incubated the EFo-Hposomes with the 
EFi-yri06C-Cy5 (and in control experiments with unlabeUed 
EFi-yri06C). After removal of the non-assembled Fi parts by 
ultracentrifugation, we obtained EFqFi in liposomes, in which 
the b-subunit in the EFq part was specifically labelled with 
TMR and the y-subunit in the Fi part was specifically labelled 
with Cy5. The ATP synthesis rate of this double-labelled 
EFqFi was 33 s~^, which corresponds to 50% of the rate 
before removal of Fi. A similar result was obtained when 
unlabeUed EFi was reassembled with EFq. Addition of 
DCCD blocked almost completely ATP synthesis (Table i). 

Using double-labeUed enzymes reconstituted in liposomes 
CFRET-labeUed EFqFi', EFQ-bQ64C-TMR-Fi-yri06C-Cy5), 
intramolecular FRET between TMR at the b-subunit and Cy5 
at the y-subunit was measured by a confocal set-up with two- 
channel detection. Photons were detected in the spectral range 
of 545-610 nm from TMR (*green channel') and at wave- 



Table 1 

Rates of ATP synthesis catalyzed by labeUed EFqFi in liposomes 



Conditions 



Rate (s-^ 



b-mutant 

EFoFi-bQ64C 61 ±3 

b-mutant after removal of Fi 

EFo-bQ64C 0 
labelled b-mutant 

EFo-bQ64C-TMR-Fi 61 ± 4 

labeUed b-mutant after Temoval of Fi 

EFo-bQ64C-TMR 0 
labeUed b-mutant after removal of Fi and reassembled with Fi 

EFo-bQ64C-TMR-Fi 29 ± 4 

b-mutant after removal of Fi and reassembled with labeUed Fi 

EFo-bQ64C-Fi-yri06C-Cy5 32 ± 6 

labeUed b-mutant after removal of Fi and reassembled with labeUed Fi 

EFo-bQ64C-TMR-Fi-yri06C-Cy5 33 ± 3 

EFo-bQ64C-TMR-Fi-yri06C-Cy5 in presence of lO^iDCCD 4±3 

EFqFi with mutations bQ64C and yri06C was reconstituted iato Uposomes and the rates of ATP synthesis were measured after energization 
by ApH and A0. Given are the mean and standard deviations firom six independent measurements. 
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Fig. 1. Photon bursts and proximity factor traces of single FRET-labelled EFqFi in liposomes. Single-molecule fluorescence of FRET-labeHed 
EFqFi reconstituted in liposomes is observed as photon bursts, when liposomes containing one FRET-labeUed enzyme traverse the confocal de- 
tection volume. Green: time traces of TMR fluorescence intensities; red: simultaneous time trace of Cy5. The blue traces in the panels above 
are the calculated proximity factors. In the presence of AMPPNP, three states with different FRET efliciencies (a-c) are observed in which the 
proximity factor remains constant during the burst, d: Reconstituted EFo-bQ64C-TMR-Fi represent enzymes with FRET donor only. e,f: In 
the presence of ATP, fluctuations of the proximity factor corresponding to the three FRET states are detected within one burst. The sequence 
of states is indicated above the proximity factor trace. 



lengths above 665 mn from Cy5 (*red channer). The average 
background signal from buffer A containing 1 mM AMPPNP 
or 1 mM ATP was approximately 3 photons/ms in both chan- 
nels. In the presence of FRET-labelled EFoFi-Hposomes in 
solution, well-separated photon bursts with count rates up 
to 150 photons/ms were observed. This indicates a single lipo- 
some migrating through the detection volume. At the begin- 



ning of the time trace a in Fig. I, both channels detect only 
background signal. At 255 ms, a liposome enters the detection 
volume and the TMR fluorescence intensity (green trace) in- 
creases to approximately 150 counts/ms. At 370 ms, the pro- 
teoliposome leaves the detection volume and and fluorescence 
decreases to background level. Simultaneously, an increase of 
Cy5 fluorescence with a count rate of 35 counts/ms is detected 
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Fig. 2. Histogram of proximity factors from single FRET-labelled 
EFoFi in the presence of 1 mM AMPPNP. L: low FRET effidency 
state; M: mediiim FRET efficiency state; H: high FRET efficiency 
state; DO: 'donor only' state from TMR-labeled EFqFi without 
FRET acceptor. The normalized distributions are calculated from 
1 ms time intervals in 10 long-lastiag photon bursts for each FRET 
state. 



(red trace). Since Cy5 is not excited directly at the laser wave- 
length and power applied, this indicates inter-subunit FRET 
from TMR to Cy5. 

There are two reasons for the fluctuations in fluorescence 
intensities. First, the proteoliposome moves through the focus 
driven by Brownian motion. Excitation power and fluores- 
cence that is detected depend on the position of the enzyme 
within the confocal volume. Second, the efficiency of FRET 
changes with inter-subunit movement. FRET efficiency was 
calculated as the proximity factor P [26] . This factor is the 
ratio of the acceptor fluorescence intensity (/a) divided by the 
sum of donor and acceptor intensities: 



P = /a/(/d + /a) 



(1) 



The blue trace in Fig, la shows the proximity factor for the 
event. Although there were large changes in the fluorescence 
intensities in both channels within the photon burst, the prox- 
imity factor remained nearly constant around P = 0.2. In the 
following, this state is called L ('low FRET state'). 

We calculated the proximity factors for different photon 
bursts and classified them into four FRET states: L (*low') 
with PwO.2, M ('medium') with P«0.5, H ('high') with 
P»0.8 and DO ('donor only') with i><0.1. Examples for 
the different states are shown in Fig, ia. ci. In FRET state 
M, fluorescence intensities of donor and acceptor are almost 
equal; in FRET state H, acceptor count rate is higher than 
donor count rate. In state DO, the acceptor fluorescence is 
near background. As a rare event, we also observed photo- 
bleaching of the acceptor fluorophore within the burst (data 
not shown) indicated by a spontaneous decrease of the accep- 
tor count rate to background level and accompanied by a 
simultaneous rise of the donor count rate. Photobleaching 



of Cy5 resulted in state DO. The same fluorescence intensity 
ratio for the DO state with P < 0. 1 was measured in control 
experiments using reconstituted EFo-bQ64C-TMR-Fi, which 
contained only the donor fluorophore on the b-subunit. The 
DO state was observed in the FRET-labelled ATP synthase 
preparation, because only 65% of EFi contained Cy5 and, as 
a consequence, 35% of the reconstituted EFqFi were lacking 
the acceptor label. Detection of some photons in the acceptor 
channel ('cross talk') is therefore due to the spectral properties 
of TMR. This cross talk leads to an * apparent proximity fac- 
tor' which does not reflect FRET. The states L, M, and H 
were attributed to diffierent conformations of the H+-ATP 
synthase. These are characterized according to the Forster 
theory of FRET [27] by different distances between the donor 
at the b-subunit and the acceptor at the y-subunit. Therefore, 
FRET state L reflects a long distance, M a medium distance 
and H a short distance. 

When single-enzyme fluorescence was measured in the pres- 
ence of 1 mM ATP, well-separated photon bursts were ob- 
served. However, the time course within one burst was quite 
different. Fig. 1 e/ show two examples during ATP hydrolysis. 
In contrast to the bursts in the presence of AMPPNP, there 
was a change in the green-to-red intensity ratio, which can be 
seen more clearly using the calculated proximity factor (blue 
trace). The proximity factors fluctuated within one photon 
burst. Three levels of the proximity factor could be distin- 
guished (L, M and H), which were comparable in magnitude 
to the three different states observed in the presence of 
AMPPNP. During ATP hydrolysis, a definite order of prox- 
imity factor transitions was observed: from the low FRET (L) 
to the medium FRET (M) to the high FRET (H) and then 
again to the low FRET (L) state. In a few cases, we observed 
within this sequence L^M->H^L also 'backward steps' 
(see Fig„ If, the high FRET state starting at 560 ms seems 
to oscillate between H and M states before reaching the L 
state at 700 ms) and in some photon bursts, the sequence 
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Fig. 3. Distribution of dweU times for the FRET states in presence 
of ATP. DweU times of 40 FRET states were measured from pho- 
ton bursts of FRET-labelled EFqFi io. liposomes, showing at least 
three distinct FRET states. First and last states of these bursts axe 
omitted for the histogram, i.e. the states during entering and leaving 
the detection volume. The monoexponential decay fit yields a mean 
lifetime of 30 ms. 
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Fig. 4. Model of the FRET-label positions in EFoFi. Left: Side view of EFqFi. The homology model by S. Engelbrecht (http://13i.l73,26,96/ 
sc/se.html) was aligned with the structure of the y- and e-subunit [30]. Amino acid position of the FRET donor TMR at the b-subunit is indi- 
cated with 'D'; amino acid position of the FRET acceptor Cy5 at the y-subunit is indicated with 'A'. Right: Cross section at the fluorophore 
level, viewed from the Fq to the Fi part. Position bQ64C is labelled 'D'; the distance between bQ64C and Yri06C is labelled L (low FRET). 
Rotation of the y-subunit by 120° steps results in two additional positions labelled M (medium FRET) and H (high FRET). The sequence 
L-»M-»H-»L corresponds to counter-clockwise rotation of the y-subunit during ATP hydrolysis. 



was reversed (L -» H -> M -» L). In all cases the movements 
occurred stepwise within our experimental time resolution 
for defining a FRET state (approximately 5 ms). 

To calculate a mean value for the proximity factor, we 
chose 10 photon bursts for each FRET state from the experi- 
ments with AMPPNP, in which the sum of donor and accep- 
tor intensities was more than 20 counts/ms. The number of 
events as a function of the proximity factors is plotted in Fig. 
2. The mean values of these distributions are P= (0.18 ±0.1) 
for L, P = (0.47 ±0.1) for M and P = (0.8 ±0.1) for H. For 
DO, an apparent value of P = (0.1 ±0.05) was obtained. 

From the experiments during ATP hydrolysis, we selected 
40 long-lasting photon bursts and measured the duration of 
the different FRET states. These dwell times were sorted into 
10 ms intervals. The distribution (Fig. 3) was fitted by a 
monoexponential decay function which yielded a lifetime of 
30 ms. This means that the average Kfetime of the FRET 
states is 30 ms. If one ATP molecule is hydrolyzed during 
one FRET state, this corresponds to a rate of ATP hydrolysis 
of 33 s~^, which agrees with the rate measured under compa- 
rable conditions in bulk experiments [28]. 

4. Discnssion 

H+-ATP synthase from E. coli carrying TMR at position 64 
of the b-subunit and Cy5 at position 106 of the y-subunit 
catalyzed ATP synthesis, when the proteoliposomes were 



energized (ApH»4 plus A0). Thus, the attachment of the 
two fluorophores did not disturb the conformational dynam- 
ics of the enzyme during proton transport-coupled ATP syn- 
thesis. By a single-molecule detection method, we measured 
the fluorescence of both fluorophores simultaneously, when 
liposomes containing one FRET labelled enzyme traversed 
tiie confocal volume. The intramolecular FRET efficiencies 
were determined by calculating proximity factor P. Three dif- 
ferent FRET states characterized by proximity factors of 
around 0.2, 0.5, and 0.8 were observed, corresponding to three 
different distances of Cy5 at the y-subunit relative to TMR at 
the b-subunit. In the presence of AMPPNP, these distances 
remained constant within the photon bursts, i.e. the enzymes 
were 'arrested' in one of the three different states. In the 
presence of ATP, the same three states were observed as char- 
acterized by their proximity factors. However, during ATP 
hydrolysis fluctuations between these three states occurred 
within one burst. The average lifetime of each FRET state 
was 30 ms. Since the time resolution is limited by the integra- 
tion time of 1 ms used in these experiments, we are only able 
to estimate the upper limit for the transition time between the 
different states, which is less than 3 ms. 

During ATP hydrolysis the distance between the fluoro- 
phores at the y-subunit and the b-subunit remains constant 
within one conformational state, until a fast transition to the 
next state occurs. In this new state, the enzyme remains again 
for a time interval followed by a fast transition to the third 
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state. From this third state, the next transition leads to the 
first conformation again. Therefore, the distances between the 
y-subunit and the b-subunit change in a cyclic manner. Based 
on the current knowledge of the structure of FqFi, we inter- 
pret these data in a molecular model (Fig. 4) which is based 
on homology modelling and alignment of structural informa- 
tion. This is shown in a side view in Fig. 4a and as a cross 
section at the fluorophore level (viewed from Fq to the Fi 
part). With b-subunits located at the interface of one ap- 
pair, the distance between bQ64C and Yri06C is 7.3 nm in 
the model. When we start at this position (L) and rotate the 
Y"Subunit in steps of 120* coimter-clockwise around an axis 
going through the center of Fi, the distances between bQ64C 
and Yn06C are approximately 5 nm (M) and 2 nm (H). These 
distances closely match the values calculated from the prox- 
imity factors, if we assume Rq — 5.3 nm [29] according to the 
Forster theory of FRET [27]. 

In some photon bursts, we observed a clockwise order of 
transitions. We attribute these events to enzymes in which the 
second b-subunit is labelled with TMR (position *D2' in Fig. 
4b). However, in most experiments, the sequence of FRET 
states was L^M^H^L. This implies that the sequence of 
states observed in the single-molecule FRET experiments is in 
accordance with a counter-clockwise rotation of the -y^subunit 
relative to the b-subunit during ATP hydrolysis. We conclude 
from this result that the y-subunit rotates stepwise relative to 
the b-subunit in fully active EFqFi in liposomes. 
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Abstract Using a confocal fluorescence microscope with an 
avalanche photodiode as detector, we studied the fluorescence of 
the tetramethyh'hodaniine labeled Fi part of the IT^-ATPase 
from Escherichia coH^ £Fi, carrying the 7X106-0 mutation 
[Aggeler, J.A. and Capaldi, RA. (1992) J. Biol. Chem. 267, 
21355-21359] in aqueous solution upon excitation with a mode- 
locked argon ion laser at 528 nm. The diffusion of the labeled 
EFi through the confocal volume gives rise to photon bursts, 
which were analyzed with fluorescence correlation spectroscopy, 
resulting in a diffusion coefGcient of 3.3 X 10"*^ cm^ s~^. In the 
presence of nucleotides the diffusion coefficient increases by 
about 15%. This effect indicates a change of the shape and/or the 
volume of the enzyme upon binding of nucleotides, i.e. 
fluorescence correlation spectroscopy with single £Fi molecules 
idlo^ the detection of coltformational chang J! 

1998 Federation of European Biochemical Societies. 

Key words: H+-ATPase; Conformational change; 
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1. Introduction 

FoFi-ATP synthases occur in bacteria, mitochondria and 
chloroplasts [1]. These enzymes catalyze ATP synthesis 
coupled with a transmembrane proton flux. They have two 
large domaias: a membrane iategrated Fq part, which is in- 
volved in proton transport, and a hydrophilic Fi part, which 
contains the nucleotide binding sites. In Escherichia coli, the 
Fi part contains five different subunits with the stoichiometry 
otsPsTSe. Based on functional studies, it was suggested that the 
catalytic nucleotide binding sites on the p-subunits operate in 
a cycHc way, which is accomplished by subsequent 'docking- 
undocking' steps of the T^subunit to the three ap pairs. The 
nucleotide binding and dissociation steps at the Fi part are 
coupled in the holoenzyme via long range conformational 
changes with proton binding, translocation and dissociation 
steps in the Fq part (for review see [2]). In this work we 
investigated conformational changes of the Fi part of the 
ATP S3mthase from E. coli, EFi, which occur upon binding 
and dissociation of nucleotides. We labeled a cysteine in the 

♦Corresponding author. Fax: (49) (761) 203-6189. 
E-mail: boeisch@ruf.iini-freibiu-g.de 

Abbreviations: EFi, hydrophilic Fi part of the proton translocatiag 
H+-ATPase from Escherichia coli; AMPPNP, P,Y-in[iidoadenosine-5'- 
triphosphate ; TMRIA, tetramethylrhodamine-5-iodoacetamide; 
MOPS, 3-[iV-morpholLQo] propanesulfonic acid; HEPES, iV-[2-hy- 
droxyethyl] piperaziae iV-[2-ethanesulfonic acid]; TMR, tetramethyl- 
rhodamine; MCS, multichannel scaler; FCS, fluorescence correlation 
spectroscopy 



y-subunit of EFi with tetramethylrhodamine and measured 
the fluorescence of single molecules by confocal laser spectros- 
copy [3]. From an analysis of the fluorescence bursts using 
fluorescence correlation spectroscopy (FCS), the characteristic 
time of diffusion of the enzyme through the confocal volume 
can be derived. We found that the diffusion time changes 
upon nucleotide binding and conclude that this reflects struc- 
tural changes in EFi. A preliminary report of this work has 
been presented before [4]. 

2. Materials and methods 

Straia pRA114/AN888, carryiag the Yri06-C mutation, is described 
in Aggeler and Capaldi [5] and was a gift from these authors. EFi was 
isolated from the mutant strain as described [6] except that the ion- 
exchange chromatography was performed on a HQ 20 FPLC column 
(PerSeptive). Isolated EFi was stored in liquid nitrogen until use. 

The y-subiuiit of EFi was selectively labded with substoichiometric 
amounts of tetiamethylrhodamiae-S-iodoacetaniide (TMRIA, Molec- 
ular Probes) m 50 mM MOPS/HCl (pH 7.0) with 10% glycerol as 
described [7\. EFi concentrations were determined by UV absorption 
using an extinction coefficient of 190730 M~^ cm~^ at 280 nm as 
described [8]. Dye concentrations were determmed usmg an extmction 
coefficient of 87000 M~^ cm~^ at 543 nm. Unbound dye was removed 
after 4 min reaction time by passing twice through Sephadex G50 
centrifuge columns. The labeling of the y-subunit was specific as de- 
termmed by the fluorescence of the gel after sodium dodecylsulfate 
polyacrylamide gel electrophoresis. The TMR labeled enzymes, TMR- 
EFi, were stored at -80°C. 

Fluorescence measurements were done in 50 mM HEPES/NaOH 
(pH 8.0) and 2.5 mM MgCl2 (buffer B) after removing fluorescent 
impurities by activated charcoal granular (1.5 mm, Merck). TMR-EFi 
aliquots were diluted to final concentrations of 10~^° M. ATP and P,y- 
imidoadenosine- 5 '-triphosphate (AMPPNP) stock solutions (10 mM) 
were prepared fresh in buffer B and were mixed with the diluted 
TMR-EFi solutions immediately before the fluorescence measure- 
ments. TMR-glutathione was prepared by mixing TMRIA with a 
100-fold excess of the tripeptide glutathione (y-Glu-Cys-Gly, Sigma) 
and diluted to a final dye concentration of 7 X 10~^^ M for single 
molecule spectroscopy. Confocal fluorescence detection was per- 
formed with an active mode locked argon ion laser (excitation wave- 
length 528 nm, repetition rate 73 MHz) (Sabre, Coherent, Palo Alto, 
CA) using an epi-illuminated confocal fluorescence microscope (Fig. 
1) similar to that described previously [9,10] with a beam-splitter at 
530 nm, a 150 {im pinhole and a dichroic band-pass emission filter 
582/50 nm (AF Analysentechmk, Tubingen, Germany). The e^anded 
laser beam was attenuated to 530 |xW and focused by a water immer- 
sion objective (UPLANAPO 60 X, NA= 1.2, Olympus) to an excita- 
tion volume of a few femtohters. The fluorescence signal was detected 
by a single photon counting avalanche photodiode (AQ 151, EG&G, 
VaudreuU, Quebec, Canada). The multiplexed detector signal was 
registered in parallel by a multichannel scaler (MCS) using a PC- 
adapter counter (CIO-CTR05, Plug-in GmbH, Eichenau, Germany) 
and by a real-time correlator card (ALV-5000/E, ALV, Langen, Ger- 
many) for FCS. The laser focus and detection volume parameters 
(radial and axial 1/e^ radius o — 0.65 |xm and zq = 3.3 \im, respec- 
tively) were determined from FCS measurements of rhodamme 6G 



0014-5793/98/$! 9.00 © 1998 Federation of European Biochemical Societies. All rights reserved. 
PII: S0014-5793(98)01247-2 



252 



Af. Barsch et allFEBS Letters 437 (1998) 251-254 



in pure water with a characteristic difiusion time of = 0.38 ms (see 
Table 1), assuming a translational diffusion coefficient of 



D = 2.8x10 ^ cm^ s ^ in water [11]. This corresponds to a detection 
volume of 



i2 s-i 

y- 1.5 



zq=1.1 fl and a focal area of A = 12 
o^ = 6.7xlO"9 cm^ [12,13]. Thus, the excitation power of P = 530 
P-W is equivalent to a quasi-continuous focal irradiance of Iq — PI 
>4 = 8 X 10^ WW (see Fig. 1). 

The samples were transferred to the microscope support by a micro- 
scope slide with a small depression and covered with a conventional 
cover glass. AH measurements were performed at room temperature. 



3. Results and discussioii 

Upon free diffusion in the solution, the TMR-labeled EFi 
enters the confocal detection volume, where the dye is excited 
by the pulsed laser light and due to the high laser repetition 
rate cycles continuously between the ground and excited states 
(see Fig. 1). Fluorescence photons are emitted as bursts when 
a smgle enzyme diffuses through the confocal volume. At the 
mean excitation irradiance of /o/2 = 4XlO^ W/cm^, the mean 
number of photons emitted by a siagle TMR fluorophore in 
water is about 10^ before photobleaching occurs [14]. Since 
the transit time of the labeled EFi is about 3 ms, the molecule 
can emit about 30 000 photons during this period [14]. Assum- 
ing a detection efficiency of about 3% [15], about 900 fluores- 
cence photons should be detected from a single molecule 
transit. Furthermore, the arrangement of the avalanche pho- 
todiode with a 150 \\m pinhole allows the selective detection 
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Fig. 1. Scheme of experimental setup for siagle molecule fluores- 
cence detection. The pulsed laser beam (wavelength 528 nm) is fo- 
cused through the microscope objective onto the sample. The fluo- 
rescence of the sample molecules entering the confocal detection 
volimie is focused on a 150 ^un pinhole and detected by a single 
photon coimting avalanche photodiode. 
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Fig. 2. Multichannel scaler trace of fluorescence signals from EFi 
labeled with tetramethyh-hodamine. TMR-EFi was diluted in 50 
mM HEPES/HCl (pH 8), 2.5 mM MgCl2 to a concentration of 
10~^° M. Data acquisition was performed at a speed of 1000 points 
per second (integration time 1 ms). The inset shows an expanded 
view of two fluorescence bursts. 



of the fluorescence from the confocal volume with a strong 
suppression of background fluorescence from other parts of 
the solution (see Fig. 1). 

Fig. 2 shows a MCS trace of TMR-EFi (10"^° M in buffer 
B) to monitor the fluorescence signals (bursts) arising from 
single TMR-EFi molecules diflusing through the confocal de- 
tection volume. The single bars of the MCS trace represent 
the number of fluorescence photons counted by the avalanche 
photodiode during an integration time of 1 ms. The trace 
shows burst rates of up to 220 kHz, while a background level 
of approximately 3 kHz is observed. This corresponds to a 
signal-to-noise ratio of 70. The background signal mainly re- 
sults from scattered Hght and partially from background fluo- 
rescence of molecules outside tiie confocal detection volume. 
Due to the spatially dependent laser excitation irradiance and 
fluorescence collection efficiency of the open detection volume, 
the randomly diffusing molecules can either cross near the 
edges or traverse the center of the volume. Therefore, an in- 
homogeneous distribution of the number of detected fluores- 
cence photons from a single molecule transit (burst size) is 
obtaiaed [15]. This can be seen ia Fig. 2, where events of 
different duration and with different count rates can be 
seen. The inset in Fig. 2 shows a magnification of a part of 
the MCS trace, showing a molecule diffusing through the 
center of the confocal detection volume (burst at 125 ms), 
and one diffusing through the edge (burst at 110 ms), where 
the excitation irradiance and, thus, the count rate is smaller. 
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Results of the fluoTescence correlation spectroscopy data analysis 


Compound 




T 


tji (ms) 


tA (ms) 


A 




Rhodamine 6G 


2.1 


0.13 


0.38 


none 


none 


5.10 


TMR-glutathione 


2.9 


0.06 


0.40 


none 


none 


0.26 


TMR-EFi 


5.7 


0.10 


3.2 


0.16 


0.03 


0.60 


TMR-EFi + ATP 


5.0 


0.11 


2.6 


0.11 


0.05 


0.37 


TMR-EFi + AMPPNP 


4.4 


0.10 


2.8 


0.20 


0.10 


0.37 



Parameteis of a fit of the PCS data to Eq. 1 : characteristic triplet correlation time, ; average fraction of molecules in the excited triplet state, T\ 
characteristic time for translational diffusion, tjy, average number of fluorescing molecules in the confocal detection volume, Ny\ characteristic 
correlation time of the 'dynamic quenched' state, tx fraction of molecules in the 'dynamic quenched* state, A. The optical parameter ( q/zq) 
was determined by fitting the PCS data of rhodamine 6G and fixed for all subsequent fits. All other parameters were allowed to vary freely. The 
estimated error for tj} is 7%. 



PCS [12,13] was used to provide precise statistical data on 
the number of molecules in the detection volume and on the 
average diffusion time of the molecules through the detection 
volume. Analysis of the PCS data was performed with a least- 
squares fit algorithm using the following expression of the 
normalized autocorrelation function, G(tc), with the correla- 
tion time tc. 



G(t^) = 1 + 



1 ^ 1 ^/ 



1 



1/2 



I+^cAd 1 + ( o/zo)\/tB 



X(l-r+rexp (-tJtT)A+A exp (-^cAa)) 



(1) 



Nf is the average number of fluorescent molecules in the 
Gaussian detection volume with the 1/e^ radii, o and zq, in 
radial and axial direction, respectively, tn = o^/4Z>, the char- 
acteristic time for translational diffusion of the observed mol- 
ecules with the diffusion coefficient, D, through the detection 
volume, A the average fraction of molecules within a so-caUed 
dynamic quenched state, tA the characteristic correlation time 
of this quenched state and T the average fraction of molecules 
in the excited triplet state with a characteristic triplet correla- 
tion time, tT [11]. 

Pig. 3 shows the measured autocorrelation curves of TMR- 



EPi in the absence of nucleotides and in the presence of 
AMPPNP together with the weighted residuals [16] of a fit 
of the data to Eq. 1, indicating a very precise description of 
the data by Eq. 1. The optical parameter ( q/zq) was deter- 
mined by fitting the rhodamine 6G PCS data and fixed for all 
subsequent fits. All other parameters were allowed to vary 
freely. The PCS data of pure rhodamine 6G and TMR-gluta- 
thione (data not shown) could be described properly by an 
autocorrelation function, G{tc), just using the diffusion and 
triplet terms of Eq. 1, i.e. the four parameters, iVp, tjy, T 
and tT. In the case of TMR-EPi the addition of the so-caUed 
dynamic term, i.e. of the parameters A and tA, is necessary. 
Obviously, this additional term, with parameters A and tA, in 
Eq. 1 is only required for an optimal fit of the TMR-EPi data 
and, therefore, it reflects a property of the fluorophore in the 
protein environment. Table 1 shows the parameters of the 
proper fits of the PCS data from measurements of different 
samples to G(tc). 

Prom the concentrations, c, of the different solutions and 
the size of the detection volume, V= 7.7 fl (see Section 2), an 
average number of molecules, N= cXV, in this volume can be 
calculated: TMR-EPi solutions: N=0.5, TMR-glutathione 
solution: JV=0.32, rhodamine 6G solution: N=4.6. These 
estimations of N are in good agreement with the average 




time / ms 



Fig. 3. Normalized fluorescence autocorrelation curves of TMR-EFi. Solid line: TMR-EFi without added nucleotides. Dashed line: TMR-EFi 
in the presence of ImM AMPPNP normalized to the curve without nucleotides. The weighted residuals [16] from the fit of the data to Eq. 1 
are shown at the bottom. Fit parameters are sunmiarized in Table 1. 
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number of fluorescent molecules, Ny, obtained from FCS (Ta- 
ble 1). 

From Table 1 and Fig. 3 it becomes obvious that the addi- 
tion of nucleotides (ATP or AMPPNP) to TMR-EFi results 
in a significant change of the characteristic difiusion time to. 
After addition of 1 mM ATP, tjy decreased by a factor of 
0.81, while addition of 1 mM AMPPNP resulted in a decrease 
of to by a factor of 0.87. The same effects could also be 
observed in additional experiments where an enlarged confo- 
cal detection volume was employed. Using the relation 
to = o^/4Z>, between the known radial 1/e^ radius, o = 0-65 
^m, and the characteristic translational diffusion time, to, the 
translational diffusion coefficient, D, of TMR-EFi can be cal- 
culated: Z)(TMR-EFi) = 3.3 X 10"^ cm^ s-\ addition of 1 mM 
ATP: Z)(TMR-EFi, ATP) = 4.1 X 10"^ cm^ s'S and addition 
of 1 mM AMPPNP: i)(TMR-EFi, AMPPNP) = 3.8 X 10"^ 
cm2 s-^ It is obvious tiiat the addition of tiie nucleotides 
increases the diffusion coefficient of TMR-EFi, and that tiiis 
effect is more enhanced with ATP. 

The diffusion coefficient, D, is related to the frictional co- 
efficient, /, via D = kT/f, where k is the Boltzmann constant 
and T the absolute temperature. At constant temperature and 
viscosity, the frictional coefficient, /, depends on the radius 
and the shape of the molecule [17]. We therefore conclude 
that binding of nucleotides (either ATP or AMPPNP) leads 
to a large conformational change in the Fi part of the en- 
zyme. 

Similar large conformational changes in the Fi part of H+- 
ATPases were observed with completely different methods. 
Treatment of mitochondrial MFi witii glycerol resulted in a 
release of tiie bound nucleotides and in a large decrease of the 
sedimentation coefficient, s, (from 11.9 S to 8.4 S) [18]. This 
large change was reversed when glycerol was removed. 

Direct evidence for a large change of shape upon nucleotide 
binding comes from X-ray analysis. The X-ray structure of 
mitochondrial Fi with tiie subunits asPs^ee containing tiiree 
AMPPNP, one ADP, and one ATP per enzyme has been 
determined [19]. Recently, the X-ray structure of the subcom- 
plex asps from the bacillus PS 3 containing no bound nucleo- 
tides has been determined as well [20]. Although the molecular 
masses of tiie subunits a and p of PS 3 are similar to that of 
mitochondrial Fi and the asPs complex does not contain the 
7-, 5-, and e-subunits, it is much larger than tiie mitochondrial 
Fi. The main difference between the two structures seems to 
be that all three nucleotide binding pockets on the p-subunits 
are empty and in the open conformation in asPs, whereas in 
MFi only one p-subunit is in such a conformation. Presum- 
ably, the large difference between the two structures resulted 
at least partly from tiie different nucleotide occupancy of botii 
enzymes. Supporting this conclusion, small-angle X-ray scat- 
tering data indicated a structural change of Fi from PS 3 
upon ADP binding which was modelled as a shrinkage of 
the six major subunits by 6% along their major axis [21]. 

Best fits to the FCS data of TMR labeled EFi-ATPases 
were achieved only after introduction of a second bunching 
term (so-called d3niamic quenched term). This term character- 
izes an equilibrium in the microsecond time range, which in- 
fluences the fluorescence of the fluorophore. Most interest- 
ingly, addition of ATP or AMPPNP affected the 'dynamic' 
time, tA, of TMR-EFi in opposite directions. With ATP the 
'dynamic' time is shortened, indicating an acceleration of tiiis 



process with respect to TMR-EFi without bound ATP. In 
contrast, with bound AMPPNP the 'dynamic' time is pro- 
longed. Due to the small amplitudes of this effect, the quanti- 
fication of tiiese processes is currentiy difficult. However, it 
might be of special interest since AMPPNP binds tightiy to 
the enzyme but, in contrast to ATP, is not hydrolyzed. There- 
fore, this 'dynamic' process might be associated with move- 
ments of the y-subunit during catalysis. 

The data in Table 1 show longer triplet correlation times 
and larger amplitudes of TMR-EFi compared to TMR-gluta- 
thione. This might be an indication tiiat diffusion of oxygen to 
quench tiie triplet state of TMR is sterically hindered when 
the fluorophore is attached to the protein. These data indicate 
that single molecule fluorescence - either combined with pol- 
arized excitation and emission [22] or combined with fluores- 
cence correlation spectroscopy analysis — can provide new 
information on conformational changes of the Fi part upon 
nucleotide binding. 

Acknowledgements: We thank Roderick A. Capaldi aad Robert Ag- 
geler for generous support and gift of the mutants makiag this work 
possible. 

References 



[1] Pfedersen, P.C. and Carafoli, E. (1987) Trends Biochem. Sd. 12, 
146-150. 

[2] Boyer, P.D. (1997) Annu. Rev. Biochem. 66, 717-749. 

[3] Keller, R.A., Ambrose, W.P., Goodwin, P.M., Jett, J.H., Martin, 

J.C. and Wu, M. (1996) Appl. Spectrosc. 50, 12A-32A. 
[4] Borsch, M., Turina, P., Eggeling, C, Fries, J.R., Seidel, C.A.M., 

Labahn, A. and Graber, P. (1998) Ital. Biochem. Soc. Trans. 11. 
[5] Aggeler, J. A. and Capaldi, R.A. (1992) J. Biol. Chem. 267, 

21355-21359. 

[6] Gogol, E.J., Luecken, U., Bork, T. and Capaldi, R.A. (1989) 

Biochemistry 28, 4709-4716. 
[7] Turina, P. and Capaldi, R.A. (1994) J. Biol. Chem. 269, 13465- 

13471. 

[8] Gill, S.C. and von Hippel, P.H. (1989) Anal. Biochem. 182, 319- 
326. 

[9] Zander, C, Sauer, M., Drexhage, K.H., Ko, D.-S., Schulz, A., 

Wolfrum, J., Brand, L., Eggeling, C. and Seidel, C.A.M. (1996) 

Appl. Phys. B 63, 517-523. 
[10] Brand, L., EggeHng, C, Zander, C, Drexhage, K.H. and Seidel, 

C.A.M. (1997) J. Phys. Chem. A. 101, 4313-4321. 
[11] Widengren, J., Mets, U. and Rigler, R. (1995) J. Phys. Chem. 99, 

13368-13379. 

[12] Elson, E.L. and Magde, D. (1974) Biopolymers 13, 1-27. 

[13] Thompson, N.L. (1991) in: Topics in Fluorescence Spectroscopy, 
Volume 1: Techniques (Lakowicz, J.R., Ed.), pp. 337—378, Ple- 
num Press, New York. 

[14] Eggeling, C, Widengren, J., Rigler, R. and Seidel, C.A.M. (1998) 
Anal. Chem. 70, 2651-2659. 

[15] Fries, J.R., Brand, L., Eggehng, C, KoUner, M. and Seidel, 
C.A.M. (1998) J. Phys. Chem. A (in press). 

[16] Koppel, D.E. (1974) Phys. Rev. A 10, 1938-1946. 

[17] Tanford, C. (1961) Physical Chemistry of Macromolecules, Wi- 
ley, Chichester. 

[18] Garret, N.E. and Penefsky, H.S. (1975) J. Supramol. Struct. 3, 
469-^78. 

[19] Abrahams, J.P., Leslie, A.G.W., Lutter, R. and Walker, J.E. 

(1994) Nature 370, 621^28. 
[20] Shirakiliara, Y., Leslie, A.G.W., Abrahams, J.P., Walker, J.E., 

Ueda, T., Sekimoto, Y., Kambara, M., Saika, K., Kagawa, Y. 

and Yoshida, M. (1997) Structure 5, 825-^36. 
[21] Furuno, T., Ikegami, A., Kihara, H., Yoshida, M. and Kagawa, 

Y. (1983) J. Mol. Biol. 170, 137-153. 
[22] Hasler, K., Engelbrecht, S. and Junge, W. (1998) FEBS Lett. 426, 

301-304. 



articles 

Resolution of distinct rotational 
substeps by submiliiseconci kinetic 
analysis of Fi -ATPase 

Ryohel Yasuda^tll, HlroyuU No|l^, Masasuke Yoshldal:^, Kazuhlko Klnoslta Jrt^ & HIroyasu ltDh$^ 

* CREST 'Genetic Programming Team 13, Teikyo University Biotechnology Center 3F, Nogawa 907, Miyamae-I^ Kawasaki 216-0001, Japan 
f Department of Physics, Faculty of Science and Technology, Keio University, Yokohama 223-8522, Japan 
t Chemical Resources Laboratory, Tokyo Institute of Technology, Yokohama 226-8503, Japan 
§ Tsukuba Research Laboratory, Hamamatsu Photonics KK, Tokodai, Tsukuba 300-2635, Japan 



The enzyme Fi-ATPase has been shown to be a rotary motor in which the central 7-subiinit rotates inside the cylinder made of 
03 Ps subunits. At low ATP concentrations, the motor rototes in discrete 120° steps, consistent with sequential ATP hydrolysis on 
the three |3-subunite. The mechanism of stepping is unlcnown. Here we show by high-speed imaging that the 120° step consiste of 
roughly 90° and 30° substeps, each teking only a fraction of a millisecond. ATP binding drives the 90° substep, and the 30° substep 
is probably driven by release of a hydrolysis product. The two substeps are separated by two reactions of about 1 ms, which 
together occupy most of the ATP hydrolysis cycle. This scheme probably applies to rotetion at full speed (-^130 revolutions per 
second at saturating ATP) down to occasional stepping at nanomolar ATP concentrations, and supports the binding-change model 
for ATP synthesis by reverse rotetion of Fi -ATPase. 



The ATP synthase is an enzyme ubiquitous in bacteria, plants and 
animals, which synthesizes ATP from ADP and inorganic phosphate 
using proton flow through a membrane^ Fi, a water-soluble 
portion of the ATP synthase, is the site of ATP synthesis, whereas 
protons flow through the membrane-embedded Fq portion. At least 
in vitrOy Fj can hydrolyse ATP to pump protons through the Fq 
portion in the reverse direction. Isolated Fi only hydrolyses ATP, and 
is called Fj -ATPase. Its subunit composition is asPs^ySe, 

The prevailing view is that ATP hydrolysis/synthesis in Fi is 
coupled to proton flow in Fq through the rotation of a common 
shaft, of which the 7-subunit of Fi is a part. This rotational coupling 
mechanism was initially proposed by Boyer^"^, and by others'*"*. 
Later, a crystal structure of Fi showed that a rod-shaped 'y-subunit is 
surrounded by a cylinder made of three a- and three P- subunits, 
arranged alternately^ (Fig. la). An analogue of ATP, ADP and none 
were bound to the three p -subunits, indicating that sequential ATP 
hydrolysis on the three p -subunits would indeed induce rotation of 
the central, asymmetrical 7- subunit. Rotation of the 7 -subunit in an 
isolated Fi during ATP hydrolysis has been demonstrated experi- 
mentally by various methods^"^°. 

We have visualized the rotation of the -y-subunit under an optical 
microscope by fixmg Fj on a surfece and attaching an actin filament 
to the 7-subunit as a marker of its orientation^". At nanomolar ATP, 
the actin filament rotated in discrete 120° steps", consistent with the 
pseudo -three- fold symmetrical structure^ of Fj. The rotation rate 
was close to one-third of the rate of ATP hydrolysis in solution, 
suggesting that one ATP molecule is consumed per 120° step". At 
high ATP concentrations, however, the actin rotation was smooth 
rather than stepwise, and the rotation was much slower than ATP 
hydrolysis. Viscous friction imposed on the actin filament prevented 
fast rotation of Fi and obscvired the stepping behaviour. Here we 
have used a smaller marker, a colloidal gold bead of 40-nm diameter, 
for which the viscous fiiction is 10"^ to 10"^ times that for actin 
(Fig. lb). In the resultant high-speed rotation, we were able to 
resolve substeps. In this study we investigate the magnitudes. 
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speeds, and timings of the substeps, and we also look at: (1) the 
motor speeds at no load; (2) whether the motor uses different rotary 
mechanisms at low and high speeds; (3) which parts of hydrolysis 
reactions drive the substeps; and (4) what structural changes may 
underlie the substeps. 

Full-speed rotation with 40-nm beads 

Bead rotation was imaged by laser dark-field microscopy^^ (Fig. Ic), 




Figure 1 Observation of F-i rotation, a, Atomic structure^ of Fi-ATPase viewed from the Fq 
side 0op in b). b, Side view of tiie observation system. ITie 40-nm bead gave a iarge 
enougii opticai signai tiiat warranted a submiiiisecond resoiution; but tiie bead was smaii 
enougii not to impede the rotation, c, i^ser dark-fieid microscopy for observation of goid 
beads. Only light scattered by the beads exited the objective and was detected. DFC, darlc- 
field condenser, d, Sequential images of a rotating bead at 2 mM ATP. Images are 
trimmed in circles (diameter 370 nm) to aid identification of the bead position; centroid 
positions are shown above the images at x3 magnification. The interval between images 
is 0.5 ms. 
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and recorded on a fast-framing charge-coupled-device (CCD) 
camera at speeds up to 8,000 frames per s. The 40-nm bead appeared 
as a spot of diffraction-Umited size (~300nm; Fig, Id). When a 
bead is attached obliquely (Fig. lb), rotation of the 'y-subunit will 
result in a circular movement of the bead image. Some beads 
showed rotation (Fig. Id; movies in Supplementary Information), 
and motions of these beads were analysed by calculating the 
centroid of the bead image". The rotation diameter of bead centroid 
ranged between 25-55 nm. Diameters up to ~60nm are possible 
for the height of Fi of ^ 10 nm and the linker lengths of ^5 nm for 
streptavidin^^ and ~10nm for BSA^^ (Fig. lb). Rotation was 
stepwise at aU ATP concentrations examined (see below). 

To see whether the friction on the 40-nm bead impeded Fi 
rotation, we varied the frictional load by attaching single or duplex 
polystyrene beads (108, 196 or 291 nm) to the 'y-subunit. At both 
2 mM and 2 fxM ATP (Fig. 2; red and blue circles, respectively), time- 
averaged rotation rates showed saturation behaviour at small fric- 
tion. Maximal rotation rate depended on ATP concentration [ATP] , 
but increasing [ATP] beyond 2 mM did not accelerate rotation (see 
below). Thus, at 2 mM ATP, the 40-nm bead rotated at the fuU speed 
of the Fi motor, which was 134 revolutions per second (r.p.s.) at 
23 °C; the bead was not an impeding load for Fi. At saturating 
speeds, aU beads rotated stepwise. On the load-dependent portions 
in Fig. 2, however, bead rotation was smooth, as was rotation of 
actin at these ATP concentrations^ \ The load dependence of actin 
rotation^ ^ (Fig. 2; triangles) is consistent with the bead assay. 

One rotary mechanism at all speeds 

The time-averaged rate of rotation showed simple Michaelis- 
Menten dependence on [ATP] (Fig. 3; the maximal rate at infinite 
[ATP], Vmax = 129r.p.s.; Michaelis constant, K^a = 15 juiM), suggest- 
ing that one mechanism accounts for rotation in the nM— mM 
range. This idea is corroborated by observations^^'^*"^' that the 
torque and its angle dependence, as well as mechanical work done 
in a 120° step, are independent of [ATP] over the nM-mM range. 
Also, the apparent rate of ATP binding, k^n^, given by 3 Vmax/^m of 
(2.6 ± 0.5) X 10^ M"^ s"^ agrees with previous estimates based on the 
analysis of step intervals at nanomolar ATP^^'^^. 

As seen in Fig. 3, the rotation rate was dose to one-third of 
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Figure 2 Reiationsiiip between rate of bead rotation and viscous friction on the bead. 
Circles, tfie average rate for a bead calculated over at least 20 consecutive revolutions; 
squares, the average over different beads (vertical envr bars indicate s.d.). ATP at 2 mM 
and 2 |jlM is indicated by red and blue colours, respectively. The abscissa is the rotational 
frictional drag coefficient | calculated as in Metfiods. Possible range of ^ for each bead is 
shown by the size of the squares or associated horizontal error bars. For comparison, 
rotation rates for an actin filament attached to the 7-subunit"''' are also plotted (triangles). 
Lines show fits with the rate expected for a motor producing a constant torque'''': 
(1/l/noioad + 2tt^/\^"'' Where /V = 40 pN nm (assumed torque) and l/noioad = 1 2.5 ± 1 .0 
r.p.s. for 2 jxlW ATP and 1 34 ± 3 r.p.s. for 2 mM ATP (s.e.m.). 



the rate of ATP hydrolysis for bead- free Fj in solution, supporting 
the contention that one ATP molecule is consumed per 120° 
rotation^*'*^*^. The hydrolysis rate, however, was lower, particularly 
around 50|xM. A probable cause is MgADP inhibition: Fj is 
stochastically inactivated during ATP hydrolysis, when it binds 
MgADP tightly^ Although we started with nucleotide-free Fi, 
some inactivation may have proceeded during the mixing time of 
'-'2 s. Indeed, the rate of inactivation increases with [ATP] and 
reaches ~0.3 s~^ at > 10 ^jiM ATP^, the position of the concavity in 
Fig. 3. Higher activity at still higher [ATP] is accounted for by 
binding of ATP to non-catalytic a-subunits, which tends to restore 
the hydrolysis activity^^, Lauryldodecylamine oxide (LDAO), a 
suppressor of the MgADP inhibition^^, produced hydrolysis kinetics 
parallel to the rotation kinetics, although V^nax/d (82 s~^) was only 
-^60% of Vinax for rotation (Fig, 3). 

The 120° step consists of 90° and 30° substeps 

At 8,000 frames per s, steps were dearly resolved in the rotation of 
40-nm beads, even at saturating ATP. At 2 mM ATP, only 120° steps 
were seen (Fig. 4a, b), whereas at 20 |uiM or 2 jxM ATP, the 120° step 
was ftirther split into roughly 90° and 30° substeps (Fig. 4c-f, where 
each panel shows a continuous record). We call the interval between 
a 30° substep and a 90° substep a *0° dwell' and the interval between 
90° and 30° substeps a '90° dweU'. In Fig. 4c-f, 0° dwells faU on black 
horizontal lines that are separated from each other by 120°, and 90° 
dwells faU on grey lines that are 30° below the black lines. The 90° 
dwells were about a few ms in duration, on average, both at 2 and 
20 iJuM ATP, whereas 0° dwells became longer at 2 jxM. The impUca- 
tion is that Fi waits for the arrival of ATP during the 0° dwell, which 
is terminated by a 90° substep induced by ATP binding. 
The subsequent 90° dwell is for a process or processes independent 
of [ATP]. This scheme predicts that, at [ATP] = 15 |ljlM, 0° 

and 90° dwells have approximately equal lengths, as was observed 
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Figure 3 Comparison of rotation and liydrolysis rates. Red circles, time-averaged rotation 
rate for individual 40-nm beads. Red squares, rotation rate averaged over different beads. 
Dark green squares, one-third of the initial rate of ATP hydrolysis. Light green circles, 
one-third of the rate of ATP hydrolysis in the presence of LDAO. Blue diamonds, 
rotation rate for an actin filament attached to the 7-subunit"''' . Standard deviations greater 
than the symbol size are shown in bars {n^2). Curves show fits with Michaells-Menten 
kinetics, V = l/max[ATP]/(A;n + [ATP]), where I/max and ^ are 1 29 ± 9 r.p.s. and 
1 5 ± 2 |jlM for bead rotation (red), 4.0 ± 0.3 r.p.s. and 0.7 ± 0.1 pM for actin rotation 
(biue), and 247 ± 9 s"** and 1 9 ± 1 |xl\1 for hydrolysis in the presence of LDAO (light 
green). Fits with two values, l/= (l/maxi^n2[ATP] + l/max2[ATP]V([ATP]' + K^zU^TP] 
+ ^mi K^), are also shown, where l/maxi = 85 ± 9 s"\ /Clni = 5.2 ± 0.7 |xM, l/max2 = 306 
± 22 s""' , and /(m2 = 393 ± 1 47 |xM for hydrolysis without LDAO (dark green), and 
l/maxi = 1 09 ± 30 r.p.s., /(n,i = 1 2 ± 4 |xM, l/max2 = 1 49 ± 32 r.p.s., and Km2 = 682 ± 
2768 |jlM for bead rotation (dashed pink), "me latter does not show improvement over the 
simple fit in red. Values are means ± s.e.m. 
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Hgure 4 Unfiltered time courses of stepping rotation of 40-nm beads at varying [ATP], 
a, b, 2 mM; c, d, 20 |jlM; e, f, 2 \xM [ATP]. Aii curves in a panei are continuous; iater 
curves are sfiffted, to save space. Grey fiorizontai iines are piaced 30° beiow biacl( iines. In 
e, some of tfie long dweiis are cut sfiort. Insets, positions of a bead wittiin 0.25-0.5 ms 



before (red) and after (green} tfie main (90*" or 1 20°) steps; runs iasting 0.5 s (2 mM) or 2 s 
(2 jjlM and 20 |jlM) were analysed. Circles indicate projectton of —0*" and —90° dwell 
points on an obliquely situated circular trajectory tfiat best fit tfie data. Angles in tfie time 
courses and in Fig. 5 are tfiose on tfie oblique circle. 



at 20 |uiM ATP. At 2 mM ATP, the expected rate of ATP binding is 
(2.6x10^ M"^ s"^ X 2 mJVI) «5xl0^ s"\ Then, 0** dwells will be 
-^0.02 ms and wiU not be detected at the current resolution. The 
absence of substeps in Fig, 4a, b is thus explained. We place dwells in 
Fig. 4a, b on grey lines, because they must be 90° dwells according to 
this explanation (see also Fig. 7). 

The 30° substeps were not always clear, but we could easily locate, 
in rotation records, steps that spanned most of a 120° interval (~90° 
or ---120° step). Positions of the bead centroid in 0.25-0.5-ms 
intervals before (red) and after (green) these main steps are shown 
in the insets of Fig. 4, At 2 |ljiM and 20 |ui.M ATP, red and green spots 
are separated by ^30^ showing the presence of substeps, whereas 
spots overlap with each other at 2 mM ATP. The traces are distorted, 
presumably because of oblique rotation on an obliquely situated Fi. 
Circles on each trace are the projection of three equally spaced pairs 
of dwell positions on a circular trajectory oblique to the surface. A 
search for the best fit with the observed traces show the separation 
between red and green circles to be 29° ± 7° (mean ± s.d.) for 13 
runs at 2, 6 and 20 imM ATP; and 4° ± 3° for 7 runs at 2 and 6 mM. 

Figure 5 shows histograms of angular positions. Separations of 
peaks at 20 |jiM ATP (crosses) averaged 35° ±13° (mean ± s.d. for 
15 peak pairs). Substeps are less obvious at 2 |uiM, but histograms for 
the 2-ms intervals before and after the main steps (green) show 
similar peak separations. Taking these and other experimental 
uncertainties into account, we estimate the substep sizes to be 
within 90° ± 10° and 30° ± 10°. 

Steps are fast 

Consecutive steps in a rotation record at 2mMATP are super- 
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Figure 5 IHistograms of angular positions over 0.5 s runs. Labels (a-f) are from records of 
whldi Fig. 4, a-f are a part. Eacfi time course was passed through a five-point median 
filter, and its histogram was calculated with 2° bins. The histogram was then averaged 
over 1 0° intervals. Green parts (e and f) indicate 2 ms before and after main steps. 
Crosses indicate peaks identified by eye. To assess the baseline noise in raw data, we also 
constructed unfiltered, unaveraged histograms at 2 |juM ATP with 2° bins (not shown). Tlie 
histograms gave three peaks, of which the half width at 1/e height was estimated by fitting 
each peak with a gaussian curve; the half widths averaged 1 8° ± 7° (mean ± s.d. for 1 5 
peaks). 
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imposed in Fig. 6a. The average (thick cyan Hne) shows that a whole 
120'* step completes within 0.25 ms (two frames) at saturating ATP. 
This value is an upper boundary because faster transients are 
unresolved with the camera that we used. Thus, time for the 
mechanical stepping (the times needed to reorientate the y- 
subunit through 90° and 30°) occupies < 10% of the ATPase cycle 
time. The 90° substeps at 2 |xM ATP are also withm 0.25 ms (Fig. 
6b). Stepping is fast, the instantaneous speed being well above 1,000 
r.p.s., whether [ATP] is above or below Kj^ for rotation. 

Substeps of 90° by ATP binding and 30° by product release 

Figure 6c shows averages of aU steps observed at indicated [ATP]. 
Presence of distinct and fast -^90° substeps is clear at all [ATP] < 
60 iulM, although whether the remaining -^30° is also stepwise 
cannot be judged in this figure. Fit with grey lines indicates the 
substep size to be 90.2° ± 0.3°. 

The averaging was made after steps other than the central one 
were eliminated from each step record, such that it started with a 
near —30° dwell and ended with a ^120° dwell. When the last dwell 
at 120° was too short to be distinguished, the previous dwell at ^^90° 
was extended to the right edge of the figure (see Methods). Thus, the 
portion of the curves between -30° and 0° reflects the distribution of 
dwell times at 0° that were started at the end of a substep from -30° 
to 0°, and terminated by a central 90° substep. The dwell is [ATP] - 
dependent, and can be explained by termination by ATP binding at 
the rate k^n^ = 3.0 x 10^ M"^ s~^ estimated from Fig. 8 (Fig. 6c; grey 
Unes). The kinetics above 90°, in contrast, is [ATP] -independent, 
except that the amplitude decreases with [ATP]. The decrease is 
accounted for by the feet that, at [ATP] > K^n, the dwell at - 120° 
becomes too short and is not represented in our averaging pro- 
cedure. As shown (grey lines), the kinetics are in accord with the 
scheme in Fig. 7a where two '^l-ms reactions govern the dwell at 
~90° (the two reactions are deduced from Fig. 8). 

Our proposed scheme is summarized in Fig. 7a. ATP binding 
drives a 90° substep (A-^B). Next are hydrolysis reactions that are 
mechanically silent (B, C). Eventually, the last hydrolysis product of 
the previously bound ATP (ADP, phosphate, or both) is released, 
accompanying a 30° substep (C^A') and resetting the system to the 
initial A state except for the 120° rotation that has taken place. At 
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Figure 6 Kinetics of substeps. a, b, Eighteen consecutive steps and tlieir average (thick 
cyan line) in a rotation record at 2 mM (a) or 2 |jlM (b) ATP (see Methods for the averaging 
procedure), c, Steps in several runs at indicated [ATP], averaged as in a and b. The 
averaging procedure retained all dwells at ~0° position In the ordinate, but some dwells 
at ~ 1 20° were converted to a horizontal line at ~90° position when the substep from 90° 
to 1 20° was contiguous (within 0.25 ms) with the next substep from 1 20° to 21 0°. This Is 
why the curves at high [ATP] do not rise much beyond the 90° line. The curves are 
reproduced on the right; superimposed grey lines represent fits with theoretical curves 
based on the scheme in Fig. 7a (see Methods). The best fit was obtained with the size of 
the 30° substep, Aao-. of 29.8° ± 0.3° (mean ± s.e.m.). 



high [ATP], another 90° substep occurs immediately, and thus the 
two substeps are not resolved (Fig. 7b) — this is why we propose that 
both of the hydrolysis products must have been released by the end 
of a 30° substep. At [ATP] ~ K^, ATP binding and hydrolysis take 
similar time, and 90° and 30° substeps are equally spaced. 

Two --l-ms reactions before a next step 

Figure 8a shows histograms of dwell times between two main steps 
(90° or 1 20°) that were easily discerned at all [ATP] . At [ATP] < 
(0.02 and 0.2 |ulM), the histograms were fitted with a single expo- 
nential (pink lines) with a rate proportional to [ATP] (Fig. 8b; open 
cirlces), indicating that ATP binding alone sets the pace of rotation 
at these [ATP] . At 2 |ljlM ATP, the histogram starts at the origin at 
zero dwell and shows a distinct peak, indicating the appearance of 
another rate-Hmiting reaction. Between 2 and 60 |xM, the histo- 
grams are explained, roughly, by the ATP-binding reaction and an 
ATP-independent, -0.5-ms"^ reaction (Fig. 8b). At [ATP] » = 
15 ixM, the histograms stiU show a distinct peak, indicating the 
presence of at least two reactions. A fit with two rate constants 
indicated both to be '^l ms (Fig. 8b; filled black cirdes), which, 
taken together, can account for the 0.5 -ms~^ reaction at intermedi- 
ate [ATP]. Thus, at least three rate constants, two for the ~l-ms 
reactions and one for ATP binding, are required to describe the 
rotation kinetics at all [ATP]. Global fit to aU histograms (green 
lines) showed k^n^ to be (3,0 ± 0.1) X 10^ M"^ s"^ (consistent with 
the estimate from Fig. 3 (above) and previous values in actin^^ and 
single-fluorophore^^ assays), and the other two rates to be 1.64 ± 
0.06 ms"^ and 0.71 ± 0.02 ms~\ Of the latter two reactions, we 
cannot discriminate which is first, and the two rates do not differ 
significantly (see also Fig. 8b). 

Thus we propose in Fig. 7a that two '^l-ms dwells separate the 
90° and 30° substeps. A simple explanation is that the first dwell is 
terminated by release of a hydrolysis product (phosphate or ADP) 
and the second by release of the other product. Alternatively, the 
first reaction may be splitting of ATP to ADP and phosphate, and 
the second reaction the release of the two. Either interpretation is 
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Hgure 7 Proposed mechanism for F-i rotation, a, Rotation scheme. ATP with asterisk 
represents ATP or ADP + phosphate; ADP (alone) may be phosphate or ADP + phosphate, 
b, Stepping time courses expected from a. c, Highly schematic diagram for the potential 
energy ^ for 7-subunit rotation. Each coloured line shows ^ in one of the four states in a. 
Tlie orientation of the 7-subunit in state A (in a) is taken as 0°. 
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Figure 8 Dwells between main steps, a, HistDgrams of dwell times between two main 

(90° or 1 20°) steps at various [ATP]. Total counts in eacli liistogram are 60, 463, 1 ,1 45, 
2,862, 631 , 2,384, 3,262 and 1 ,457 In the order of 0.02-6,000 jxM. Histograms for 
Individual runs are distinguished by colours; they are added to constitute a whole 
histogram. Pink lines at 0.02 and 0.2 |xM ATP are single-exponential fits, 
constant-exp(-/cO, with /c shown in open circles in b. Pink lines at other [ATP] are fits with 
two rate constants, constant-[exp(-/QiO-exp(-Ab$], with and shown (filled black 
circles In b). Green lines show the result of a global fit to the Individual histograms {n= 38; 
equal weight for each counQ with sequential reactions fig. 7e) starting with ATP binding 
at the rate /5J''[ATP] and two ATP-lndependent reactions with rates and k2: = (3.0 
±0.1)x10^M""'s"\ k\ = 1.64 ± 0.06 ms~\ and /& = 0.71 ±0.02 ms"^ (s.e.m.). 
b, ATP dependence of the rate constants. Blue circles show the total rate, kor kJ(J 
{ko^+k^ for the individual fits. Green lines show the rate constants obtained in the global 
analysis, and the total rate, [/cj"'' + Z^"* + (/^ri''[ATP]) is shown by the blue line. 



consistent with the biochemical evidence that the releases of ADP 
and phosphate occur at similar rates^\ 

Discussion 

In our rotation scheme (Fig. 7a) one or two of the three catalytic 
sites are fiUed at any time with ATP or its product(s) of hydrolysis. 
This is the so-called bi-site mechanism^'^'*^, which is the norm at 
least at submicromolar [ATP]. Previously we have demonstrated 
rotation at [ATP] as low as 20 nM, indicating that bi-site hydrolysis 
accompanies rotation^ ^ and that bi-site is the fundamental mode of 
rotation^^. Present results suggest bi-site to be the norm also at 
physiological (mM) [ATP], Alternation between two-filled and 
three-fiUed states (the tri-site mechanism) has been proposed for 
hydrolysis at high [ATP], from non-Michaelis-Menten kinetics^^*^ 
as in Fig. 3, and from quenching of tryptophan fluorescence in the 
active sites^. These results indicating the tri-site mechanism may 
have been influenced by the MgADP inhibition^. In contrast, ATP 
synthesis by ATP synthase is insensitive to the inhibition and seems 
to proceed by a bi-site mechanism^^. Our rotation assay focuses only 
on active Fi, and is unaffected by the inhibition. The rotation rate 
was higher than one-third of the hydrolysis rate at aU [ATP] (Fig. 3), 
implying that part of Fj in solution was not fuUy active already at 
several seconds after exposure to ATP. We do not necessarily deny 
possible occupancy of three sites, but we daim that filling aU sites 
does not significantly accelerate, nor add power to, rotation. 

Presence of the substeps indicates the appearance of two meta- 
stable structures during rotation, in which the equilibrium positions 
of the "Y-subunit differ by -^90** (Fig. 7a: A versus B or C). If our bi- 
site interpretation is correct, the crystal structure in Fig. la probably 
corresponds to the two-nucleotide structure B (or C). The other 
central structure with one nucleotide (A) is yet to be solved. A 
structure^^ reported recently indicates that the protruding portion 



of the 'Y-subunit, where the bead was attached, is torsionaUy flexible, 
therefore it is possible that the substeps revealed here might be an 
artefact: intrinsic steps of the ^-subunit are always 120°, whereas the 
bead is somehow obstructed at 90° and lags behind for a few ms. 
We dismiss this possibility because (1) the bead rotated over 120° 
within 0.25 ms at 2m]VIATP; (2) the dwell time at 90° was 
mdependent of the radius of bead rotation, which presumably 
reflects differences in bead attachment; and (3) larger beads did 
not show any sign of obstruction. 

In Fig. 7a, binding of ATP drives the 90° substep. By reciprocity, 
the affinity of the p-subunit for ATP, on the left of the arrow on the 
central 7-subunit, must increase as the 'y-subunit rotates^'^*'*^'^^. 
Likewise, because release of ADP (or phosphate, or both) drives the 
30° substep, the affinity of the p-subunit for ADP, on the right of the 
arrow, must decrease as the 'y-subunit rotates over the last 30°. The 
magnitude of these affinity changes can be estimated as below. We 
have shown that, at least under a high load, Fi does 80-90 pN nm of 
mechanical work per 120° step" and that the torque it produces is 
nearly independent of the rotation angle^*"^^ (the potential energy 
for 7-subunit rotation is linearly downhill). Rotation of 291 -nm 
beads gave similar resvQts (R.Y., unpublished observations). These 
results, combined with Fig. 7a, suggest the diagram in Fig. 7c for H^: 
in state A in Fig. 7a, is minimal at 0°; ATP binding produces state 
B where !Pb is linearly downhiU toward its minimum at ~90°; is 
also minimal at ~90°; product release recovers state A' where Y'a' is 
again linearly downhill toward 120°. Solid lines have the constant 
slope of 80 pN nm of work per 120° as indicated by experiment, 
whereas dashed lines are drawn arbitrarily to provide minima at 
experimental positions. Affinity for ATP is proportional to e3q> [ ( 
- ^b)/^bT1» where k^T '-'4.1 pNnm is the thermal energy at room 
temperature^^. Thus, rotation from 0° to 90° accompanies an 
increase in the affinity of more than exp[(60pNnm)//:B^ '^2x 
10^ The affinity for ADP, ocexp[(¥'A'- ^cV^b^. decreases more 
than exp[(20pNnm)/fcB'r] - 10^ 

These affinity changes account for ATP synthesis by forced 
clockwise rotation of the 'y-subunit. Starting from state A' in 
Fig. 7a, the affinity for ADP of the empty p-subunit on the right 
of the arrow increases as the 'y-subunit rotates clockwise, and this p- 
subunit will pick up ADP from the medium. Further rotation 
decreases the affinity for ATP of the p-subunit carrying the 
previously synthesized ATP (on the left of the arrow), and this 
ATP will eventuaUy be released. In this simple scheme, it is the 
central, asymmetrical -y-subunit that dictates which of the three p- 
subunits should change its affinity, and to what extent. A coroUary 
of the ^-y-dictator modeF is that tri-site operation may take place 
during hydrolysis if ADP release is somehow slowed down, for 
example by partial inhibition. After a 90° substep, the 'y-subunit is 
already pointing close to the next p-subunit and signals this P- 
subunit to bind ATP. ATP binding may therefore occur without 
waiting for the 30° substep accompanying ADP release; however, 
this would be an inefficient tri-site mechanism with no additional 
power and speed, as we claim above. 

Synthesis (or hydrolysis) of ATP on a p-subunit does not require 
much energy, because ATP and ADP + phosphate are in equilibrium 
on the p-subunit^^"^^. This is in accord with our observation that the 
hydrolysis reaction, which we presume to occur during the two ~ 1- 
ms reactions, is mechanically silent. The principal function of 
hydrolysis is to allow the release of bound ATP by converting it to 
products, thereby resetting the machine for the next round of 
stepping. Complete mechanical silence, however, is unfavourable 
for efficient ATP synthesis, because clockwise rotation should 
release ATP, not ADP + phosphate. Indeed, it has been indicated 
that the equilibrium shifts toward ATP during synthesis^. One such 
mechanism (a switch-less model for Fj motor^*'^^) is suggested in 
Fig. 7c, where the minimum in Y'c (after hydrolysis) is placed 
slightly to the right of the minimum in Y'b (before hydrolysis). Such 
a slight shift would be undetectable at the present resolution. 
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Our work is essentially a consolidation and embodiment of the 
binding-change modeF"^. Similar to Fi-ATPase, myosin hosts ATP 
and its hydrolysis product at near equilibrium^*'^". This may be the 
general tactics adopted by many AlP-dependent molecular 
machines. If so, much of the free-energy drop accompanying ATP 
hydrolysis occurs in the ATP-binding step^^'^*'^*', ATP binding maybe 
the principal source of power in these molecular machines. Indeed, 
ATP binding induces large conformational changes in molecules 
such as myosin^^ kinesm'^'^ and chaperonin^. □ 

Methods 

Proteins 

A mutant (a-C193S, 3-Hisio at amino terminus, 7-S107C, -y-OlOC) aaPa^y subcomplex 
(referred to here as Fi) derived from a thermophilic BaciUus PS3 vras biotin^ated at two 
cysteines (7-1070 and 7-2100) by incubation with fourfold molar esrcess of 6-N'-[2-(^- 
maleimido)eth^]-N-piperazinylamidohexyl-D-biotinamide for 1 h at 23 "C. Unbound 
biotin was removed with PDIO (Amersham Pharmacia). Two biotin moieties per protein 
were found in an assay using 4-hydroxya5!;obenzene-2-carboxylic acid^. Eight molar excess 
of streptavidin (Pierce) was added to the biotinylated Fi and purified on Superdex-2Q0HR 
(Amersham Pharmacia). 

Beads 

Colloidal gold (diameter 40 nm; British BioCell International) was coated with biotiny- 
lated BSA by incubating 0.2% colloidal gold in 2 mM potassium phosphate pH 7.0 with 
0.4mgnil~^ BSA and 0.2mgml~^ sulphosuccinimidyl-N-[iS/^'-(D-biotinyl)-6-aniinohexa- 
noyl]-6'-aminohexanoate (biotin- (AC5)2 sulpho-OSu; Dojin) for 1 h at 23 °C. We stored 
modified gold partides in a solution containing 2 mM potassium phosphate and 0.05% 
polyetfa^ene glycoL Polystyrene beads were biotin^dated as follovra: 2.7% amino 
beads (108 nm; Fblyscience) in 20 mM potassiimi phosphate pH 7.0 was incubated with 
1 mgml"^ of biotin- (AC5)2 sulpho-OSu for 1 h at 23 °C. We incubated 2.5% of carboxy 
beads (196 or 291 nm; Bangs) in 20 mM potassium phosphate with 3.6 mM 5-(((N- 
(biotinoyl)amino)hexanoyl)amino)pentylamine trifluoroacetate salt (Molecular Probes), 
1% l-ethyl-3-(3-dimethylaminopropyl)-carbodiiniide hydrochloride and 1% i^-hydro- 
xysulphosuccinimide for 1 h at 23 °C. 

Rate of ATP hydrolysis 

Nucleotide- depleted Fi was prepared^*, and its ATPase activity was determined at 23 °C 
with an ATP- regenerating system*^'^ containing 1 mM phosphoenolpymvate, 200 figml"^ 
pyruvate kinase, 100 |xgml~^ lactate dehydrogenase, O.lSmMNADH, and indicated 
MgATP in buffer A (50 mM KCl, 2 mM MgClj, 10 mM 3- [JV-moipholino] propane- 
sulfbnic add-KOH, pH 7.0). The initial hydrolysis rate was determined ftom the slope of 
absorbance decrease at 340 mi^ the average slope was estimated for the period of 2-5 s 
(2mM-12 jjlMATP), 2-7s (6-2 jjuM) or 2-12 s (0.2-0.06 jJuM) after the start of the 
reaction. In the presence of LDAO, hydrolysis kinetics showed a lag, and thus the rate was 
estimated at steady state during 400-500 s (2mM-60 |jlM), 1,000-1,500 s (20-6 |JuM) or 
4,500-5,000 s (2-0.6 ^M). 

Microsoopy 

Beads of 196 imi and 291 nm were observed with transmitted light on an Olympus DC- 70 
microscope. We observed 40-imi and 108-imi beads with laser dark-field microscopy*^ 
(Fig. Ic). A laser beam (532 nm, diameter 3 mm, 200 mW; Millenia II, Spectra Physics) was 
introduced into a dark-field condenser (mmierical aperture, NA 1.2-1.4; Olympus) to 
illimiinate the specimen obliquely. Light scattered by beads was collected with a xlOO 
objective (NA1.35; Olympus) with its iris diaphragm set to NA —1.1 to block the direct 
ray. The field of view was ~50 (xm. To confirm that we observed single 40-nm beads and 
not their aggregates, vre measured the intensity of beads with regular dark-field micro- 
scopy with a halogen lamp, which provided homogeneous illumination. The intensity 
distribution had a single, large peak and a second small peak at four times the intensity of 
the first peak. Because an object smaller than the wavelength scatters light in proportion to 
the square of its volume'^, these peaks should correspond to single and duplex beads. Most 
beads were thus single. Twenty-fold reduction in the laser intensity did not affect the 
rotation speed, indicating that heating by the laser was insignificant. 

Bead rotaiioii assay 

A flow cell was made of two KOH-cleaned coverslips separated by two spacers with 50 |xm 
thickness". Mixture of 0.1- 1 nM beads and Fi at 10- 100 times the bead concentration in 
buffer B (buffer A plus 5 mgml"^ BSA) was applied to the flow cell. Unbound beads 
were washed out with buffer B, and then buffer B plus MgATP (Mg^"^ 2 mM in excess), 
0.1 mgml~^ creatine kinase, and 1 mM creatine phosphate was infiised. Bead images were 
recorded as an eight-bit AVI file with a fast fiaming CCD camera (HiD-Cam, Nac) at 8,000 
frames per s at [ATP] ^ 2 |jlM; 125 frames per s at 0.2 |jlM; and 60 frames s~^ at 0.02 |jlM. 
The temperature was 23 "C From each unmodified image, the bead centroid was 
calculated as Exi(/i-/th)/2(^i-^th)> where (or/j) is the pixel coordinate, Ii the pixel 
intensity, a threshold value, and the summation was for Ii ^ 

The rotational frictional drag coefficient ^ for the beads was calculated as follows: for a 
single bead of radius a rotating in water with viscosity ij ( = 10~* pN nm~^ s), minimal ^ 



is given by STnja v^en the rotation axis is at the bead centre, and maximal ^ by 
S^rqa^ + 6iriia^ = when the axis is at a bead edg^ for a bead duplex, minimal $ is 

given by 2 X ^w^c^ = ISmia^ for a vertical duplex rotating around its centre, vHiereas 
maximal ^ is 2 x Smic^ + 67nia^ + 6iriia {3af = 76irqc^ for a horizontal duplex rotating 
around an edge. 

Analysis off substep kinetics 

We superimposed and averaged time courses of individual steps as in Fig. 6a, b. First, we 
identified all main (90° or 120°) steps in a continuous run by eye. Then, individual steps 
were aligned on the time axis by positioning the midpoint of each step at time 0 (to within 
± one frame). Vertical alignment was made by shift, by multiples of 120°. Next, steps other 
than the central one were eliminated from individual step records: the part earlier than 

0. 25 ms after the preceding 90° substep (or 120° step vHien substeps were contiguous) was 
replaced with a horissontal line at the angle at the 0.25-ms point (for example, a red line in 
Fig. 6a), and the part later than 0.25 ms before the following 90** substep (or 120** step) was 
replaced with the angle at that point (Fig. 6a; red line). Finally, we averaged all step records 
(Fig. 6a, b; thick cyan lines). 

Step records averaged over several runs (Fig. 6c; coloured lines) were fitted with 
theoretical kinetics (Fig. 6c; grey lines) on the basis of the scheme in Fig. 7a; rate constants 
in the scheme were fixed to the values determined in Fig. 8. For the [ATP] -dependent 
kinetics for the rise from -30° to 0", the grey lines show A3o° exp(A^^[ATP] t) - \i^ere 
A30- is the size of the 30° substep, 1^ = 3.0 x 10^ M~^ s"\ and t (< 0) is time from the 
central 90° substep. For the rise finm 90° to 120°, two reactions with rates ki = 1.64 ms~^ 
and ^2 — 0.71 ms~^ are assumed, giving A3o° B {l — [fc2exp( — Jfc,t) — fciexp( — fc2 01/(^2 — 
fci)} + (120° - A3o°) where B = exp(-^;^J^ [ATP] 0.25 ms) accounts for the loss of some of 
the 30° substeps in the averaging process . Global fit to all e3q>erimental curves (allowing a 
vertical shift for each curve) yielded Ajo" = 29.8 ± 0.3°. 
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Abstract The proton translocating ATP syntliase is conceived 
as a rotatory molecular engine. ATP hydrolysis by its headpiece, 
CFi, drives the rotation of subunit y relative to the hexagonally 
arranged large subunits, (aP)3. We investigated transition states 
of the rotatory drive by polarized confocal fluorometry 
(POCOF) as applied to single molecules of engineered, 
immobilized and load-firee spinach-CFi. We found that the 
hydrolysis of ATP caused the stepped and sequential progression 
of subunit y through three discsrete angular positions, with the 
transition states of y being too shortlived for detection. We also 
observed the stepped motion of e, whereas 5 was immobile as 
(aP)3. 

1998 Federation of European Biochemical Societies. 
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1. Introduction 

The proton translocating ATP synthase is conceived as a 
rotatory molecular engine [1—3]. In the isolated headpiece, Fi, 
and under conditions of ATP hydrolysis the rotation of sub- 
unit y relative to the hexagonally arranged large subunits, 
(aP)3, has been established experimentally by chemical cross- 
Knking [4], by polarized absorption recovery after photo- 
bleaching [5] and by video-microfluorometry [6]. In contrast 
to the first technique, the latter two are time- resolving (see [3] 
for a review) with the following virtues and drawbacks: The 
about 2 |im long fluorescent actin filament in [6] allows the 
direct observation in the microscope of the unidirectional ro- 
tation of subunit y in single molecules of Fi. Because tiie 
filament represents a heavy viscous load on subunit y the 
maximum torque of this motor has been revealed, but finer 
motional details were damped out. The photoselection tech- 
nique in [5], on the other hand, operates with a single dye 
molecule on subunit y, i.e. at a negligible load, but tiie ensem- 
ble averaging hinders a rigorous discrimination between ran- 
dom and unidirectional rotation in some special cases. We 
treated this difficulty by developing a theory of molecular 
stepping motors [7] that has led to the proposal that the 
data in [5,8] might indicate a three-stepped rotation of subunit 
y In the present article we describe a new technique, polarized 
confocal fluorometry (POCOF), that overcomes the remaining 
ambiguity by combining the benefit of single-molecule spec- 
troscopy with the one of a low load on the enzyme. Within 
noise limits it revealed a three-stepped rotation of subunit y 
and a first assessment of the residence times of subunit y in the 



♦Corresponding author. Fax: +49 (541) 969 1221. 
E-mail: junge@uni-osnabrueck.de 



three metastable and tiie more transient interstitial angular 
positions (transition states). 

2. Materials and methods 

A photon-coiJiiting fluorescence microscope, the ConfoCor (Carl 
Zeiss Jena/Evotec, Hamburg, Germany), was used to investigate tran- 
sients of the angular position of the fluorescent probe on a single 
molecule of immobilized Fi. The instrument counts fluorescence 
quanta emitted by single molecules contained in a very small volume 
element ( < 10~^^ 1) as defined by the focus of a laser beam (diameter 

0.5 nm). Fig. IB schematically iUustrates the focus of the ConfoCor 
with a y-labetled molecule of Fi immobilized at the glass/water inter- 
face of a chamber sUde. 

The headpiece of the ATP synthase from spinach, CFi (for chloro- 
plast Fi), was covalently and specifically labelled by a fluorescent 
probe, a maleimide derivative of tetramethyhhodamine, TMR. The 
target cysteine residue was located either on subunit y, e or 5. The 
particular residue on y (C322) was iotrinsic [9,10] whereas the residues 
on e (S85C) and 5 (SIOC) were engineered [11,12]. In the latter two 
cases the engineered subunits were first labelled and then reassociated 
with imlabefled CFi(— 6, e). The labelling did not impair the ATPase 
activity of solubilized CFi. Labelled and solubilized CFi was inuno- 
bilized by adsorption at the water/glass interface of a chamber slide. 
Fig. 1 A shows the respective positions of the label in a current struc- 
tural model of CFi. Fig. 2 shows the selective modifications of the 
subunits as detected by the fluorescence of the bound dye. 

2.1. Preparation of CFj and labelling 

CFi lacking subunits 5 and e (CFi(— 5, e)) was prepared by aaion 
exchange rechromatography [13] of CFi prepared from EDTA ex- 
tracts of spinach diloroplasts [14]. The specific activity of the enzyme 
was 23 U/mg (measured in 50 mM Tris/HCl, pH 7.8, 5 mM ATP, 
2 mM MgCl2, 20% (v/v) methanol, 10 mM Na2S03, 5 min, room 
temperature). CFi(— 5, e) was specifically labelled at y-Cys-322 by 
treatment at room temperature for 10 min at pH 7.0 (25 mM 
MOPS/NaOH) with 20 ^lM tetramethylrhodamine-5-maleimide 
(TMR-5-M, Molecular Probes, Eugene, OR, USA). The reaction 
was terminated by addition of 1 mM iV-acetyl-cysteine. Free dye 
was separated from labelled protein by gelfiltration (Pharmacia 
NAP5/25 mM Tris/HCl, pH 7.8). Recombinant subunits 5 (SIOC) 
and e (C6S, S85C) were over-expressed in E. coli [15] and purified 
by anion exchange and hydrophobic interaction chromatography (5) 
[11] or cation exchange chromatography (e) [12]. Both proteins were 
labefled with tetiamethylrhodamine-5-maleimide similarly as CFi(— 5, 
e) but at 200 TMR and overnight. CFi(— 5, e) was substituted 
witii either 5 or e or both subunits by addition of molar amounts of 5 
and a 6-fold molar excess of E. In the latter case care was taken to 
keep the total amount of urea below 800 mM (e was dissolved in 8 M 
urea). AH samples (CFi(-5, e)+TMR^, CFi(-5, 8)+TMR-5, 
CFi(— 5, e)+TMR-e) were getfiltrated through Pharmacia Superose- 
12 immediately prior to measurements in the ConfoCor. The column 
was cleaned by 6 M guanidinium-HCl weekly. Samples (100 jil), irre- 
spective of the content of small subunits, eluted at Kc = 9.6 ml. The 
specific Mg-ATPase activities remained unchanged upon labelling of 
either subunit y or after substitution with TMR-5 or TMR-e. Ca- 
ATPase activity was decreased from 20 U/mg to 4 U/mg upon 
substitution with TMR-e as expected. 

2.2. The immobilization of labelled CFi 

The immobilization of labelled CFi by adsorption to the glass- 
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buffer interface was monitored by fluorescence autocorrelation spec- 
troscopy, FCS, with the exciting hght beam of the ConfoCor focussed 
into the solution, about 200 ^im above the surface of the chamber 
slide. The concentration was determined by fluorescence autocorrela- 
tion analysis using the algorithms, hard- and software by Evotec 
(Hamburg, Germany). For immobilization of CFi a droplet contain- 
ing 30 |il of 5 nM CFi lacking subunits 5 and 8 and labeUed with 
TMR-5-M within subimit y (CFi(— 5, e)+TMR^) in the presence of a 
100-fold excess of unlabeUed enzyme was deposited on the chamber 
shde. It contained 25 mM Tris/HCl, pH 7.8 and either 1 mM Mg- 
ADP in case of the inhibited enzyme or 20% (v/v) methanol for suc- 
cessive measurements under hydrolyzmg conditions. After about 1 
min the solution was removed from the slide, exchanged by the re- 
spective enzyme-free buffers and left to settle for another 10 min. 
After this procedure, there were no free CFi molecules detectable in 
the focus spot at about 200 jxm depth in the buffer, indicating that the 
immobilization of the enzyme by adsorption to the glass/water inter- 
face was complete. 




confocal volume 




polarized 
exciting light 




Fig. 1. A: A current model of CFi [21]. The arrows indicate the re- 
spective cysteine residues in subunits y, 5, and e that were modified 
with TMR-5-M. B : The focal volume in polarized confocal fluorom- 
etiy. The barrel on the surface of the chamber sUde denotes immo- 
bilized and y-labeUed Fi (not drawn to scale). The arrow denotes 
the exciting laser beam and its plane the E-vector of the hght. 
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Fig. 2. SDS-gel of CFi demonstrating by fluorescence the specific 
labelling of subunits y and 5 and e, respectively. 

2.3. Intersubunit rotation 

Intersubunit rotation was detected with the exciting Hght beam 
directly focussed onto the surface of the chamber slide. TMR was 
excited by continuous linearly polarized focussed argon ion laser, 
514 nm (10 mW). Because of the linear polarization of the exciting 
hght, the count rate varied as a function of the orientation of the dye 
relative to the polarization vector of the beam. This holds true even 
though the bound dye molecule carries out rapid rotational diffusion 
(in nano- to microseconds) around its bond axis relative to the pro- 
tein. This motion diminishes the polarization anisotropy only gradu- 
ally as demonstrated elsewhere [8]. 

3. Results and discussion 

In a first set of experiments the laser focus was positioned 
into a droplet containing solubilized CFi which was deposited 
on a chamber sHde. The concentration of labelled CFi was 
determined by the ConfoCor operating in the FCS mode (flu- 
orescence correlation spectroscopy) [16]. To avoid too tight a 
packing of labeUed molecules at the interface, we used a 100- 
fold dilution of labelled in unlabelled CFi molecules in the 
starting solution. The specific Mg-ATPase activity of immobi- 
lized CFi was decreased 10-fold in comparison with the one of 
the solubilized enzyme. 

After 10 min of incubation the laser focus was readjusted to 
hit immobilized CFi at the glass/water interface. The Confo- 
Cor was then used in a new mode, coined POCOF. The prin- 
ciple of POCOF, which stands for polarized confocal fluo- 
rometry, is illustrated in Fig. IB. When the laser beam 
impinges on the immobilized enzyme molecule in the dou- 
ble-conical focal volume, the angular position of the transition 
moment of the probe (shown as a double arrow) relative to 
the linear polarization (flat arrow) of the laser beam deter- 
mines the probabiHty of excitation (cos2(a) law) and thereby 
the rate of emission of fluorescence quanta. This is why the 
rotation of the label is expected to cause transients of the 
fluorescence count rate. Such transients were recorded over 
an illumination interval of about 5 s. Within this period of 
time the dye molecule in the focus was irreversibly bleached. 
The laser beam was then shut oflf, the focus was horizontally 
shifted (typically by about 10 ^un), and the shutter was re- 
opened for a recording at another spot containing another 
molecule of labelled CFi with a Tresh' dye molecule. The Hght 
stabiHty of rhodamines, the favored dyes in such studies, is of 
the order of 10^-10^ absorbed quanta [17]. 

3.1. Transients with y-labelled CFi 

Transients with y-labelled CFi are depicted in Fig. 3. Fig. 
3A shows the typical result when the hydrolysis activity was 
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Fig. 3. Transients of fluorescence upon excitation with linearly, po- 
larized light of a single, y-labelled and immobilized molecule of 
CFi(— 5, e). A shows a typical transient imder inhibition of the en- 
zyme by ADP. B shows active CFi in the presence of ATP: after 
adsorption of the enzyme to the chamber slide and a 10 min incu- 
bation period in the presence of 50 jxl 25 mM Tris/HCl, pH 7.8 and 
20% (v/v) methanol, 450 ^il reaction mixture for ATP hydrolysis 
(50 mM Tris/HCl, pH 7.8, 5 mM ATP, 2 mM MgClg, 20% (v/v) 
methanol, 10 mM Na2S03) were added to the sample. Note the 
three distinct levels (PI to P3) of the count rate of fluorescence (as 
impUed by the three thin lines) which indicate three discrete angular 
positions of the labelled y relative to (aP)3 during the catalytic 
cycle. The final level, BL, results from the bleaching of the molecule 
in the very focus. 



blocked by the presence of ADP. At zero time, when the 
photoshutter was opened, the count rate started at a level of 
5 kHz. It continued at this level over about 1.5 s, to suddenly 
break-down to a background level of about 2 kHz. From 
there on the count rate continuously decayed over more 
than 10 s. We attributed the count rate of 3 kHz over a back- 
ground of 2 kHz to a particular single molecule of CFi in the 
very focus and the sudden break-down to the irreversible 
bleaching of the covalently attached TMR molecule (see [17] 
for a study on the light stability of TMR). The slower and 
gradual bleaching of the background level was attributed to 
molecules in a more unfavorable position at the fringes of the 
focal volume. The constancy of the count rate until the sud- 
den bleaching occurred, demonstrated that subunit y was not 
rotating under these conditions. This interpretation was 
corroborated by two observations: (1) When the same focal 
spot was reilluminated after a period of darkness, the high 
count rate was never revived. (2) When the microscope table 
was laterally moved to bring a Tresh' molecule into focus, 
high count rates with a sudden break-down were observed 
again. This pattern of a high count rate to start with and a 
sudden drop to a lower level was observed in over 200 shots 
with y-labelled CFi. Under the same conditions we never 
observed a sudden rise of the count rate. The only other 
behavior was a background rate right from the begmning 
because there was no favorably oriented, labelled enzyme in 
the focus. 

Fig. 3B documents the completely different behavior when 
y-labelled CFi was hydrolyzing ATP. The count rate started 
at a low level, PI, to suddenly raise to a higher level, P2, and 
by another jump to P3, then dropping down to PI, and fur- 
ther up to P2, again. Then the level dropped to a basal level, 



BL. The latter was attributed to the background after the 
bleaching of the particular dye molecule in the very focus. 
We interpreted this behavior as indicating the stepwise pro- 
gression of subunit y through three discrete angular positions, 
PI, P2, and P3. A stepped motion of subunit y in the presence 
of ATP was observed in more tiian 200 experiments, a clearly 
three-stepped one only in three. This is not astounding in 
terms of the angular resolution and the noise limits. It is note- 
worthy that e.g. two angular positions that are symmetrically 
placed to tiie interface are indistinguishable by their resulting 
count rates. This will be discussed elsewhere. The noise level 
was strictly proportional to the square root of the counts per 
time channel. Attributable to quantum or shot noise it was 
not to be lessened, because of the limited Hght stability of the 
dye. 

The time base was digitized at intervals of 10 ms. At this 
time resolution tiie transients were instantaneous. The number 
of three steps conformed with the expectation based on (a) the 
number of three reactive sites on the hexagon formed by the 
large subunits of CFi [1,18], and (b) the interpretation of our 
previous experiments with a large ensemble of labelled CFi [5] 
in tiie Ught of a tiieory of molecular stepper motors [7]. It is 
an important new result tiiat tiie angular positions of tiie 
transitions between two sites were so shortiived. Witiiin sto- 
chastic limits, the average persistence on levels PI, P2, and P3 
was much longer than the transit time between them. It was 
compatible with the observed average turnover time of 85 ms 
for ATP hydrolysis by immobilized CFi. 

3,2. Transients of CFi labelled at subunits 5 and e 

Transients of CFi labelled at subunits 5 and e are shown in 
Fig. 4A and B, respectively. In both cases the enzyme was 
hydrolyzing ATP. With the label on subunit 5 we consistently 
(200 cases) observed a behavior with a high count rate at the 
beginning and a sudden drop to a bleached level that was 
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Fig. 4. Transients of fluorescence upon excitation of a siogle immo- 
bilized molecule of CFi with linearly, polarized Ught under ATP hy- 
drolysis (ATP but no ADP added). TMR was covalently attached 
to either subunit 5 (A), or e (B). AU other conditions were as de- 
scribed in Fig. 3. Note that the signal with TMR-labeUed e displays 
the same major characteristic of the signal recorded with active CFi 
which was labelled hi y* The signal with TMR-labeUed 5 resembles 
the signal obtaiaed with y-labeUed but ADP-inhibited CFi (cf. Fig. 
3A). 
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attributable to a non-rotating element (see also Fig. 3A). Con- 
trastingly there were uprising count rates when the label was 
attached to subunit e. Subunit e did not reveal the regularity 
of stepping through three positions as subunit y. Altiiough tiie 
stochastic character of the stepped rotation hampers a direct 
comparison of tiie rotations of e and j, these data demon- 
strate that e belongs to the rotor (with y) and 5 to the stator 
of the enzyme (with (ap)3) (see [21] for an up-to-date struc- 
tural model). The same assignment has been inferred from the 
different effects on the ATPase activity of crosslinldng of 5 
with (ap)3 [11] and of e with y [12,19], respectively. 

Which reaction provides the major driving force for a rota- 
tional step is a key question to understand the driving mech- 
anism. The hydrolysis of ATP is expected to proceed in the 
following sequence: binding of ATP, its sequestration ('tight 
binding'), the hydrolysis of sequestered ATP, and the release 
of ADP and Pj. Since our experiments were carried out under 
ATP saturation, tiie binding of ATP was not rate-limiting and 
certainly not responsible for tiie sudden steps. The hydrolysis 
of ATP by CFi is generally conceived as the inversion of ATP 
synthesis by FqFi, and the release of tightiy bound ATP as 
the major energy requiring step [20]. When Fi operates in the 
reverse direction the sequestration of (loosely) bound ATP 
may be tiie major driving step. It is conceivable tiiat it pulls 
the DELSEED 'lever* at the 'loose' binding site (see [1]) in- 
wardly towards the crank of subunit y to click it round to the 
next angular position. 
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ABSTRACT 

Fluorescence lifetime and anisotropy measurements were 
made on tiie red fluorescent protein (DsRed) from tropical 
coral of the Discosoma genus, both at single-molecule and bulk 
concentrations. As expected from previous work, the fluores- 
cence lifetime of DsRed in solution is dependent on laser 
power, decreasing from an average fluorescence lifetime in the 
beam of about 3.3 ns at low power (3.5 ns if one extrapolates to 
zero power) to about 2.1 ns at 28 kW/cm^. At the single- 
molecule level, exciting with 532 nm, 10 ps laser pulses at 80 
MHz repetition rate, DsRed particles entering the laser beam 
initially have a lifetime of about 3.6 ns and convert to a form 
having a lifetime of about 3.0 ns with a quantum yield of 
photoconversion on the order of 10~^ (calculated in terms of 
photons per DsRed tetramer). The particles then undergo 
additional photoconversion witli a quantum yield of roughly 
10~^y generating a form with an average lifetime of 1.6 ns. 
These results may be explained by rapid photoconversion of 
one DsRed monomer in a tetramer, which acts as an energy 
transfer sink, resulting in a lower quantum yield for photo- 
conversion of subsequent monomers. Multiparameter corre- 
lation and selective averaging can be used to identify DsRed in 
a mixture of fluorophores, in part exploiting the fact that 
fluorescent lifetime of DsRed changes as a function of excita- 
tion intensity. 

INTRODUCTION 

Intrinsically fluorescent proteins such as those derived from green 
fluorescent protein (GFP) have been widely used as intracellular 
probes (1). The widespread use of GFP is attributed to its visible 
excitation and emission bands, mutability to create color-shifted 
variants and the case with which GFP can be fused to other 
proteins and expressed in a variety of genetic backgrounds. More 
recoitLy, another intrinsically fluorescent protein, red fluorescoit 
protein (DsRed) from tropical coral of the Discosoma genus, has 
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been discovered and represents an alternative to GFP (2—4). Like 
some GFP variants, DsRed has a high quantum yield of 
fluorescence (0.7) and a high extinction coefficioit (75000 
cm~^). In addition, the emission from DsRed is substantially 
ledshifted, and it is more resistant to photobleaching than known 
GFP variants (5,6). The low rate of photobleaching and redshifted 
emission should facilitate single-molecule measurements to be 
made in living ceUs because this allows longer signal averaging 
while using excitation and detection conditions that give rise to less 
autofluorescence. 

Widespread use of DsRed as a fluorescent reporter has prompted 
investigation into the biochemical and photophysical properties of 
the protein. It has been shown that DsRed forms a stable tetramer at 
nanomolar concentrations (5-8). The absorbance spectrum peaks 
at 558 nm and the fluorescence spectrum at 583 nm (5). The 
fluorescence lifetime of the native complex has been measured as 
3.6 ns (6). Interestingly, on continued illumination, photoconver- 
sion occurs from this spectral form to redshifted form emitting at 
595 nm and with a variable fluorescence decay time between 1.5 
and 2.7 ns (9). The chemical nature of this photoconverted form is 
not known. Recent single-molecule measurements have been per- 
formed that follow the time course of the photochemical 
conversion process, in terms of both the shift in the emission 
spectmm and of the fluorescence lifetime. In some cases, the 
process appears to occur in steps, possibly because of successive 
conversion of the monomers within the tetramer (9). 

Another interesting aspect of DsRed photochemistry is that the 
anisotropy in bulk solution has been shown to be about 0.23, 
independent of the solvent viscosity (6,10). From the crystal 
structure of the DsRed tetramer, the alignment of the chromo- 
phores indicates that energy transfer may occur on time scales 
roughly an order of magnitude faster than the excited state lifetime 
(10). Rapid energy transfer would depolarize the fluorescence, 
resulting in a low anisotropy for DsRed tetramers, although the 
rotation time in solution of the complex is much longer than the 
fluorescence lifetime (6). Thus, energy transfer between monomers 
in tiie tetramer is apparently an important aspect of DsRed 
photophysics. 

In the measurements that follow, the photophysical and 
photochemical properties of DsRed are further investigated. Most 
of the measurements described are single-molecule fluorescence 
studies. Unlike previous single-molecule work on DsRed, these 
measurements arc performed on molecules free in solution, thus 
avoiding potential surface effects. To obtain a reasonable signal-to- 
noise ratio during the brief transit of the molecules through the 
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Figure 1. Schematic of the microscope showing the microscope objective 
(a), dichroic mirror (b), emission filter (c), pinhole (d) and broadband 
polarizing beam splitter (e). The synchronization detector and electronics 
(Sync.) are used to generate a synchronization signal for the time-correlated 
photon-counting board (Timeharp). The router multiplexes the input signals 
from different detectors (Dl and D2) so tiiat tiiey are stored in different time 
ranges (time relative to the synchronization pulse in nanoseconds). 

laser beam, selective averaging methods are used, which allow one 
to correlate processes such as photoconvoision with time sp^t in 
the beam and fluorescence lifetime while maintaining knowledge 
of the distributions of the parameters involved. Finally, the ability 

to distinguish between DsRed and other fluorophores in a mixture 
using the property of photoconversion as a distinguishing 
characteristic for DsRed is explored. 

MATERIALS AND METHODS 

Instrumentation. Single-molecule spectroscopy was performed using 
a home-built, confocal microscope system and a high-fiequency, mode- 
locked laser as the excitation source to &cilitate time-correlated measurements 
of fluorescence lifetime. A schematic of this system is shown in Fig. 1. The 
excitation laser was a fi^uency-doubled, pulse-conqxressed, mode-locked 
Nd: YAG laser (532 nm, 10 ps pulse diuation) with a repetition rate of 82 MHz. 
All fitting of fluorescence decays was performed by convolution of a single- 
exponential decay with tiie instrument response function of the system. The 82 
MHz repetition rate gives a small degree of "wrap-around" fluorescence, but 
the contribution of photons that are detected more than 1 2 ns after the excitation 
pulse does not change the fluorescence lifetime of the samples studied in this 
work. The error in lifetime calculations is primarily due to the poor signal-to- 
noise ratio inherent in single-molecule fluorescence measurements. To ensure 
proper beam quality and polarization, the laser light was passed through 
a single-mode, polarization-preserving glass fiber (F-SPA; Newport, Irvine, 
CA) and a polarizing beam splitter (0SBC1SFH.3; Newport). The laser ligjit 
was delivered to an inverted, confocal microscope and reflected up toward the 
microscope objective witii a dichroic mirror (Q570LP; Omega Optical, 
Brattieboro, VT). The sample, a 50 microliter droplet, was spread onto a glass 
coverslip (22 X 50 mm, no. 1.5; VWR, West Chester, PA). The same objective 
(lOOX PlanApo 1.4NA; Olympus, Tokyo, Japan) used to focus the laser also 
collected the fluorescence. The collected fluorescence passed through the 
dichroic mirror and was focused onto a 100 |4.m diameter pinhole (9 lO-PHlOO; 
Newport). The fluorescence was then split by a polarizing beam splitter 
(05FC16PB.3; Newport), sending photons polarizedparallel to thelasertoone 
detector (SPCM-AQR-12; Perkinr-Ehner, Fremont, CA) and photons 
polarized perpendicular to the laser to a second detector. To reduce Raleigh 
and Raman scattering, die fluorescence was passed througfh a custom emission 
filter (Omega Optical). The filter specifically blocks 532nmlightand the water 
Raman scattering from 532 nm light. 

The signal from the detectors and a synchronization signal from each 
laser pulse were sent into a multiplexer (router), which was designed and 
built locally. The multiplexer sends a start signal and a stop signal into 
a Timeharp time-correlated single-photon counting (TCSPC) card 
(Picoquant, Berlin, Germany). The multiplexer separates the signals from 
the two detectors in regions of time defined by the repetition rate of the laser 
such tiiat signal of each detector occupies a different 12 ns region of time 
recorded by the Timeharp. 



Sample preparation. Single colonies of DH5a Escherichia coli cells 
expressing DsRed from the plasmid pDSRED (Qontech, Palo Alto, CA) 
were used to inoculate twenty 3 mL tuttes of L-Broth media (by weight: 1% 
bactotryptone and 1% NaCl and 0.5% Bacto yeast extract) and were grown 
for 48 h at 3TC with shaking. Eighteen of the 20 tubes turned red due to 
DsRed expression and subsequent assembly and chromophore formation. 
The media from these tubes, 54 mL in total, were spun in an SA 600 rotor 
(Sorvall) at 8000 rpm for 8 min at 20°C. The cells were washed with buffer 
(10 mM Tris-HCl, pH 7.4, 100 mM NaCl, 1 mM MgCl2, 10 mM 
dithiothreitol), respun and suspended in 10 mL of the same buffer. The cells 
were then broken by passing them through a French pressure cell at 
pressures greater tiian 20000 psi. The efQuent from the pressure cell was 
spun at 8(XX) rpm for 8 min in SA 6(X) rotor. The supernatant was collected 
and respun at 1 5 000 rpm for 1 5 min in S A 6(X) rotor, and the supernatant of 
this spin was collected. To remove permanent aggregates, the product was 
filtered twice using a 100 kDa, nominal molecular weight cutoff membrane 
filter (Centricon) spun at 2500 rpm in SA 600 rotor for 30 min. Aliquots of 
this sample were stored at — 80°C until use. The sample was diluted for 
single-molecule measurements in 10 mM Tris, 1 mM ethylenediaminetetra- 
acetic acid buffer to approximately 0.1 nM. Higher concentrations were 
used for measuring bulk properties of the protein. 

Note that in these samples, most of the DsRed do not pass through the 
100 kDa cutoff membrane, as would be expected for tetramer DsRed 
complexes. This filtration was performed to remove the very large 
complexes from the mixture (anytiiing significantiy larger tiian a tetiamer 
is very unlikely to pass through), although botii the single-molecule 
(analyzed in terms of fluorescence correlation spectroscopy and diffusion 
times) and dynamic light scattering measurements indicate that larger 
complexes reform dynamically after filtration (data not shown). 

RESULTS 

The power dependence of the ensemble-averaged 
fluorescence lifetime 

For measuring ensemble-averaged fluoiescence lifetime as a func- 
tion of laser intensity, 200 nAf DsRed was used. The isotropic 
lifetime of the excited state, as measured by TCSPC and fit with 
single-exponential least squares curve fitting, is strongly dependent 
on the intensity of the 532 nm excitation light, changing from 3.3 
ns at the lowest power measured to about 2.1 ns at 28 kW/cm^ 
(Fig. 2). The focal path has a diameter of 0.8 \xm\ consequently, 
very high light intensities are achieved with relatively low laser 
powers. Extrapolating to zero power gives a lifetime of 3.5 ns, in 
good agreement with the previously measured value of 3.6 ns at 
low power (6). The steady-state anisotropy for the ensemble- 
averaged measurements remained 0.22, independent of laser 
intensity (data not shown). 

Distinguishing between single-particle fluorescence 
and background using multiparameter measurements 

The remaining measurements were performed at the single- 
molecule (single-particle) leveL A technical issue in single-molecule 
spectroscopy is distinguishing between photons generated by 
fluorescence from individual particles in the beam and those 
from other sources. This is perhaps most critical in experiments 
such as these, where diffusion of particles in and out of the beam 
plays a major role in the dynamics. Typically, the average amount 
of fluorescence that can be detected from a single dye molecule 
when it diffuses through the probe volume under saturating ex- 
citation, using the instrument described here, is between 30 and 
50 photons, although this depends on the fluorescence yield as 
well as on the dark-state lifetime and 3deld. In addition to the 
fluoiescence from single particles, there are background signals 
from Raman scattering by the solvent and electronic noise in the 
detection electronics. 
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Figure 2. The fluorescence lifetime measured in bulk samples depends on 
the intensily of Hght used for excitation. Fluorescence lifetimes were 
measured as described in tfie text for 70 nM DsRed on the same microscope 
described for the single-molecule work. 



The multidimensional nature of the data collection procedure 
used in this work makes it possible to resolve the single-molecule 
or single-particle fluorescence fiom the background signal most of 
the time by sorting the photons or photon buists based on their 
measured characteristics. In these measurements, three types of 
information were recorded for each photon: the time of arrival in 
the lab frame, the time of arrival relative to the excitation laser 
pulse (i.e. the lifetime of the excited state that gave rise to the 
photon) and the polarization of the photon. 

The most obvious discriminating characteristic between back- 
ground fluorescence and fluorescence from a single fluorophore in 
the probe volume is the number of photons detected within a certain 
amount of time (in the lab time frame). This can be characterized 
by the average closed time between photons (the inverse photon 
density [IPD]). Photons making up the background signal will have 
more time between them than photons making up the single- 
molecule events. Consequentiy, the IPD for photons in the 
background will be larger than for photons that derive from 
singile-molecule events. 

The IPD for each detected photon was calculated by determining 
the average difference betweoi the arrival time of a given photon 
and the arrival times of the nine photons closest to it in time. For 
the level of background and burst sizes in these data, 10 photons 
were sufficient to give a sharp change in density when a molecule 
was in the beam and were large enough so that background 
photons did not give a density high enough to be considered 
a single-molecule burst 
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where Tq is the arrival time of the photon under consideration, and 
Tj is the arrival time of the nine closest photons in time. 

The distribution of IPD is given in Fig. 3e. There are two clear 
populations of photons defined by tiiis distribution. Most of the 
photons observed have a high inverse density (a long time between 
photons) on the order of 0.5 ms. A smaller population of photons has 
a low inverse density (a short time between photons), mostiy less 
than 100 ^s. By determining the average values of otfao: measurable 
parameters for photons as a function of IPD, it is possible to tell 
which IPD region contains the single-molecule fluoresc^ce signals 
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Figure 3. Selectively averaged fluorescence parameters are shown as 
a function of the IPD (all photons within a narrow IPD range are used to 
determine each parameter shown in panels a-d). The IPD is the average 
amount of time el^sed between detected photons. An IPD of greater than 
0.1 ms was used as the threshold for single-molecule events in later 
analyses. As a function of the IPD, the amount of uncorrelated signal (the 
constant term in a fit of the fluorescence decay) (a), exponential 
fluorescence decay lifetime (b), pre-exponential amplitude of fluorescence 
decay component from the fit (c), anisotropy (d) and relative number of 
photons detected (e) are shown. 

of interest and in tills way rejectmost of tiie background signal. Todo 
this, the value of a particular parameter is detramined for all photons 
having a particular range of IPD values. The fluorescence anisotropy, 
R, is calculated according to the equation 



R = (GI, - 1, )/(GI, + 21, ) 



(2) 



where G is a correction factor for balancing the parallel and 
perpendicular channel detection efficiencies, I|| is tfie number of 
photons having a polarity parallel to the laser and Ij, is the number 
of photons having a polarity perpendicular to the laser. The result is 
shown in Fig. 3d. 
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It is also possible to measure the arrival time of each photon 
relative to the mode-locked laser pulse and for a particular IPD 
range and to determine the average lifetime of the excited state that 
gives rise to those photons. Figure 3a,c shows the nondecaying 
amplitude, lifetime and pre-exponential amplitude of a least 
squares curve fit of the isotropic fluorescence decay (that is, the 
sum of the parallel and twice the perpendicular components of the 
fluorescence decay) to the equation 

F(t) = Ae"*/*^ + B (3) 

where F(t) is the number of detected photons at time t measured 
relative to the laser pulse, A is the pre-exponential amplitude, x is 
the decay time constant and B is the constant (nondecaying) term 
accounting, in these data, for photons not correlated with the laser 
pulse. 

Photons due to Raman scattering should have a high fluores- 
cence anisotropy and arrive within the laser pulse. Background 
signals caused by electronic noise are not temporally correlated 
with the laser pulse and should have a zero anisotropy. Looking 
first at the fluorescence lifetime as a function of IPD (Fig. 3b), it is 
clear that the dominant fluorescence lifetime becomes very short 
for the photons in the main (near 0.6 ms) region of the inverse 
photon distribution. This is consistent with a fast component in this 
part of the IPD curve that arises finom Raman scattering. The same 
region of the curve also shows a large amplitude for the constant 
term in the fit of the isotropic fluorescence decay (Fig. 3a). This is 
consistent with a fraction of the photons with this inveree density 
arising from electronic noise that is not correlated with the mode- 
locked laser pulse. Thus, in the broad region of high IPD centered 
around 0.6 ms, there appears to be a mixture of photons due to 
Raman scattering and electronic noise. For this reason, the average 
anisotropy value observed for this part of the distribution is 
intermediate between the high value expected for scattering and the 
zero value expected for signal due to uncorrelated electronic noise. 

The photons with inverse densities of less than 100 \is have 
much longer fluorescence lifetimes on average (Fig. 3b), as would 
be expected for the DsRed particles (9). In addition, the value of 
the constant term in the fluorescence decay fit is almost zero (veiy 
little uncorrelated electronic noise at this IPD). The anisotropy for 
IPD values less than 100 ^s varies around 0.2. Although this is not 
very different fi^m the anisotropy of the noise, it is also about what 
is expected for the anisotropy of DsRed (6,10). Although 
fluorescence correlation spectroscopy results obtained from this 
sample can be explained in terms of a single size of diffusing 
particles (data not shown), the variation in anisotropy between 
0 and 300 us IPD may be due to heterogeneity in DsRed particle 
size and composition. Generally, the greater the number of 
energetically equivalent DsRed molecules in the aggregate the 
lower the fluorescence anisotropy is likely to be because the energy 
will equilibrate betweoi many different monomers with different 
transition dipole orientations. Larger aggregates, or aggregates that 
lack a monomer that has undergone photoconversion to a short 
fluorescence lifetime form, would be likely to give more photons 
per unit time and thus lower IPD. Thus, one might predict that 
within the low-IPD population, the lowest IPD should be 
associated with the lowest anisotropies. 

By analyzing the IPD as a function of the other measured 
parameters as described above, it is possible to discriminate 
between the single-particle fluorescence signals and much of the 
background noise using an IPD threshold. An IPD threshold of 



0.1 ms was used in all the studies reported here. The low— invCTse- 
density photons thus selected are then reassembled into discrete 
single-molecule events. This is performed by chronologically 
scanning tiie remaining photons in tiie data set and considering any 
photon with no photon before it for at least 2 ms to be the start of 
an event and any photon with no photon for at least 2 ms after it to 
be the end of an event All singile-particle bursts identified in this 
way that contained at least 40 photons were included in the 
analyses described below. The use of these parameters for single- 
molecule event identification is not new (11); however, they have 
not been applied to the analysis of DsRed. 

Monitoring the photochemical conversion 

of DsRed at the single-molecule level in solution 

Previous work has shown that as DsRed is exposed to Hght, it 
undergoes a photochemical conversion (9,10,12). To foUow this 
process in single molecules in solution, photons in the bursts were 
biimed according to the time that elapsed from the initiation of the 
single-particle fluorescence burst. Then, the photons from all 
events with a particular bin time were averaged together, giving 
a representation of the fluorescence as a function of time after 
a single particle entered the beam. This is sirmlar to performing an 
average measurement in a bulk experiment as a function of time 
after the onset of illumination, except that tiie experiment can be 
performed continuously without having to synchronize the opening 
of the shutt^ with the data collection, and the distributions of the 
measured parameters and their average are available. Figure 4 
shows the result of such an experiment and analysis monitoring the 
average fluorescence lifetime as a function of time in the beam. 
These data were collected at a series of diff<mnt laser intensities 
from 4 to 40 kW/cm^ (using a 80 MHz mode-locked, doubled 
Nd:YAG laser at 532 nm for excitation). At the lowest power 
(bottom panel. Fig. 4), the lifetime is initially near 3.6 ns, as has 
been measured at low laser power in bulk solution. This rapidly 
decreases to about 3.0 ns over a few hundred microseconds. At 
highor powers (middle and upper panels. Fig. 4), the lifetime is 
shorter to begin witii and then drops with time. At the highest 
power, tiie fluorescence lifetime is initially about 2.2 ns and 
decreases with time to about 1.6 ns. The drop in fluorescent 
lifetime of DsRed as a function of time in the beam is consistent 
with previous single-molecule measurements performed on 
immobilized samples and thus likely represents the photoconver- 
sion process (9). 

Detecting DsRed particles in a mixture using 
selective averaging and multiparameter correlation 

One long-range goal of this work is to be able to better foUow 
DsRed in complex mixtures, such as cells. As shown above, 
selective averaging of collected photons can be used to distinguish 
between signal and background or to correlate one parameter, such 
as fluorescence lifetime, with another, such as time spent in the 
beam, in order to determine the time depcnidence of photochem- 
istry. To explore the uses of selective averaging and multiparam- 
eter correlation in order to distinguish DsRed from other 
fluorescent molecules, mixtures of DsRed and rhodamine-6G (R- 
6G) dye molecules, both of which at low power have similar 
fluorescence lifetimes, were used as a model system. 

To selectively average fluorescent burst parameters and to 
perform multiparameter correlations, it is necessary to determine 
the values of these parameters for individual fluorescent bursts. 
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Figure 4. Photons from individual single-particle bursts are collected into 
groups as a function of time relative to the start of die event For each group 
of 1000 photons (the first 1000, second 1000, etc.), the fluoiescenoe decay 
lifetime is calculated (y-axis) and plotted as a function of tiie average arrival 
time of the photons in that group relative to the point tiie DsRed particle 
entered the beam (x-axis). At different excitation intensities, the rate of 
photoconversion varies. While in the beam, the particles convert into a form 
having a shorter lifetime. For each plot, tiie excitation intensity used is 
shown in tiie top right comer of tiie panel. 

The estimation of single-moleciile fluoiesc^ce lifetimes of in- 
dividual bursts was performed using the maximum-likelihood 
estimator (MLE) optimization technique. MLE is an exceptionally 
good, curve-fitting technique for analyzing small numbers of 
photons (13). The lifetimes are calculated from the isotropic decay 
(calculated as the sum of the time-correlated photons in the vertically 
polarized channel plus twice the photons in the horizontally 
polarized channel) to avoid polarization distortion of tiie decay 
times. The steady-state anisotropy is calculated for each burst 
according to Eq. 3. The transit time through the beam for a particular 
single-particle event is the difference in arrival times between the last 
and the first photon in a burst. The average IPD is the arithmetic 
average of the IPD for the individual photons in the burst. 

Figure 5 shows the results of a selective averaging experiment in 
which the average values of the anisotropy, burst duration and IPD 
(milliseconds between photons) are shown for aU bui^sts with 
a particular fiuoiescence lifetime (panels a-c). Figure 5d shows 
a simple histogram of the numbers of bursts with different 
fluorescence lifetimes. For the data shown here, a high excitation 
power was used that had the effect of shifting the average 
fluorescence lifetime of the DsRed bursts to shorter times, whereas 
having no effect on the R-6G lifetime. Two peaks in the lifetime 
histogram are observed, near 2 and 4 ns. Looking at the selectively 
averaged parameters (Fig. 5a— c), the one most obviously correlated 
with the two lifetime populations is the anisotropy. R-6G has an 
anisotropy near zero, whereas, as discussed above, the anisotropy 
of the DsRed on average is near 0.23 (6,10). 

Direct correlation between the lifetime and anisotropy in- 
formation for individual bursts allows yet better separation of the 
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Figure 5. The average anisotropy (a), burst duration (b) and IPD (c) are 
shown as a function of fluorescence decay lifetime for a mixture of DsRed 
particles and R-6G. The value of each average parameter was determined 
by summing the photons from all events having the fluorescence decay 
lifetime shown on tiie x-axis. Also shown is the number of events recorded 
as a function of fluorescence lifetime (d). 

two populations of fluorophores in the mixture. In Fig. 6, the 
anisotropy and the fluorescence lifetime of each burst are used to 
generate a two-dimensional contour plot, with the height of the 
contour representing the number of bursts with a particular lifetime 
and anisotropy. The first two panels show the result of separate 
two-dimensional correlations between anisotropy and lifetime for 
DsRed and R-6G. The final panel shows the mixture. When both 
parameters are simultaneously taken into account, the bursts can be 
cleanly separated into the two populations with almost no overlap. 

DISCUSSION 

Applying selective averaging to single-molecule spectroscopy 

In general, there are at least two ways to use information fnom 
single-molecule fluorescence studies. Perhaps, the most obvious is 
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Figure 6, MultqKinuneter contour plots 
of the number of events with a particular 
anisotropy and fluorescence lifetinie for 
R-6G, DsRed and a mixture of the two. 
The two-component histograms allow 
better resolution of the components in 
the mixture. 
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to follow the dynamics of individual molecules for periods of time, 
studying dynamics directly. Examples include the movemoit of 
molecular motors or following the individual steps of a multistep 
reaction such as RNA pol3mierization or protein folding (14-18). 
Although very powerful, such methods are limited at the single- 
molecule level in terms of signal to noise, time resolution (it is 
difficult to collect more than about 100 photons/ms from a single 
molecule) and usually involve subjective decisions about the 
molecules that should be used to make mechanistic conclusions. 

There is another broad set of analytical methods that can be 
applied to single-molecule data. This includes various forms of 
single-molecule correlation techniques and other selective averag- 
ing techniques. These types of analyses take advantage of the fact 
that when fluorescence data are collected at the single-molecule 
level, it is left to the investigator to decide what parameters to 
average or correlate and what parameters to consider independently 
for each molecule in the population. For example, one can 
determine the average anisotropy of all the bursts that have 
a particular fluoresc^ce lifetime or determine the fluorescence 
amplitude as a function of time for only those bursts that have 
a certain anisotropy. Single-molecule fluorescence correlation 
spectroscopy is another common example of this type of approach. 
Such analyses reveal correlations between observable parameters 
that are lost in an ensemble average. Although any averaging 
results in some loss of information, selective averaging can play 
a very important role in error reduction, particularly in cases where 
the amoimt of signal from each individual molecule is low or 
where one wants to look for mechanistically important trends that 
apply to a subpopulation of molecules or states of molecules in 
a mixture. 

Selective averaging and multiparameter correlation techniques 
have been used in several different ways in this work. First, 
selective averaging of certain fluorescent properties (hfetime, 
anisotropy, etc.) was performed as a function of IPD. In tiiis way, it 
was possible to objectively determine which photon densities 
corresponded to single-particle fluorescence signals versus back- 
ground scattering and electronic noise (see Results). Photons were 
then selected for further analysis based on this crit^on. The 
second example involved generating the time dependence of the 
photoconversion process using events that occun-ed at random 
intervals due to diffusion through the laser beam. By determining 
when each burst initiates and thus when each molecule enters the 
excitation volume, it is possible to line up aU the events in time and 
look for the time-dependent processes that occur on average. The 
final example involved correlating multiple parameters from 
individual bursts with one another to determine the average 
properties of the individual subpopulations within a mixture, in this 



case a mixture of DsRed and the dye, R-6G. The results from the 
second and third examples will be discussed in more detail below. 

Photoconversion of DsRed particles in solution 

The fluorescence lifetime of DsRed in solution depends on the 
laser power used for excitation (Fig. 2). One might think that this 
could be because of energy transfer and singlet— singlet axmihilation 
among closely spaced DsRed monomers (rapid energy transfer 
presumably accoimts for the low anisotropy of the DsRed tetramer 
[6]), but previous work has shown that preillumination of DsRed 
samples also results in shorts fluorescence lifetimes (9). This 
clearly indicates that photoconversion of DsRed between two or 
more forms is involved, as concluded by others. Past measure- 
ments have shown that extended irradiation of DsRed results in an 
apparently irreversible conversion of the 583 nm fluorescing form 
witii a 3.6 ns fluorescence decay time to a form tiiat fluoresces with 
a peak at 594 nm and has a rather heterogeneous fluorescence 
decay time between about 1.5 and 2.7 ns (9). The heterogeneity of 
the final state has been shown by single-molecule spectroscopy. 
Immobilized DsRed particles that were followed individually 
during the photoconversion process arrived at different final 
fluorescent lifetimes (final in the sense that they preceded 
photobleaching of the particle) and sometimes showed intermedi- 
ate fifetime st^s during the photoconversion process (9). 

Evidence for a multistep photoconversion process is also given 
in Fig. 4. In this case, the DsRed sample was free in solution. This 
has the disadvantage that it is more difficult to follow the dynamics 
of individual molecules because of rapid diffusion but the 
advantage that it avoids effect due to the interaction of the DsRed 
with a surface, which could itseif generate lifetime heterogeneity. 
In addition, solution measurements provide a simple means of 
studying microsecond to millisecond photodynamics. For the data 
shown in Fig. 4, fluorescence was recorded as a function of time 
after particles entered the beam and monitored for as long as they 
remained in the beam (typically himdreds of microseconds to 
rmlliseconds). Under high levels of illumination (about 40 kW/ 
cm^. Fig. 4, top panel), one can see that within the time resolution 
of the measurement, photoconversion of the 3.6 ns, native DsRed 
molecule to a form with a hfetime of about 2.2 ns had already 
taken place (the earliest lifetimes measured were roughly 2.2 ns). 
Subsequentiy, a much slower photoconversion takes place over at 
least 5 ms, resulting in a 1.6 ns form. Thus, there are clearly at least 
two processes at work here. A high^-yield photoconversion step 
that results in an intermediate-lifetime species and then a lower- 
yield photoconversion resulting in an additional lifetime decrease. 
The initial photoconvmion process can be seen more clearly at 
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low power. Using 4 kW/cm^ (Fig. 4, bottom panel), one sees that 
the lifetime of DsRed starts at essentially the low-power fluo- 
rescence lifetime value (nearly 3.6 ns) and then decreases on a 
time scale of roughly a few hundred microseconds to about 3.0 ns. 

It is possible to estimate the quantum yields of the fast and 
slow photoconversion steps from the data presented here. From 
the extinction coefficient of DsRed (75000 cm~^ at 560 nm 
or about 38000 Af~* cm~* at 532 nm), one can calculate an 
absoibance cross section at 532 nm of about 1.4 X 10~** cm^ for 
a monomer or 5.7 X 10~ cm^ for a tetramer. If we take the initial 
photoconversion time at 4 kW/cm^ intensity to be about 300 p.s 
(Fig. 4, bottom panel), then this implies that about 6.8 X 10~*^ J of 
light is absorbed by a tetramer during the initial photoconversion 
period, which corresponds to about 1.8 X 10^ photons at 532 nm 
absorbed during this time. Thus, the quantum yield of the initial 
photoconversion in a tetramer is on the order of 10~^. The 
subsequent photoconversion event(s) occurs on the millisecond or 
longer time scales at 3- to 10-fold higher Hght intensity, suggesting 
that the average quantum yield for subsequent photoconversion 
steps is one to two oniers of magnitude less tiian tiiat for the initial 
pholoconvereion process. (The estimation of the range of photo- 
conversion efficiencies for later steps ignores any effect of the 
changing absorbance cross section of the tetramer during the 
photoconversion process. This could become significant especially 
after two or three of the molecules in the tetramer have undergone 
photoconversion.) 

One possible explanation of the multistep photoconversion 
process is as follows. Initially, tiie DsRed tetramers consist only of 
583 nm fluorescing monomers and possibly immature green forms, 
depending on the growth and isolation procedures used (5,19). 
Energy transfer between monomers assures that the energy is 
r^idly localized on a 583 nm fluorescing monomer, which has an 
excited state lifetime of 3.6 ns. This long-lived excited state is 
reactive and can undergo photoconversion with a quantum yield of 
nearly 10 . After photoconversion of the first monomo: in the 
tetramer to the 595 nm fluorescing form, a tetramer is created in 
which subsequent excitation results in rapid energy transfer from 
the 583 nm fluorescent species to the photoconverted 595 nm 
fluorescent species. The 595 nm species has a shorter lifetime 
(perhaps about 1.5 ns), and thus the total fluorescence decay time 
consists of the energy transfer time plus the decay time of the 
photoconverted species, giving rise to a fluorescoice decay for the 
tetramer that is int^mediate between the 3.6 ns decay time of 
native DsRed and that of the pure photoconverted form. Because of 
this energy transfer to the photoconvCTted form, the excited states 
of the 583 nm fluorescent monomers are now shorter lived and thus 
less likely to undergo photoconversion themselves (there are also 
fewer of them to absorb light). Therefore, the rate of the 
photoconversion process for the remaining 583 nm fluorescing 
monomois is slower than for the initial monomer. As photo- 
conversion slowly occurs and more and more of the 583 nm 
monomors are converted to monomers fluorescing at 595 nm, the 
overall fluorescence lifetime continues to decrease because more 
and more of the photons are directiy absorbed by the photo- 
converted form, and thus the lifetime approaches that of the pure 
photoconverted DsRed tetramer. Of course, as the number of 583 
nm fluorescing monomers decreases, the total amount of time that 
these monomers are likely to spend in the excited state becomes 
less, both because of a decreased absorbance cross section and 
because of the &ct that there are a greats numbo: of long- 
wavelength monomers to rapidly transfer energy to. The result is 



that the overall yield for continued photoconversion from the 583 
nm (long fluorescence lifetime) to 595 nm (short fluorescence 
lifetime) form in the tetramer continues to decrease. 

Using the photochemical and photophysical 
properties of DsRed to distinguish it 
from other fluorophores in solution 

One of the obvious uses of DsRed is as an intracellular probe in the 
same sense that GFP is used widely for this purpose. It has been 
suggested that DsRed may have some advantages in this regard, 
other than simply providing a long-wavelength version of an 
intrinsically fluorescent protein. DsRed has a lower rate of 
photobleaching than GFP as well as a high extinction coefficient 
and large quantum yield of fluorescoice (5,10). Some of its other 
photophysical properties may also prove useful in this regard. The 
fact that it undergoes photoconversion could be used as a tool in its 
identification. During a confocal scan of a cell, for example, setting 
the excitation power effectively sets the average fluorescence 
lifetime of DsRed. Thus, one could either set the lifetime to a value 
different from contaminating fluorophores or perform differential 
scans witii low and tiien high power, looking for tiie part of tiie 
signal that changes lifetime between the scans. Of course, it is often 
the case for ceU imaging that fluorescent proteins are expressed as 
fusion proteins. Under these conditions, aggregation to form the 
tetramer may not occur or the aggregation state may be different 
from what is observed for DsRed alone. However, there is nothing 
to suggest that the photoconversion event depends on the 
aggregation state (other than in terms of cross section as described 
above), and thus although the details of the power dependence of 
tiie fluorescence Ufetime may change, tiie abihty to tune tiie 
fluorescence lifetime with laser power should be qualitatively 
similar. 

An example of setting the DsRed lifetime and performing 
measurements at the single-molecule level in mixtures is shown in 
Fig. 6. Here, two fluorophores (DsRed and the laser dye R-6G) 
with similar fluorescence lifetimes at low power were mixed in 
solution. By using a high-intensity laser beam to probe the system, 
it was possible to tune the fluorescence lifetime of DsRed away 
from that of R-6G, using the photoconversion process described 
above. The fluorescence lifetime and the anisotropy (that should 
also be diffaient for these two molecules because of the size 
difference) were measured for each particle as it passed through the 
beam. By correlating these two parameters, it was straightforward 
to assign individual bursts as either DsRed or R-6G, with nearly 
complete confidence. In this case, the measurement of either 
anisotropy or fluorescence lifetime alone would allow event 
identification most of the time, but by combining the two 
properties, the two subpopulations become almost completely 
resolved. As shown in Fig. 5, both burst duration and IPD are also 
correlated with the change in lifetime, and thus correlations with 
these two parameters could have been used to further distinguish 
between the two subpopulations. The analytical methodology used 
for such multiparameter correlations has been worked out 
previously (20). As described above, this is a case where single- 
molecule detection techniques become extremely useful even when 
average information is the desired goal. By performing the 
measurements at the single-molecule level, it is possible to 
correlate parameters with one another and in this way sort 
fluorophores into subpopulations, facilitating the determination of 
average properties of multiple subpopulations within a larger 
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population. This is not easily done in most cases from ensemble- 
averaged measurements. 
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The atomic force microscope (AFM) is a powerful tool for imaging 
individual biological molecules attached to a substrate and placed 
in aqueous solution. At present, however, it is limited by the speed 
at which it can successively record highly resolved images. We 
sought to increase markedly the scan speed of the AFM, so that in 
the future it can be used to study the dynamic behavior of 
biomolecules. For this purpose, we have developed a high-speed 
scanner, free of resonant vibrations up to 60 kHz, small cantilevers 
with high resonance frequencies (450-650 kHz) and small spring 
constants (150-280 pN/nm), an objective- lens type of deflection 
detection device, and several electronic devices of wide band- 
width. Integration of these various devices has produced an AFM 
that can capture a 100 x 100 pixeP image within 80 ms and 
therefore can generate a movie consisting of many successive 
images (80-ms intervals) of a sample in aqueous solution. This is 
demonstrated by imaging myosin V molecules moving on mica (see 
http://www.s.kanazawa-u.ac.jp/phys/biophys/bmv_movie.htm). 

One of the advantages of the atomic force microscope (AFM) 
(1) is its capacity to image individual biomolecules in, say, 
a buffered solution containing ions at physiological concentra- 
tions (2, 3). Such capacity suggests that the instrument can be 
used to record the dynamic behavior of such molecules. In 
practice, however, only very slow processes can be recorded (2, 
4-6), because commercially available AFMs require minutes to 
form an acceptable image, and many interesting biological 
processes occur at much higher rates. To understand, and 
overcome, the factors that limit the scanning rate of an AFM, we 
begin by considering relations between the characteristics of the 
constituting components. 

We consider only the "tapping mode" of AFM operation 
(Digital Instruments, Santa Barbara, CA). This is the mode 
suitable for imaging biological macromolecules, because vertical 
oscillation of the cantilever at (or near to) its resonance fre- 
quency reduces lateral forces between the tip and the sample (7). 
The oscillating tip briefly taps the surface at the bottom of each 
swing, resulting in a decrease in oscillation amplitude. During the 
x-y scan of the sample stage a feedback loop (see below) keeps 
this decrease (and hence the tapping force) constant; this is 
necessary for minimizing the deformation of soft samples. The 
error signal — the difference between a preset signal and the rms 
amplitude of the cantilever — is fed into a proportional-integral- 
differential (PID) feedback circuit. The PID output is amplified 
and then sent to the z-piezo actuator; this is repeated until the 
error signal returns to zero. For the three-dimentional move- 
ment of the sample stage to follow the sample topography 
accurately, the bandwidth of the feedback loop should be 
comparable to, or larger than, the frequency determined by the 
x-y scan velocity and the apparent width of the features on the 
surface. To increase the imaging bandwidth, all elements in the 
feedback loop have to be optimized. The first element in this 
loop is the cantilever. A high resonance frequency in liquid is 
required to generate a high tapping frequency; also required is 
a low spring constant, so as to minimize the tip-sample interac- 
tion force. These conflicting requirements can be met only with 
small cantilevers (8, 9). Every circuit in the feedback loop must 



have a wide bandwidth. The last element, the z-scanner, is a 
piezo actuator with a high resonance frequency, which in as- 
sembly must react very quickly without unwanted mechanical 
vibrations. Viani et al (10, 11), using the tapping mode, have 
used small cantilevers with resonance frequencies ^15Q kHz to 
image DNA and GroEL-GroES complexes on mica, in liquid. 
Their frame rate was less than 1 s~^ because they used a 
conventional piezo-tube with a resonance frequency of «*15 kHz 
in the z direction, together with feedback electronics of a limited 
bandwidth. Suldiek et aL (12), using a cantilever with an 
integrated zinc oxide piezo actuator (13), achieved an imaging 
bandwidth of 38 kHz and frame rate of «*4/sec. However, tibeir 
cantilever was too stiff to image soft samples. In the AFM 
assembly described below we have tried to optimize each ele- 
ment for its use in imaging samples in solution, using fast tapping 
operation, and an objective lens type of deflection detection 
device suitable for small cantilevers. Assembly of these selected 
components results in an imaging bandwidth of ^60 kHz (fastest 
among the high-speed AFMs that use the tapping mode and are 
suitable for examining soft samples in liquid). Combination of 
this high speed with zooming-in capability (reduction of the scan 
area to 240 X 240 nm^) allows viewing of liquid samples of 
proteins on mica at a rate of 12.5 s~^. 

Test Sample. Components, and Imaging 

Sample Preparation. For test purposes, myosin V was extracted 
from chick brains and purified as described (14). Myosin V 
stored at 0°C in buffer A (25 mM KCl/25 mM Imidazole, pH 
7.6/2 mM MgCl2/l mM EGTA/0.2 mM DTT) was diluted to 7 
nM with buffer A. A drop (^2 /xl) of myosin V solution was 
placed on freshly cleaved mica (1 mm in diameter) for 3 min, 
rinsed with buffer A to remove unattached myosin V, and then 
imaged in 2 mM ATP plus buffer A without DTT. 

Cantilevers. The cantilevers, made from silicon nitride using 
micromachining techniques, were 140 nm thick, 2 fjm wide, and 
9-11 fjum long. The rear side of the cantilevers was coated with 
gold 20 nm thick. All surfaces of the cantilevers were further 
coated with osmium, about 2 nm thick. The tips were grown by 
electron-beam deposition (15), with growth rate of about 5 nm/s. 
The tip length was adjusted to '=«1.0 fjum. The radius of the tip end 
was 6-8 nm. The mechanical properties of the cantilevers were 
tested by measuring the spectra of their thermal motion. The 
resonance frequencies were 1.3-1.8 MHz in air, and 450-650 
kHz in water, and the spring constants were estimated to be 
150-280 pN/nm. 

Deflection Detection System. We used an optical lever method to 
measure the deflection of a cantilever. The optics must generate 
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Fig. 1. Schematic drawing of the AFM head integrated with an inverted type 
of optical microscope, (a) A/4 Plate, {b) Dichroic mirror, (c) Polarization beam 
splitter. (cO Focusing lens, (e) Split photodiode. {i) Collimator lens, (g) Laser 
diode, {h) Illumination lamp. (/) Collimator lens, (j) Half mirror, {k) Viewing 
system of the optical microscope. (/) Screw for adjusting the cantilever height, 
(m) Piezo for exciting cantilever, (n) Objective lens, (o) Cantilever, (p) Sample 
stage, (g) Sample stage of the optical microscope. All of the components (a-g 
and n) of the deflection detection system are hung down from the sample 
stage of the optical microscope. The scanner (see Fig. 2 also) is mounted on the 
same sample stage. The collimated laser beam is reflected up, incident on the 
objective lens at an off-centered position. The outgoing beam thus tilted 
about 10° from the vertical plane is focused onto the cantilever whose plane 
is tilted about 10° from the horizontal plane. The beam reflected at the 
cantilever is collimated by the objective lens, separated from the incident 
beam by the polarization beam splitter with A/4 wave plate, and reflected 
onto the split photodiode. The optical microscope allows us to view the 
cantilever and the focused laser spot. The specimen is su pported at the bottom 
of the sample stage (p), and an inverted cantilever (o) probes the specimen 
from below. 



a very small incident beam spot, because our cantilevers for rapid 
scanning are much smaller than commercial^ available cantile- 
vers. Schaffer et al. (16) designed an AFM using optics for small 
cantilevers, wherein a laser beam reflected from the rear side of 
the cantilever is collected and collimated with the same lenses as 
are used for focusing the incident laser beam onto the cantilever. 
We used the same method, but instead of using single lenses we 
used an objective lens (CFI Plan Fluor ELWD 20xQ Nikon) with 
a long working distance of 8 mm. As shown in Fig. 1, our AFM 
is integrated with an inverted type of optical microscope (1X50, 
Olympus, Tokyo). The focusing objective lens also is used to view 
the cantilever and the focused laser spot with the optical 
microscope. The focused spot is 2-3 /Ltm in diameter. In the 
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Fig. 2. Scanner assembly. {Lower) The side view when looked at the scanner 
from +x to —X. The scanner has a two-layered structure that guarantees that 
the x-scan and the y-scan do not interfere with each other. For structural 
details, see the text. 



design of Schaffer etaL, the entire optical train is tilted from the 
vertical so that the incident beam is normal to the plane of the 
cantilever. Instead of this, we used a much simpler design 
wherein the vertical laser beam is incident on an off-centered 
position of the objective lens to make the outgoing beam axis 
normal to the plain of the cantilever. The incident and reflected 
light beams were separated by polarization using a beam splitter 
and a quaterwave plate. The laser beam, reflected, collimated, 
and then slightly converged, is fed into the two closely spaced 
photodetectors (S2721-02, Hamamatsu Photonics, Hamamatsu 
City, Japan) whose photocurrents are sent to a differential 
amplifier with a wide bandwidth (2.5 MHz). 

Scanner. The scaimer is shown in Fig. 2. Stack piezoelectric 
actuators (AE0203D04, Tokin, Tokyo) are used in this scanner. 
They have the resonance frequency of 260 kHz in free oscillation, 
their maximum displacement is 4.5 /Lim, and their capacitaace is 
90 nF. The y-piezo displaces base-1 on which the x- and 
z-scanners are mounted, while the x-piezo displaces base-2 on 
which the z-scanner is mounted. The z-scanner has two z-piezo 
actuators placed in the opposite direction to one another (the 
reason is described below). Base-1 is glued to the y-piezo, and 
base-2 is glued to the x- and z-piezo actuators. A sample stage 
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Fig. 3. Circuit for fast amplitude measurement. The output sinusoidal signal 
from the split-photodiode amplifier is fed to this circuit. The output of this 
circuit provides the amplitude of the sinusoidal input signal at periodicity of 
the input signal. For details, see the text. 



is attached to one of the z-piezo actuators via a tliin layer of 
vacuum grease. When the z-piezo is displaced quickfy, hydro- 
dynamic force is generated as the reaction from the sample 
solution to the sample stage. To minimize this reactive force, a 
glass of a circular-trapezoid shape with a small top surface of 
1-mm diameter is used as the sample stage. To hold the mass 
balance, a dummy stage that has the same mass as the sample 
stage was attached to the other z-piezo actuator. Base-1 sits on 
three steel ball bearings, and its top plane is pushed down by 
plate-1 through three steel ball bearings. The balls are 1 mm in 
diameter. Base-2 is also held between base-1 and plate-2 in a 
similar manner. This way of holding base-1 and base-2 between 
two flat surfaces reduces vibrations of these bases in the 
z-direction and allows smooth displacements in the desired 
directions. However, these arrangements are not sufficient to 
minimize such vibrations. Quick displacements of a z-piezo exert 
impulsive forces on base-2, which causes vibrations of base-2, and 
in turn vibrates the z-piezo itself. Such impulsive forces are 
countered by the simultaneous displacements of the two z-piezo 
actuators of the same length, in the counter direction. 

Electronics. The custom-made piezo-drive amplifiers have a 
1-MHz bandwidth for a 90-nF capacitive load, in 2 Vp-p output. 
These amplifiers are used to drive the z-piezo actuators and to 
drive another piezo actuator (custom-made; with a 2-MHz 
resonance frequency) that excites the cantilever in the tapping 
mode AFM operation. The x- and y-piezo actuators are driven 
by amplifiers (HJPZ-0.15Px3, Matsusada Precision, Tokyo) with 
a 25-kHz bandwidth for a 90-nF capacitive load, in 100 Vp-p 
output. The PID and two types of rms-to-dc circuits were 
custom-made. The bandwidth of the PED circuit is 300 kHz. An 
analog type of rms-to-dc converter with a fourth-order low-pass 
filter can accept up to a 1-MHz input signal and requires five 
waves of the signal to output an accurate rms value; in other 
words, five tip-sample contacts are needed to measure the 
tapping amplitude. To increase the conversion rate, we made a 
novel circuit with 1-MHz bandwidth, which requires only one 
wave to generate the amplitude signal (Fig. 3). The input 
(source) signal is divided into two. One is sent to two sample/ 
hold (S/H) circuits. The other is sent to a 90° phase shifter, and 
its output is fed to a zero-cross comparator to generate rectan- 
gular wave signals. These rectangular signals are sent to a 
mono-stable multivibrator to generate pulses at the top and 
bottom peaks of the input source signal. The pulses from the 
mono-stable multi vibrator are used as triggers to hold the source 
signal by the S/Hs. The outputs from the two S/Hs are sent to 



a differential amplifier to generate output that represents the 
amplitude of the source signal. The signal generation for the x-y 
scan, and the data sampling are made in different ways. One is 
made by using a digital-signal-processor (DSP) control board 
integrated with A/D and D/A converters. This allows flexible 
controls for AFM imaging, but its operation rate is limited 
because of our slow DSP system. The other is made by using an 
external circuit that generates analog signals for the x-y scan and 
timing pulses for data acquisition, and by using a fast A/D board 
for data sampling. This system allows a fast operation (the pixel 
clock can go up to 1 MHz). 

AFM Head Assembly. The sample stage of the inverted type of 
optical microscope was replaced with a stainless steel plate 10 
mm thick. Every mechanical and optical component of the AFM 
is mounted on or hung from this plate, lliis design reduces 
relative movement between these components, despite vibra- 
tions of the optical microscope itself, wMch is susceptible to the 
ambient sounds or mechanical noises because of its low reso- 
nance frequency. Because we use conf ocal optics for focusing the 
laser beam onto a cantilever, as well as for viewing the cantilever 
and the incident spot with the optical microscope, the rear side 
of the cantilever must face downward (and hence, the tip must 
stick upward). This setup is opposite to that widely used in 
commercial AFMs. The cantilever is attached to a transparent 
glass plate via a short cantilever holder glued to the glass plate. 
The scanner, with a sample stage facing downward, is mounted 
over the cantilever (see Fig. 1). 

Imaging Bandwidth. When a piezo actuator is oscillated with its 
one end being fixed, its resonance frequency is half the resonance 
frequency when it is oscillated with bolji ends free. So, the 
z-piezo actuators must have a resonance frequency of 130 kHz, 
without considering the mass attached to the free end of the 
actuators. The z-scanner in the whole scanner assembly must 
have a lower mechanical resonance frequency, considering the 
mass (of the order of g), the size (of the order of mm), and the 
Young's modulus (of the order of W^-W^ N/m^) of the 
materials such as super invar steel and resin used in the piezo 
actuators. Therefore, the z-scanner is the component that is most 
difficult to optimize for high-speed scanning. We tested the 
z-scanner's ability to follow a frequency modulation of the 
driving signal, without a feedback system, while the tip of a 
cantilever remained in contact with a mica surface. The sensor 
output is shown as a function of the driving frequency (Fig. 4). 
When the z-piezo for the counterbalance was not operating, the 
z-scanner showed resonance vibrations around 8.5 kHz, 34 kHz, 
and 100 kHz (Fig. 4a). The resonance at 100 kHz is caused by the 
resonance of the z-piezo actuator itself. The extra amplitudes 
gained by resonance at the lower frequencies are, however, 
relatively small compared with the oscillation amplitude without 
resonance, indicating the value of our scanner design in sup- 
pressing the resonance vibrations. When the counterbalancer 
was operating (Fig. Ab), the extra amplitudes at the lower 
resonance frequencies were reduced, and we further noticed 
much less variation in the amplitude over frequencies lower than 
60 kHz, compared with those without counterbalance. 

The bandwidth of the whole feedback loop system was as- 
sessed by measuring the system's ability to follow a modulation 
of the set point while the tip was in contact with a mica surface 
(12). Sinusoidal signals were sent to the set point input, and the 
output of the PED was fed to a lock-in amplifier to measure the 
amplitude and the phase relative to the input sinusoidal signals. 
The output amplitude was almost constant up to 60 kHz, and 
then gradually decreased with increasing frequency. This grad- 
ual decrease is caused by the resonance oscillation of the 
z-scanner around 100 kHz. The phase difference was zero over 
the frequencies tested (up to 100 kHz). These results indicate 
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Fig. 4. Oscillation of the z-scanner as a function of the driving frequency. The 
amplitude of the driving signal was kept constant while sweeping its fre- 
quency. The data were obtained with (b) or without (a) operating the z-piezo 
for the counterbalance. 



that the component that has the lowest bandwidth in the Sfystem 
is the z-scanner, and therefore, the bandwidth of the system is 
about 60 kHz. 

The minimum imaging time (t) per frame with pixels of size 
p is given by 



[1] 



where fs is the bandwidth of the feedback loop, and w is the 
apparent width of features on the surface, w is given by convo- 
lution of the tip with radius R and the sample with radius r, and 
is approximated asw = \'rR (17). A scan range of 200 nm and 
100^ pixels, thus a pixel size of 2 nm, can be scanned in 33 ms for 
iv = 10 nm and fs = 60 kHz. Because the corresponding pixel time 
is 1.65 /xs, this scan condition requires the tapping and data 
sampling frequencies to be equal to, or larger than 606 kHz (the 
reciprocal of the pixel time). One step of the x-y scan has to be 
completed within this pixel time. 

Imaging. We examined whether imaging can really be made at (or 
near to) the maximum rate expected above. Myosin V directly 
attached to mica, in solution containing 2 mM ATP, was imaged 
successively (240-nm scan range; 100^ pixels) for 4 sec (50 
frames). The scan rate was 1.25 kHz, corresponding to a tip speed 
of 0.6 mm/s, and the frame rate was 12.5 s~^. In the myosin V 
molecule located on the center of the images (Fig. 5), its two 
head/neck regions, one long tail, and the globular tail end are 
clearly resolved. The tail seems to be about 100 nm long, longer 
than a previous estimate based on its primary structure and the 
electron microscopic images (18). The angle of the long tail 
relative to the head/neck regions changed between the eighth 
and ninth frames (see the third and fourth images in Fig. 5). 
After this quick change, the tail and the tail end are slowly 



Fig. 5. Successive images of myosin V on mica in buffer solution. The same 
area of 240 x 240 nm^ was imaged 50 times with 100 x 100 pixels. Only nine 
successive images are shown. The tip speed is 0.6 mm/s (scan rate, 1.25 kHz), 
and the frame rate is 1 2.5/s. The tapping frequency is 620 kHz, the amplitude 
of the cantilever's free oscillation is 12 nm, and the set point is about 11.5 nm. 
In the first image, the head/neck regions, the long tail, and the globular tail 
end are marked with arrow 1, arrow 2, and arrow 3, respectively. A mirror 
image of one of the head region is marked with arrow-4 in the second image. Ili^l 
The still panels here appear noisier than the dynamic images in the movie on ll^ii 
the web site http://www.s. kanazawa-u.ac.jp/phys/biophys/ bmv_movie.htm. $i$xS 
The two head/neck regions are overlapped in the z direction, and therefore, SSiSs 
they are not well resolved and do not show a typical "Y" shape. Dynamic iiiiii^ii: 
images of myosin V that show "Y" shape are also presented at http:// 
www.s.kanazawa-u.ac.jp/phys/biophys/bmv_movie.htm. 



moving. The reconstructed movie and the other movies that 
dramatically show the changes in the locations of the myosin 
heads, neck, and tail regions can be viewed at http:// 
www. s .kanazawa-u .ac.jp /phys/biophy s /bmv_movie .h tm . We 
tried faster imaging (25 frames/sec for 240-nm area with 100 X 
100 pixels), but it was not successful because of the inertia of the 
mass attached to the free end of the x-piezo. This inertia 
produces mirror images on the left half area of the total image. 
Note that even in Fig. 5 mirror images are noticed (indicated by 
arrow-4 in the second image). 

Discussion 

The components that limit the imaging rate of our AFM are first 
the x-scanner, and second the z-scanner and the cantilever. An 
upper limit of the imaging rate that is reachable by improving 
these components can be estimated. We used super invar steel 
for base-2. To minimize the inertia of the mass, we may be able 
to use hard ceramics for base-2, which can reduce the mass to a 
quarter, and therefore can increase the imaging rate four times. 
Cantilevers smaller than those we made can be constructed. For 
instance, a silicon nitride cantilever, 42 nm thick, 6 fjum long, and 
2 fjum wide, may be possible, although this will be probably the 
smallest size, considering how it has to be used in practice. This 
small cantilever with a gold coat will have the likely resonance 
frequency of ^2.5 MHz in air and the spring constant of ^200 
pN/nm. The resonance frequency in water must be around 1.2 
MHz, about two times larger than that of our present cantilevers. 
In the tapping mode operation with this cantilever the pixel 
imaging time can be reduced to 0.83 /lis, which corresponds to 17 
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ms of the imaging time, for 100^ pixels. If the second harmonic 
oscillation can be used, the imaging time can be reduced to 8.3 
ms. The speed limit of our z-scanner arises from the resonance 
frequency of the piezo actuator. Piezo actuators with higher 
resonance frequency are available, but their maximum displace- 
ment is shorter. Considering the necessary displacement size, 
say, 1 /Ltm, the maximum resonance frequency will be about 1 
MHz when one end of such a piezo actuator is fixed. This is about 
seven times larger than that used in this study. Considerations 
made here suggest that an upper limit of the imaging time seems 
to be 10—20 ms, as long as we use the optical deflection detection 
method as well as the tapping mode operation. 

The long tail of a myosin V molecule changed its orientation 
in less than '=^30 ms (just before the change, the tail end located 
aty ^ 214 nm (i.e., at almost the end of the frame scan), while 
just after the change the tail end located at y 64 nm). On the 
other hand, after the quick change the tail moved slowly. This 
marked contrast suggests that this quick orientational change 
may be driven by ATP. However, at present we must reserve 
judgment, because we could not observe repeated orientational 
change with the same myosin V molecule. To decide this, we 



have to record images for a much longer time, because the 
ATPase rate of myosin V alone is very small ('=«0.05/s). At 
present our AFM system cannot store many images at once. In 
the near future we may be able to resolve this problem. 

In life science, it has been a dream to view the nanometer-scale 
dynamic behavior of individual biopolymers in solution. The 
capacity of acquiring successive images every 80 ms will allow a 
large expansion in the scope of biological processes that can be 
examined in real time. In the near future we should be able to 
see the behavior of processive motors [such as kinesin (19) and 
myosin V (20)] walking along their tracks, or of molecular 
chaperones assisting a polypeptide chain to fold (21). Such direct 
observations will be a great help in understanding the mecha- 
nisms by which biomolecular machines operate. 

We thank Keicfai Mukai and Ryuki Tsunekawa for their help in making the 
AFM head, Takeshi Sakamoto and Tatsi^ Kinosita for preparing myosin 
V, and Professor Manuel Morales (University of the Pacific) for his 
continuous encouragement and for reviewing this manuscript. This work 
was supported by the Proposal-Based New Industry Creative Type Tech- 
nology R&D Promotion Program from the New Energy and Industrial 
Technology Development Organization (NEDO) of Japan to T.A 
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Abstract Surface plasmon resonance (SPR) biosensors 
are optical sensors exploiting special electromagnetic waves 
— surface plasmon-polaritons — to probe interactions be- 
tween an analyte in solution and a biomolecular recogni- 
tion element immobilized on the SPR sensor surface. Ma- 
jor application areas include detection of biological ana- 
Ij^es and analysis of biomolecular interactions where SPR 
biosensors provide benefits of label-free real-time analyt- 
ical technology. This paper reviews fundamentals of SPR 
afOnity biosensors and discusses recent advances in de- 
velopment and applications of SPR biosensors. 

Keywords Optical sensor ■ Biosensor • Affinity 
biosensor • Immunosensor • Surface plasmon resonance 



Introduction 

The last two decades have witnessed remarkable progress 
in the development of affinity biosensors and their appli- 
cations in areas such as environmental protection, bio- 
technology, medical diagnostics, drug screening, food safety, 
and security. An affinity biosensor consists of a transducer 
(electrochemical [1], piezoelectric [2], or optical [3]) and a 
biological recognition element which is able to interact 
with a selected anal3^e. Various optical methods have 
been exploited in biosensors including fluorescence spec- 
troscopy [4], interferometry (reflectometric white light in- 
terferometry [5] and modal interferometry in optical 
waveguide structures [6]), spectroscopy of guided modes 
of optical waveguides (grating coupler [7] and resonant 
mirror [8]), and surface plasmon resonance (SPR) [9, 10]. 
Fluorescence-based biosensors offer high sensitivity but, 
due to the use of labels, they require either multi-step de- 
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tection protocols or delicately balanced affinities of inter- 
acting biomolecules for displacement assays, causing sen- 
sor cross-sensitivity to non-target analytes [11]. Sensors 
such as optical interferometers, grating coupler, resonant 
mirror, and SPR rely on the measurement of binding-in- 
duced refractive index changes and thus are label-free 
technologies. This paper focuses on SPR biosensor tech- 
nology, reviews fundamentals of SPR sensing, and dis- 
cusses advances in development and applications of SPR 
biosensors and emergmg trends in SPR biosensing. 



Fundamentals of surface plasmon resonance (SPR) 
biosensors 

Surface plasmon-polariton 

A surface plasma wave (SPW) or a surface plasmon-po- 
lariton is an electromagnetic wave which propagates along 
the boundary between a dielectric and a metal, which be- 
haves like quasi-free electron plasma [12, 13]. An SPW is 
a transverse-magnetic (TM) wave (magnetic vector is par- 
allel to the plane of interface) and is characterized by the 
propagation constant and electromagnetic field distribu- 
tion. The propagation constant of an SPW, jS, can be ex- 
pressed as: 



(1) 



D 



where CO is the angular frequency, c is the speed of light in 
vacuum, and and are dielectric functions of the di- 
electric and metal, respectively [12, 13]. This equation de- 
scribes an SPW if the real part of £m is negative and its ab- 
solute value is smaller than Sq, At optical wavelengths 
this condition is fulfilled for several metals of which gold 
is most commonly employed in SPR biosensors. The real 
and imaginary parts of the propagation constant describe 
spatial periodicity and attenuation of an SPW in the direc- 
tion of propagation, respectively [12]. The electromag- 
netic field of an SPW is confined at the metal-dielectric 
boundaiy and decreases exponentially into both media. 
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Fig. 1 Distribution of the magnetic field intensity for an SPW at 
the interface between gold and dielectric (refractive index of the 
dielectric 1.32) in the direction perpendicular to the interface cal- 
culated for two different wavelengths 



Fig. 1 . For an SPW at the boundaiy between gold and a di- 
electric with a refractive index of 1 .32 the penetration depth 
(the distance from the interface at which the amplitude of 
the field falls to \le of its value at the metal surface) into 
the dielectric is typically 100— 500 nm in the visible and 
near infrared regions [10, 13]. 



Concept of surface plasmon resonance biosensing 

Owing to the fact that the vast majority of the field of an 
SPW is concentrated in the dielectric, the propagation 
constant of the SPW is extremely sensitive to changes in 
the refractive index of the dielectric. This property of 
SPW is the imderlying physical principle of affinity SPR 
biosensors — biomolecular recognition elements on the sur- 
face of metal recognize and capture analyte present in a 
liquid sample producing a local increase in the refractive 
index at the metal surface. The refractive index increase 
gives rise to an increase in the propagation constant of 
SPW propagating along the metal surface (Fig. 2) which 
can be accurately measured by optical means. 



^^^^^ -^ii 








Fig. 3 Surface plasmon-poleiriton probing: (a) biomolecul£ir inter- 
action occurring within a short distance from metal surface, and 
(b) biomolecular interaction occurring within the whole extent of 
the field of a SPW 



The magnitude of the change in the propagation con- 
stant of an SPW depends on the refractive index change 
and its distribution with respect to the profile of the SPW 
field. There are two limiting cases: 

1. analyte capture occurs only within a short distance 
from the metal surface (Fig. 3a), and 

2. analyte capture occiurs within the whole extent of the 
SPW field (Fig. 3b). 

Perturbation theory [14] suggests that if the binding oc- 
curs within the whole depth of the SPW field (Fig. 3b), the 
binding-induced refractive index change. An, produces a 
change in the real part of the propagation constant, AjS, which 
is directly proportional to the refractive index change: 



Re{A)5}= kdm 



(2) 



where k denotes the free-space wavenumber [15]. If the 
refractive index change is caused by a binding event oc- 
curring within a distance from the surface d, much smaller 
than the penetration depth of the SPW, the corresponding 
change in the real part of the propagation constant can be 
expressed as follows: 



Re{A>5}= 42lM^A« = Fk^n 



(3) 



Fig. 2 Principle of SPR biosensing 



where % and denote the refractive index of the biolayer 
and the refractive index of the background dielectric me- 
dium (sample), respectively. The binding-induced change 
in the propagation constant of the SPW is proportional to 
the refractive index change and the depth of the area 
within which the change occurs. The factor F (F<V) ac- 
counts for the fact that the interaction occurring within a 
thin layer is probed by only a fraction of the field of the 
SPW. 



Excitation and interrogation 
of surface plasmon-polaritons 

In SPR sensors, an SPW is excited by a light wave and the 
effect of this interaction on the characteristics of the light 
wave is measured. From these meastirements, changes in 
the propagation constant of the SPW can be determined. 
Excitation of an SPW by light can occur only if the com- 
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Fig. 4 Excitation of surface 
plasmon-polaritons: (a) by a 
light beam via prism coupling, 
(b) by a guided mode of opti- 
cal waveguide, and (c) by light 
diffraction on a diffraction 
grating 



a) ^ 
wave 
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ponent of the light's wave vector that is parallel to the 
metal surface matches that of the SPW. This can be achieved 
by means of prism coupling, waveguide coupling, and grat- 
ing coupling. 

In configurations based on prism coupling a light wave 
passes through a high refractive index prism and is totally 
reflected at the prism— metal layer interface generating an 
evanescent wave penetrating the metal layer (Fig. 4a). This 
evanescent wave propagates along the interface with a prop- 
agation constant which can be adjusted to match that of 
the SPW by controlling the angle of incidence. This method 
is referred to as the attenuated total reflection (ATR) method 
[13], The process of exciting an SPW in an optical wave- 
guide-based SPR structure (Fig. 4b) is similar to that in 
the ATR coupler. The light wave is guided by an optical 
waveguide and, when entering the region with a thin 
metal layer, it evanescently penetrates through the metal 
layer exciting an SPW at its outer boundary. Altematively, 
an SPW can be excited by diffraction on a grating. Fig. 4c. 
The component of the wave vector of the diffracted waves 
parallel to the interface is diffraction-increased by an amount 
which is inversely proportional to the period of the grating 
and can be matched to that of an SPW [16]. The interac- 
tion of a light wave with an SPW can alter light's charac- 
teristics such as amplitude, phase, polarization and spec- 
tral distribution. Changes in these characteristics can be 
correlated with changes in the propagation constant of the 
SPW. Therefore, bmding-induced changes in the refrac- 
tive index at the sensor surface and, consequently, the 
propagation constant of the SPW can be determined by 
measuring changes in one of these characteristics. Based 
on which characteristic is measured, SPR biosensors can 
be classified as angle, wavelength, intensity, phase, or po- 
larization modulation-based sensors. In SPR sensors with 
angular modulation the component of the light wave's 
wavevector parallel to the metal surface matching that of 
the SPW is determined by measuring the coupling strength 
at a fixed wavelength and multiple angles of incidence of 
the Ught wave and determining the angle of incidence 
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Fig. 5 Reflectivity and phase for light wave exciting an SPW in 
the Kretschmann geometry (SF14 glass prism — 50 nm thick gold 
layer - dielectric) versus (a) the angle of incidence for two differ- 
ent refractive indices of the dielectric (wavelength 682 nm), and 
(b) wavelength for two different refractive indices of the dielectric 
(angle of incidence 54°) 
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yielding the strongest coupling (Fig. 5a, upper plot). In 
SPR sensors with wavelength modulation the component 
of the light wave's wavevector parallel to the metal surface 
matching that of the SPW is determined by measuring the 
coupling strength at a fixed angle of incidence and multiple 
wavelengths and determining the wavelength yielding the 
strongest coupling (Fig. 5b). In SPR sensors with intensity 
modulation, the change in the intensity of the light wave in- 
teracting with the SPW is measured at a fixed wavelength 
and angle of incidence (Fig. 5b). In SPR sensors with phase 
modulation, shift in phase of the Ught wave interacting 
with the SPW is measured at a fixed wavelength and angle 
of incidence (Fig. 5a, lower plot). In SPR sensors with po- 
larization modulation, changes in the polarization are mea- 
sured at a fixed wavelength and angle of incidence. 



Performance characteristics 

The main performance characteristics relevant for SPR 
biosensors include sensitivity, accuracy, precision, repeata- 
bility, and the lowest detection limit. Sensor sensitivity 5, 
is the ratio of the change in sensor output, P (e.g. angle of 
incidence, wavelength, intensity, phase, and polarization 
of light wave interacting with an SPW) to the change in 
measurand (e.g. analyte concentration, c). SPR biosensor 
sensitivity can be decomposed into two components — 
sensitivity to refractive index changes produced by the 
binding of analyte to biomolecular recognition elements 
on the sensor surface 5ri, and the efficiency with which 
the presence of analyte at a concentration c is converted 
into the change in the refractive index ni 



<3 — — ojy^ 

dn dc 



(4) 



The efficiency E depends on the properties of the biomo- 
lecular recognition elements and the target analyte. The 
refractive index sensitivity can be decomposed into 
two terms: 



3Re{>5} dn 



— SiS2 



(5) 



The first term depends on the modulation method and 
the method of excitation of an SPW [17, 18, 19, 20]. The 
5^2 term is independent of the modulation method and the 
method of excitation of the SPW and describes the sensi- 
tivity of SPW's propagation constant to the refractive in- 
dex change, Eqs. (2) and (3). 

Accuracy describes the degree to which a sensor output 
represents the true value of the measurand (analyte con- 
centration). Accuracy is often confused with precision which 
refers to the way in which repeated measurements con- 
form to themselves without a reference to any true value. 
Repeatability refers to the capacity of a sensor to repro- 
duce output reading under the same measurement condi- 
tions over a short interval of time. The lowest detection 
limit describes the lowest concentration of analyte that 
can be measured by the sensor. 



SPR biosensing formats 

An interaction between a biomolecular recognition ele- 
ment on an SPR sensor surface and analyte in a liquid 
sample is governed by kinetic equations. In order to illus- 
trate fundamental properties of the interaction, we shall 
discuss the pseudo-first-order kinetic equation: 



6t 



= Kc{l-R)-k^R 



(6) 



where R is the relative amoxmt of boimd analyte, c is ana- 
lyte concentration, t is time, and and are the associa- 
tion and dissociation kinetic rate constants, respectively 
[21]. This interaction model assumes 1:1 binding, rapid 
mixing of the analyte from the bulk phase to the sensor 
surface layer, and single-step binding. Observed binding, 
however, may deviate from this simple model due to more 
complex mechanisms of interaction and mass transport 
limitations [22]. Equation (6) yields for R: 



R(t) = 



(l_g-(V+Ad)/)+i?Q 



(7) 



where R^ denotes the initial amount of analyte bound at 
the time /^O [21]. 



Analyte 





Sefisof surface coated with antibodies 

Fig. 6 Direct detection 
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Fig. 7 Direct detection. Binding between antibody and analyte cal- 
culated for four different concentrations of analyte, k^3^1Qf^ mol"* 
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Various measurement formats have been adopted in 
SPR biosensing to ensure that the monitored binding event 
produces a measurable sensor response. The most fre- 
quently used measurement formats are direct detection, 
sandwich assay, and inhibition assay. In direct detection 
format, analyte in a sample interacts with a biomolecular 
recognition element (antibody) immobilized on the sensor 
surface. Fig. 6. The resulting refractive index change is di- 
rcctly proportional to the concentration of analyte. 

Figure? which shows a kinetic model of the interac- 
tion between antibody and analyte suggests that the bind- 
ing between the target analyte and antibody is fast ini- 
tially. As the interaction progresses, the binding rate grad- 
ually decreases and eventually reaches a state in which the 
association and dissociation processes are in equilibrium. 
The time required for the interaction to reach the equilib- 
rium depends on the concentration of analyte and is longer 
for lower concentrations of analyte. 

Figure 8a shows dependence of the relative binding at 
equilibrium on the concentrations of analyte. At low ana- 
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Fig. 8 Direct detection: (top) relative equilibrium binding as a 
function of analyte concentration; (bottom) initial binding rate as 
a function of analyte concentration, A;a=3xl0^moH Ls~^; k^=5x 
lO^s-i 



lyte concentrations the equilibrium binding increases lin- 
early with analyte concentration. At higher analyte con- 
centrations, the binding sites provided by the biomolecu- 
lar recognition elements are saturated and a further in- 
crease in the analyte concentration produces a smaller in- 
crease in the amount of bound analyte. The initial binding 
rate dR/dt (t=0) is directly proportional to the association 
rate constant and analyte concentration (Fig. 8b). Both the 
amount of analyte bound at eqmlibrium and initial bind- 
ing rate can be used to determine analyte concentration. 
The measurement of the bmding rate is faster and offers a 
larger dynamic range than measurement of the equilib- 
rium binding. In the sandwich assay format the measure- 
ment consists of two steps. In the first step, sample con- 
taining analyte is brought in contact with the sensor and 
the analyte molecules bind to the antibodies on the sensor 
surface. Then the sensor surface is incubated with a solu- 
tion containing „secondary" antibodies. The secondary 
antibodies bind to the previously captured analyte further 
increasing the number of bound biomolecules (Fig. 9) and 
thus also the sensor response. 

The inhibition assay is an example of a competitive as- 
say. In this detection format, a sample is initially mixed 
with respective antibodies and then the mixture is brought 
in contact with the sensor surface coated with analyte mol- 
ecules, so that the unoccupied antibodies could bind to the 
analyte molecules (Fig. 10). 

The amoimt of boimd analyte versus time may be esti- 
mated by calculating the equilibrium concentration of an- 
tibody which did not bind to the analyte in the sample and 
then simulating the interaction betwin the unbound anti- 
body and the analyte-derivatized surface. Figure 11 shows 
the equilibrium binding and the initial binding rate as a 
function of analyte concentration, assuming that antibody 
at a concentration of 0.1 nmolL"^ was incubated with sam- 
ple and the mixture was provided with enough time to 
reach equilibrium. The amount of bound antibody and the 
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Fig. 11 Inhibition assay: (a) relative equilibrium binding of the 
antibody as a function of anal3^e concentration; (b) initial binding 
rate as a function of analj^ concentration; A:a=3xl0^mol-^ Ls"^; 
k^=5xl(y^ s"^; antibody concentration 0. 1 nmol L-^ 



initial binding rate are inversely proportional to analyte 
concentration. Fig. 11. 



Features and challenges 

SPR biosensor technology exhibits various advantageous 
features. These include in particular: 

1. Versatility — generic SPR sensor platforms can be tai- 
lored for detection of any analyte, providing a biomo- 
lecular recognition element recognizing the analyte is 
available; analyte does not have to exhibit any special 
properties such as fluorescence or characteristic absorp- 
tion and scattering bands. 

2. No labels required — binding between the biomolecular 
recognition element and analyte can be observed di- 
rectly without the use of radioactive or fluorescent la- 
bels. 



3. Speed of analysis - the binding event can be observed 
in real-time providing potentially rapid response. 

4. Flexibility — SPR sensors can perform continuous mon- 
itoring as well as one-time analyses. 

SPR biosensors exhibit two inherent limitations: 

1. Specificity of detection — specificity is solely based on 
the ability of biomolecular recognition elements to rec- 
ognize and capture analyte. Biomolecular recognition 
elements may exhibit cross-sensitivity to structurally 
similar but non-target molecules. If the non-target mol- 
ecules are present in a sample in a high concentration, 
sensor response due to the non-target analyte molecules 
may conceal specific response produced by low levels 
of target analyte. 

2. Sensitivity to interfering effects - similar to other affin- 
ity biosensors relying on measurement of refractive in- 
dex changes, SPR biosensor measurements can be 
compromised by interfering effects which produce re- 
fractive index variations. These include non-specific 
interaction between the sensor surface and sample (ad- 
sorption of non-target molecules by the sensor surface), 
and background refractive index variations (due to sam- 
ple temperature and composition fluctuations). 



Advances in SPR biosensor teciinoiogy 

This paper follows up an SPR sensor technology review 
paper published in Sensors and Actuators B in 1999 [9], 
and therefore focuses primarily on recent advances in SPR 
biosensor technology. In this section the following areas 
are discussed in particular: SPR sensor platforms, data 
analysis, and biomolecular recognition elements. 



Instrumentation 

In recent years, research into optical platforms for SPR 
biosensors has been increasingly application-oriented, tar- 
geting specific appHcation areas and providing solutions 
meeting imique requirements of specific applications. Two 
important representatives of this trend are: development 
of laboratory SPR sensor platforms with a large number 
of sensing channels for high-throughput screenmg appli- 
cations and development of mobile SPR sensor platforms 
for analysis of complex samples in the field. 

In traditional multichaimel SPR sensors, SPWs were 
excited via a prism coupler in multiple areas which were 
arranged perpendicularly to the direction of propagation 
of SPWs; angular [23] or spectral [24] distribution of re- 
flected light was analyzed to yield information about the 
measurand in each channel. While this spectroscopic ap- 
proach led to development of high-performance SPR sens- 
ing devices [23], the number of sensing channels which 
could be realized using this approach was rather limited 
(<10). In order to increase the number of sensing chan- 
nels, various SPR sensor platforms have been proposed 
[25, 26, 27, 28]. One approach is based on SPR imaging 
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in which a coUimated light beam from a polychromatic 
light source passes through a prism and is made incident 
on an SPR-active metal layer; the reflected light is de- 
tected with a CCD camera after passing through a narrow- 
band interference filter [25]. This approach has been 
demonstrated for monitoring of adsorption of a single- 
stranded DNA-binding protein on to single-stranded DNA 
pattemed into an array of 500x500-micron squares [29], 
The choice of operatiag wavelength for imaging SPR sen- 
sors have been discussed by Johansen et al. [30]. An alter- 
native approach [26] is based on detection of spatial 
changes in the phase of light exciting an SPW and inter- 
ferometry. Two interferometric schemes have been pro- 
posed. In the Mach— Zehnder interferometer-based scheme 
monochromatic light was split into reference and signal 
beams; the signal beam passed through a prism and, after 
reflection from an SPR-active metal layer, was recom- 
bined with the reference beam producing an interference 
pattem on a CCD camera [28] . In the TEr-TM polarization 
interferometer TE and TM polarized beams passed through 
a prism and, after reflection from an SPR-active surface, 
were shifted with respect to each other and recombined by 
means of a polarizer producing an interference pattem on 
a CCD camera [26]. Another interesting approach is based 
on SPR microscopy [31] and uses surface scanning and 
SPWs excited by means of an objective lens [27]. Most re- 
cently a new approach has been proposed which is based 
on spectroscopy of SPW on an array of diffraction grat- 
ings [32]. 

A great deal of research has been focused on develop- 
ment of mobile SPR sensor platforms with referencing ca- 
pabilities enabling applications of SPR biosensors in out- 
of-laboratory enviroimients and for analysis of complex 
samples. The traditional approach consisted in simultane- 
ous SPR measurements in two sensing channels containing 
different biomolecular recognition element coatings, one 
with (signal channel) and one without (reference channel) 
affinity for the analyte, and subtraction of the reference 
channel response from that of the signal channel [24] . Re- 
cently, a new approach to multichannel SPR sensing has 
been developed which is based on excitation of surface 
plasmons in different sensing channels at different wave- 
lengths and encoding information from different sensing 
channels into different regions of the light spectrum. This 
can be accomplished in a sensing element of special de- 
sign in which light is made incident on the SPR-active 
metal at different angles of incidence (Fig. 12a) [33], or by 
employing a thin dielectric overlayer which shifts the res- 
onant wavelength for a part of the sensing surface to 
longer wavelengths (Fig. 12b) [34]. This approach has been 
demonstrated to have capacity to discriminate effects oc- 
curring in the proximity of the sensor surface (specific 
binding, non-specific adsorption) from those occurring m 
the whole medium (interfering background refractive in- 
dex fluctuations) which is a prerequisite for advanced ref- 
erencing [35] . Referencing approaches have been studied 
[36, 37]. It was found that a residual error for compensa- 
tion of interfering background refractive index variations 
is typically 1—3% of the total refractive index change; 
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Fig. 12 Dujil-surface-plasmon spectroscopy, (a) Geometry with 
two different singles of incidence [33]. (b) Geometry with a high 
refractive index overlayer [34] 



compensation for temperature variations is less accurate 
with an error of 5—10% of the total response due to the 
temperature change [37]. 

There is continuing effort to develop miniature SPR 
platforms based on miniaturized prism couplers [38] and 
optical fibers. The first fiber optic SPR sensors were based 
on wavelength-modulation in multimode optical fibers with 
partly removed cladding and a metal film deposited sym- 
metrically around the exposed section of fiber core [39] or 
on intensity-modulation in single-mode optical fibers which 
were side-polished and coated with a thin metal film [40]. 
These SPR sensors suffered from modal noise (multimode 
fiber-based sensors) and polarization instability (single- 
mode fiber-based sensors). Two methods of overcoming 
the limitation originating from the need to precisely con- 
trol polarization of light in the single-mode optical fiber- 
based SPR sensors have been developed. In the first method, 
Ught from a polychromatic light source passes through a 
Lyott fiber optic depolarizer which produces unpolarized 
light which is then coupled into a fiber optic SPR sensing 
element. The transmitted light is analyzed with a spectro- 
graph [41], The second method uses a polarization-main- 
taining fiber to control polarization of light in the fiber- 
optic SPR sensing element. Fig. 13. This approach pro- 
vides best suppression of polarization effects due to fiber 
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Fig. 13 SPR sensor based on a 
polarization-maintaining sin- 
gle-mode optical fiber 
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deformations resulting in enhanced-stability fiber optic SPR 
sensors with performance comparable to bulk-optic SPR 
sensors [42] . 

Research is also carried out to improve detection capa- 
bilities of SPR biosensors by exploiting a special type of 
an SPW, the so-called long-range surface plasmon. The 
use of long-range surface plasmons provides two benefits 
- increased sensitivity and very narrow angular or spectral 
dips, which makes it possible to determine the spectral po- 
sition of the SPR dip with a high accuracy [43], As long- 
range surface plasma waves penetrate deeper into the 
probed medium, their use benefits especially biosensors 
with extended matrices of biomolecular recognition ele- 
ments [44], 

Sensor data analysis 

Detection limits of SPR biosensors are ultimately con- 
strained by the noise-based precision of the SPR instru- 
ment itself. The precision depends on noise contributions 
made by individual components of an SPR sensor and 
data processing method. Therefore, study of noise in SPR 
sensors and development of optimized algorithms for pro- 
cessing data from SPR sensors have received much atten- 
tion lately. Earlier approaches to understanding the per- 
formance of SPR sensor data analysis have included stud- 
ies of the effects of noise and comparisons between algo- 
rithms. Locally-weighted parametric regression and other 
methods were compared for low- and high-noise detec- 
tors, demonstrating that a small number of low-noise de- 
tector pixels may ou^erform a larger but noisier detector 
array [45]. Linearization of data processing algorithms 
and an optimal linear data analysis method were proposed 
as a means of optimizing algorithm parameters [46]. Con- 
tribution of analog-to-digital converter resolution and num- 
ber of pixels in the detector array to the performance of a 
number of SPR data-analysis algorithms was studied [47]. 
Also, effects of sensitivity deviations on concentration 
analyses and kinetic studies have been investigated [48] . 
Most recently, sources of noise were investigated for a 
wavelength-modulated SPR sensor. Shot noise of the de- 
tector was found to be the dominant source of noise and 
an analytical formula was derived which allows predic- 
tion of the noise of the sensor output based on the detec- 
tor noise [49]. 



Biomolecular recognition elements 
and their immobilization 

The main types of biomolecular recognition elements used 
in SPR biosensors include antibodies, nucleic acids and 
biomimetic materials. Antibodies are used most frequently 
because of their high affinity, versatility, and commercial 
availability. Various immobilization chemistries have been 
developed to attach antibodies to SPR-active gold films. 
Traditional approaches include formation of streptavidin 
layer on the gold surface followed by attachment of bi- 
otinylated antibodies [50], use of self -assembled alka- 
nethiol films with suitable reactive groups [51], use of a 
hydrogel matrix composed of carboxyl-methylated dex- 
tran chains which can be modified allowing antibodies to 
be attached via surface-exposed amine, carboxyl, sulfhydryl, 
and aldehyde groups [52]. Altematively, SPR sensing sur- 
faces may be functionalized by thin polymer films to 
which antibodies are coupled via amino groups [53]. Pro- 
tein contact printing has been examined for spatially-con- 
trolled attachment of bovine serum albumin (BSA) and 
dinitrophenylated BSA onto adjacent reference and signal 
channels of a dual-channel SPR sensor [54]. DNA can be 
immobilized on gold SPR sensor surfaces by formation of 
a streptavidin layer on a gold surface followed by attach- 
ment of biotinylated DNA [55] . A multistep surface mod- 
ification based on alkanethiol self -assembled monolayers 
has also been used to attach DNA to gold surfaces [56]. 
Recentiy, biomimetic materials consisting in molecularly 
imprinted polymers (MIPs) have been exploited in SPR 
biosensors [57]. 



Applications of SPR biosensors 

Two major application areas for SPR biosensing are in de- 
tection and identification of biological analytes and bio- 
physical analysis of biomolecular interactions. This re- 
view focuses on applications for detection and identifica- 
tion of biological analytes; recent advances in SPR-based 
biomolecular interaction analysis can be foimd in Ref. 
[58]. 

Numerous SPR biosensors have been developed for 
detection and identification of specific analytes. These 
biosensors use a number of platform designs, biomolecu- 
lar recognition elements and detection formats. The choice 
of detection format for a particular application depends on 
size of target anal3rte molecules, binding characteristics of 
biomolecular recognition element, and the range of ana- 
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lyte concentrations to be measured. Direct detection is 
usually preferred in applications where direct binding of 
analyte of concentrations of interest produces a sufficient 
response. If necessary, the lowest detection limits of the 
direct SPR biosensors can be improved by using sandwich 
assay. The secondary antibodies may also be coupled to 
large particles such as latex particles [59] and gold beads 
[60] to further enhance the SPR sensor response. Smaller 
analytes (molecular weight <1000) are usually measured 
using inhibition assay. 



Detection of small analytes 

SPR biosensors have been demonstrated for small ana- 
lytes relevant to enviroimiental protection (simazine and 

atrazine), medicine (morphine, methamphetamine, and 
theophyline), and food safety (fumonisin Bl, sulfameth- 
azine, and sulfadiazine). Minunni and Mascini used an 
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SPR sensor and binding inhibition assay to detect atrazine 
[61]. Monoclonal antibodies against atrazine were mixed 
with a sample containing atrazine and free antibody con- 
centration was determined by exposing the sample to an 
atrazine derivative-coated SPR biosensor. A detection limit 
of 0.05 ngmL"^ was achieved [61]. Figures 14 and 15 show 
typical sensorgrams for binding inhibition assay detection 
of atrazine using a wavelength-modulated SPR sensor; the 
antibody concentration was 4|igmL-i. 

Detection of the triazine herbicide simazine in water 
samples was demonstrated using an integrated optical SPR 
sensor and binding inhibition assay involving anti-simazine 
IgG antibodies or anti-simazine Fab fragments [62, 63]. 
The lowest detection limit for simazine was determined 
to be 0.16ngmL-i for anti-simazine IgG antibodies and 
O.llngmL"^ for anti-simazine Fab fragments [63]. Mor- 
phine detection based on an SPR sensor and binding inhi- 
bition assay was reported by Miura et al. [64], who de- 
tected morphine at concentrations down to O.lngmL"^. 
Sakai et al. developed an SPR biosensor for detection of 
methamphetamine using binding inhibition assay; the low- 
est detection limit was determined to be 0.1 ng mL-^ [65]. 
Direct detection of fvunonisin B 1 in aqueous samples was 
demonstrated by Mullet et al. using a laboratory SPR sys- 
tem with angular modulation and polyclonal antibodies 
[66]. The detection limit was determined to be 50ngmL-i. 
Sulfamethazine in milk was detected using a commercial 
SPR system based on angular modulation (Biacore) and 
binding inhibition assay involving anti-sulfamethazine 
polyconal antibodies. The lowest detection limit for sul- 
famethazine was determined to be 1 ng mL-i [67] and 2 ng 
mL~^ [68]. Baxter et al. demonstrated SPR biosensor- 
based detection of streptomycin in milk using a commer- 
cial SPR biosensor (Biacore 2000) and inhibition assay 
involving polyclonal anti-streptomycin antibodies. The 
lowest detection limit of the SPR biosensor for strepto- 
mycin was determined to be 4ngniL"^ [69]. Elliott et al. 
demonstrated detection of sulfadiazine in pig bile using a 
commercial SPR biosensor (Biacore) and inhibition assay 
with a detection limit of 20 ng inL~^ [70] . A laboratory 
SPR biosensor with angular modulation was combined 
with a molecularly imprinted pol3niier for detection of 
theophylline by Lai et al. For aqueous samples the lowest 
detection limit was estimated at 0.4mgnLL"^ [57]. 



^ Detection of medium-size analytes 
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Fig. 15 Detection of atrazine using inhibition assay. Equilibrium 
response and initial binding rate as a function of atrazine concen- 
tration 



Examples of medium-size analytes detected by SPR 
biosensor technology include food safety-related analytes 
such as staphylococcal enterotoxin B, botulinimi toxin, 
and E. coli enterotoxin. Choi et al. demonstrated a direct 
SPR biosensor for botulinum toxin using a commercial 
SPR biosensor (Biacore X) and monoclonal antibodies; 
this biosensor was able to detect botulinvun toxin in buffer 
at a concentration of 2.5 |LLgmL"n71]. Spangler et al. 
demonstrated direct detection of E. coli enterotoxin in 
aqueous solutions; the lowest detection limit was estab- 
lished to be 6 (LigmL"^ [72]. Detection of staphylococcal 
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enterotoxin B was performed using a commercial angular 
sensor (Biacore 3000) [73] and laboratory wavelength 
modulation-based SPR sensors [74] . The lowest detection 
limits were 1-10 ng mLr^ (in milk and meat) [73], 5 ng mLr^ 
(direct detection in buffer), and 0.5 ng mL-^ (sandwich as- 
say in buffer and miUc) [74], A typical sensorgram for de- 
tection of Staphylococcal enterotoxins B (SEB) is shown 
in Fig. 16 [74], Figure 17a shows the equilibrium sensor 
response for both the direct capture of SEB and amplifi- 
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Fig. 18 Direct detection of Salmonella enteritidis. Kinetic re- 
sponse to Salmonella at a concentration of lO^cfumLr^ 



cation by secondary antibodies as a function of SEB con- 
centration. The use of secondary antibodies provides 10-fold 
increase in the sensor response. In Fig. 17b, initial binding 
rates for SEB and secondary antibody concentration are 
shown for several SEB concentrations [74], 



Fig. 16 Detection of Staphylococcal enterotoxin B using sand- 
wich assay. Kinetic response [74] 
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Fig. 17 Detection of Staphylococcal ent^otoxin B. (a) Equilib- 
rium response for different concentrations of SEB. (b) Initial bind- 
ing rate for different concentrations of SEB [74] 



Detection of large analytes 

Representatives of large analytes targeted by SPR biosen- 
sor technology include bacteriaL pathogens such as Esche- 
richia coliy Salmonella enteritidis^ and Listeria monocyto- 
genes. Detection of Escherichia coli 0157:H7 was per- 
formed by Fratamico et al. who used an angular-modula- 
tion commercial SPR sensor and sandwich assay [75]. 
They used monoclonal antibodies immobilized on the 
sensor surface for capturing E, coli and polyclonal second- 
ary antibodies for enhancing the specific sensor response. 
The lowest detection limit for E. coli was established at 
SxlO^cfumL-^ Direct detection of Salmonella enteritidis 
and Listeria monocytogenes at concentrations down to 
lO^cfiimL-^ was demonstrated by Koubova et al. using a 
laboratory wavelength-modulated SPR sensor and mono- 
clonal antibodies [76]. A sensorgram illustrating direct de- 
tection of Salmonella enteritidis is shown in Fig. 18. 

The main challenges in detecting bacterial analytes at 
lower levels include the low concentration of a particular 
antigen relative to total cellular material and slow diffu- 
sion of bacterial cells to sensor surface [77]. 



Commercializalion of SPR biosensors 

The first commercial SPR biosensor was launched by Bia- 
core Intemational AB in 1990. In the following decade, 
Biacore has developed a range of laboratory SPR instru- 
ments (Biacore 1000, Biacore 2000, Biacore 3000, Bia- 
core C, Biacore X, Biacore J, Biacore Q) [78]. Most re- 
centiy the Biacore S51 has been developed which offers 
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higher sensitivity and throughput. Other SPR sensors have 
been developed by British Windsor Scientific (IBIS) [79], 
Nippon Laser and Electronics Laboratory (SPR-670 and 
SPR-CELLIA) [80], Texas Instruments (Spreeta) [81], and 
Analytical |Li-Systems (BIO-SUPLAR 2) [82], 



Oiitlook 

Over the past ten years, surface plasmon resonance (SPR) 
biosensor technology has made great strides, and a large 
number of SPR sensor platforms, biomolecular recogni- 
tion elements, and measurement formats have been devel- 
oped. There has been growing interest in commercializa- 
tion of SPR biosensor technology leading to a number of 
systems available on today's market. SPR biosensors have 
played a significant role in research into biomolecules and 
their interactions and have been increasingly used for de- 
tection and identification of chemical and biological sub- 
stances. SPR biosensors have been particularly successful 
for detection of small and medium size analytes, a number 
of which have been detected at practically relevant levels. 
Detection limits for large analytes such as bacteria and 
viruses still need to be improved to meet today's needs. 

Undoubtedly, future development of SPR biosensors will 
be driven by the needs of the consumer. SPR biosensor 
technology has the potential to benefit numerous impor- 
tant fields including pharmaceutical research, medical di- 
agnostics, environmental monitoring, food safety, and se- 
curity. Applications in these areas present unique chal- 
lenges and impose special requirement on analytical tech- 
nologies. Pharmaceutical research, which was fast to adopt 
optical biosensors, needs laboratory-based high-through- 
put systems with high sensitivity to facilitate parallel 
screening. Medicine would benefit from high-throughput 
diagnostic tools for centralized laboratories and from com- 
pact diagnostic systems dedicated to selected diagnostic 
applications which could be used at clinics. In addition, 
there is a growiag interest in tools for home medical (self-) 
diagnostics. Mobile analytical systems enabling rapid de- 
tection of food-borne pathogens in food would be impor- 
tant for food producers, processors, distributors and regu- 
latory agencies and thus benefit the food safety. Environ- 
mental monitoring would benefit from detection systems 
which could be deployed in the field for continuous mon- 
itoring and from mobile systems enabling fast identifica- 
tion of environmental threat. SPR biosensors could also 
play an important role in defense, where fast, portable and 
rugged units are needed for early detection and identifica- 
tion of biological warfare agents in the field. Develop- 
ment of these systems will require significant advances in 
miniaturization of SPR biosensing platforms, develop- 
ment of robust biomolecular recognition elements witii 
high specificity and long storage life, integration of SPR 
sensor platforms with microfluidic devices, and applica- 
tion-specific sampling systems. 

Acknowledgement This work was supported by the Grant Agency 
of the Czech Republic under contracts 102/03/0633, 303/03/0249, 
and 203/02/1326. 



Iteforencss 

1. Ghindilis AL, Atanasov P, W^iUdns M, Wilkins E (1998) 
Biosens Bioelectron 13:11 3—13 1 

2. Chu X, Un ZH, Shen GL, Yu RQM (1995) Analyst 120: 
2829^2832 

3. Gauglitz G (1996) Opto-chemical and opto-immuno sensors, 
sensor update, vol 1. VCH, Weinheim 

4. Rowe-Taitt CA, Hazzard JW, Hoffinan KE, Cras JJ, Golden 
JP, Ugler FS (2000) Biosens Bioelectron 15:579-589 

5. Piehler J, Biecht A, GaugHtz G (1996) Anal Chem 68:139-143 

6. Heidonan RG, Kooyman RPH, Greve J (1993) Sens Actuators 
B 10:209-217 

7. Clerc D, Lukosz W (1994) Sens Actuators B 19:581-586 

8. Cush R, Cronin JM, Stewart WJ, Maule CH, Molloy J, God- 
dard NJ (1993) Biosens Bioelectron 8:347-353 

9. Homola J, Yee S, Gauglitz G (1999) Sens Actuators B 54:3-15 

10. Homola J, Yee S, Myszka D (2002) Surface plasmon biosen- 
sors. In: Ligler FS, Taitt CR (eds) Optical biosensors: present 
and future. Elsevier 

11. Rabbany SY, Lane WJ, Marganski WA, Kusterbeck AW, 
Ligler FS (2000) J Immunol Methods 246:69-77 

12. Boardman AD (1982) (ed) Electromagnetic surface modes. 
Wiley, Chichester 

13. Reather H (1983) Surface plasmons on smooth and rough sur- 
faces and on gratings. Springer tracts in modem physics, vol 
111. Springer, BerHn Heidelberg New York 

14. Snyder AW, Love JD (1983) Optical waveguide theory. Chap- 
man and Hgdl, London 

15. Parriaux O, Voirin G (1990) Sens Actuators A 21-23:1 137 

16. Hutley MC (1982) Diffration gratings. Academic Press, Lon- 
don 

17. Kooyman RPH, Kolkman H, van Gent J, Greve J (1988) Anal 
Chim Acta 213:35-45 

18. Yeatman EM (1996) Biosens Bioelectron 11:635-649 

19. Homola J (1997) Sens Actuators B 41:207-211 

20. Homola J, Koudela I, Yee S (1999) Sens Actuators B 54:16-24 

21. Edwards PR, Leatherbarrow RJ (1997) Ansd Biochem 246:1-6 

22. Vijayendran RA, Ligler FS, Leckband DE (1999) Anal Chem 
71:5405-5412 

23. LofSs S, Malmqvist M, Ronnb^g I, StenbCTg E, Liedberg B, 
Limdstrom I (1991) Sens Actuators B 5:79 

24. Nenninger GG, Clendenning JB, Furlong CE, Yee S (1998) 
Sens Actuators B51:38 

25. Jordan CE, Frutos AG, Thiel AJ, Com RM (1997) Anal Chem 
69:4939 

26. Nikitin PI, Beloglazov AA, Kabashin AV, Valeiko MV, 
Kochergin V E (1999) Sens Actuators B54:43 

27. Kano H, Knoll W (2000) Optics Communications 182:11-15 

28. Kabashin A V, Nikitin P I (1998) Opt Commun 150:5-8 

29. Nelson BP, Fmtos AG, Brockman JM, Com RM (1999) Anal 
Chem 71:3928-3934 

30. Johnasen K, Arwin H, Lundstrom I, Liedberg B (2000) Rev Sci 
Instrum 71:3530-3538 

31. Rothenhausler B, Knoll W (1988) Nature 332:688 

32. Dostalek J, Homola J, Miler M (2002) EUROPT(R)ODE VI, 
Manchester, UK, Book of Abstracts, 265 

33. Homola J, Dostalek J, Ctyroky J (2001) Proc SPIE 4416:86 

34. Homola J, Lu HB, Yee S (1999) Electr Lett 35:1105 

35. Homola J, Lu HB, Nenninger GG, Dostalek J, Yee S (2001) 
Sens Actuators B76:403-410 

36. Ober RJ, Ward E S (1999) Anal Biochem 271:70-80 

37. Homola J, Dostalek J, Piliarik M, Yee S (2002) EUROP- 
T(R)ODE VI, Manchester, UK, Book of Abstracts, 71 

38. Stemmler I, Brecht A, Gauglitz G (1999) Sens Actuators B 
54:98-105 

39. Jorgenson RC, Yee S (1993) Sens Actuators B12:213 

40. Homola J (1995) Sens Actuators B 29:401^05 

41. Homola J, Piliarik M, Slavlk R, Ctyroky J (2001) SPIE Proc 
4416:82-85 

42. Piliarik M, Homola J, Manikova Z, Ctyroky J (2003) Sens Ac- 
tuators B 90:236-242 



539 



43. Nenninger GG, Tobiska P, Homola J, Yee S (2001) Sens Actu- 
ators B74: 145 

44. Liedberg B, Lundstrom I, Stenberg E (1993) Sens Actuators 
Bll:63 

45. Johnston KS, Booksh KS, Chinowsky TM, Yee S (1999) Sens 
Actuators B 54:80-88 

46. Chinowsky TM, Jung LS, Yee S (1999) Sens Actuators B 54: 
89^97 

47. Johansen K, Stdlberg R, Lunstr^m I, Liedberg B (2000) Meas 
SciTechnol 11:1630-1638 

48. Johansen K, Lundstr3m I, Liedberg B (2000) Biosens Bioelec- 
tron 15:503-509 

49. Nenninger GG, Piliarik M, Homola J (2002) Meas Sci Technol 
13:2038-2046 

50. Morgan H, Taylor DM, (1992) Biosens Bioelectron 7:405^10 
51.1>uschl C, Sevin-Landais A, Vogel H (1995) Biophys J 70: 

1985 

52. Lofas S, Malmqvist M, Ronnberg I, Stenberg E, Liedberg B, 
Lundstrom I (1991) Sens Actuators B 5:79 

53. Nakamura R, Muguruma H, Ikebukuro K, Saseiki S, Nagata R, 
Karube I, Pedersen H (1997) Anal Chem 69:4649 

54. Lu HB, Homola J, CampbeU CT, Nenninger GG, Yee S, Rat- 
ner BD (2001) Sens Actuators B 74:91-99 

55. Watts HJ, Yeung D, Parkes H (1995) Anal Chem 67:4283- 
4289 

56. Brockman JM, Frutos AG, Com RM (1999) J Am Chem Soc 
121:8044-8051 

57. Lai EPC, Fafara A, VanderNoot VA, Kono M, Polsky B (1998) 
Can J Chem 76:265-273 

58. Rich RL, Myszka DG (2002) J Mol Recognit 15:352^376 

59. Severs AH, Schasfoort RBM (1993) Biosens Bioelectron 8:365 

60. Leung PT, Pollard-Knight D, Malan GP, Finlan MF (1994) 
Sens Actuators B22:175 

61. Minunni M, Mascini M (1993) Anal Lett 26:1441 

62. Mouvet C, Harris RD, Maciag C, Luff BJ, Wilkinson JS, 
Piehler J, Brecht A, GaugUtz G, Abuknesha R, Ismail G (1997) 
AdsH Chim Acta 338:109 



63. Harris RD, Luff BJ, Wilkinson JS, Piehler J, Brecht A, 
Gauglitz G, Abuknesha RA (1999) Biosens Bioelectron 14:377 

64. Miura N, Ogata K, Sakai G, Uda T, Yamazoe N (1997) Chem 
Lett 8:713 

65. Sakai G, Nakata S, Uda T, Miura N, Yamazoe N (1999) Electro- 
chim Acta 44:3849 

66. MuUett W, Edward PC, Yeung MJ (1998) Anal Biochem 258: 
161-167 

67.Stemesjo A, Mellgren C, Bjorck L (1996) ACS Synq) Ser 

621:463^70 

68. Gaudin V, Pavy M-L (1999) JAOAC Int 82:1316-1320 

69. Baxter GA, Ferguson IP, O'Conner MC, EUiott CT (2001) 
J Agric Food Chem 49:3204-3207 

70. Elliot CT, Baxter GA, Crooks SRH, McCaughey WJ (1999) 
Food Agric Immunol 11:19—28 

71. Choi K, Seo W, Cha S, Choi J (1998) J Biochem Mol Biol 
31:101-105 

72. Spangler BD, Wilkinson EA, Murphy JT, Tyler BJ (2001) Anal 
Chim Acta 444:149-161 

73. Rasooly A (2001) J Food Prot 64:37-43 

74. Homola J, Dostalek J, Chen S, Rasooly A, Jiang S, Yee S 
(2002) Int J Food Microbiol 75:61-69 

75. Fratamico PM, Strobaugh TP, Medina MB, Gehring AG (1998) 
Biotechnol Tech 12:571-576 

76. Koubova V, Brynda E, Karasova L, Skvor J, Homola J, Dos- 
talek J, Tobiska P, Rosicky J (2001) Sens Actuators B 74: 
100-105 

77. Perkins EA, Squirrell DJ (200O) Biosens Bioelectron 14:853 

78. Biacore websitewww.biacore.com 

79. IBIS Technologies website www.ibis-spr.nl 

80. Nippon Laser and Electronics Laboratory website www.nle-lab. 
co.jp/EnglishyZO-HOME.htm 

81. Texas Instruments website www.ti.com/sc/docs/productsAnsp/ 
control/spreeta 

82. Anal3^cal |l-Systems website www.micro-systems.de 



